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19  ABSTRACT  (Continue  on  reverse  if  necessary  eiod  identify  by  block  number) 

—The  purpose  of  this  research  is  to  develop  a  stochastic  model  of  distributed  radar  targets,  especially 
ships,  that  directly  incorporate  target  structure  and  motion.  The  model  was  required  to  be  especially  useful  as 
a  tool  in  the  stochastic  simulation  and  analysis  of  tracking-radar  signals  over  time  intervals  that  are  short  rela¬ 
tive  to  the  time  constants  of  the  target  motion.  The  model  is  based  on  the  observation  that  distributed  targets 
often  appear  to  radars  as  being  composed  of  several  dominant  scatterers.  A  concept  (unit-scatterer)  is  intro¬ 
duced  that  quantifies  this  observation  and  that  leads  to  a  useful  model  of  distributed  targets.  Based  on  this 
concept  and  assuming  the  presence  of  over-water  multipath,  analytical  representations  of  radar  cross  section 
and  glint  are  developed,  and  implications  of  the  small  time-interval  requirement  are  investigated.  Using  these 
representations,  a  simulation  is  developed  and  used  to  investigate  the  stochastic  properties  of  both  radar  cross 
section  and  glint  for  an  example  ship  target.  Simulation  outputs  are  presented  and  analyzed  to  illustrate  the 
implications  of  the  model  given  variations  in  the  significant  parameters. 
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*>-The  model  was  developed  to  incorporate  the  major  strengths  of  the  existing  deterministic  and  stochastic 
models:  ability  to  account  directly  for  target  structure  and  motion,  and  ease  of  obtaining  target-signature  time 
series,  respectively.  The  deterministic  models,  though  precise,  are  extremely  inefficient  in  generating  these 
time  series  and  require  a  great  deal  of  information  about  target  structure  and  motion.  The  stochastic  models 
do  not  directly  account  for  target  structure  and  motion  and  rely  heavily  on  target  measurements  .^The  model 
presented  here  efficiently  generates  target-signature  time  series  given  information  about  target  structure  and 
motion.  The  existing  deterministic  and  stochastic  models  can  be  viewed  as  limiting  cases  of  this/new  model 
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MODELING  OF  ELECTROMAGNETIC 
SCATTERING  FROM  SHIPS 


1.  INTRODUCTION 

Scattering  of  electromagnetic  (EM)  energy  has  been  a  subject  of  interest  to  the  scientific  com¬ 
munity  since  the  nineteenth  century.  The  first  and  most  extensively  studied  scattering  object  has  been 
the  sphere.  However,  very  few  scattering  problems  have  solutions  that  are  known  and  exact  although 
many  problems  have  asymptotic  solutions  (lj.  The  problems  that  have  been  solved  have  been  prob¬ 
lems  of  scattering  from  simple  geometric  shapes  [1,2].  The  modeling  of  complex  targets  is  an 
extremely  difficult,  inexact  task  [3] ;  quoting  from  Ref.  4,  "An  exact  solution  to  the  scattering  problem 
for  a  complex  target  is  out  of  the  question.  .." 

Most  of  the  interest  in  EM  scattering  has  come  from  the  field  of  radar.  Radars  obtain  information 
about  targets  by  illuminating  the  target  and  then  measuring  the  EM  field  scattered  by  the  target  and 
environment.  Radars  have  been  used  for  this  purpose  since  World  War  II.  Since  then  their  perfor¬ 
mance  capabilities  have  improved  rapidly.  Paralleling  the  improved  radar  performance  has  been  an 
increasing  need  for  better  understanding  of  target  scattering  properties.  This  need  has  directed  most  of 
the  scattering  research  since  World  War  II.  Almost  all  targets  of  interest  for  military  and  civilian  radars 
have  physical  dimensions  much  larger  than  the  radar  wave  length  (scattering  in  the  optical  region); 
that  is,  they  are  distributed  targets.  Clearly,  any  scattering  model  that  is  to  be  useful  in  radar  applica¬ 
tions  must  be  limited  in  its  ability  to  fully  describe  target  scattering. 

Although  EM  scattering  is  the  phenomenon  which  makes  radars  possible,  by  producing  scattered 
fields  for  the  radar  to  sense,  it  is  well  known  that  it  is  only  necessary  to  understand  two  projections  of 
these  fields  to  predict  the  performance  of  most  radar  systems.  These  projections  are  radar  cross  section 
(RCS)  and  glint.  RCS  represents  the  apparent  size  of  the  target,  and  glint  represents  the  apparent  loca¬ 
tion,  in  angle,  of  the  target,  as  sensed  by  the  radar.  Though  generally  analyzed  as  separate  phenomena, 
it  is  well  known  that  RCS  and  glint  are  different  manifestations  of  the  same  target-induced  effects. 
Most  of  the  research  on  scattering  since  the  development  of  radar  has  been  directed  to  understanding 
RCS  and  glint  rather  than  scattering  per  se. 

Distributed  targets  are  typically  modeled  as  a  finite  number  of  individual,  point-source  scatterers 
whose  characteristics  are  determined  by  the  target's  structure.  (We  refer  to  such  models  as  N-source 
models.)  Both  RCS  and  glint  can  be  represented  by  using  such  models.  The  various  N-source  models 
that  have  arisen  differ  depending  on  how  the  individual  scatterers  are  characterized.  Generally,  the 
scatterer  characterizations  that  are  made  depend  on  the  intended  application  of  the  target  model.  Build¬ 
ing  on  these  earlier  N- source  models,  this  report  develops  a  new  approach  to  modeling  the  scattering 
from  ships. 

1.1  The  Scattering  Problem 

Practical  solutions  to  the  modeling  of  distributed  targets  are  approximations  that  range  from 
extremely  complex  models  that  describe  many  aspects  of  scattering  (and  are  difficult  and/or  expensive 
to  use)  to  simple  models  that  ignore  many  aspects  (but  are  easy  and/or  inexpensive  to  use).  Each 
model  is  developed  subject  to  constraints  imposed  by  application  requirements.  We  remark  that  no 
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model  can  completely  describe  all  aspects  of  a  physical  phenomenon  and  very  "crude"  models  are  often 
the  most  useful;  modeling  a  resistive  element  in  an  electric  circuit  by  E— IR  is  often  quite  useful 
though  such  a  model  ignores  known  resistive  element  dependencies  on  temperature,  power,  frequency, 
etc.,  and  hence  is  not  useful  in  predicting  failure  modes,  noise  properties,  etc.,  of  the  element. 

/.  1. 1  Scattering  Models 

The  models  that  have  been  developed  to  represent  complex  target  scattering  can  be  described  as 
being  one  of  three  types  (or  a  combination  thereof):  deterministic,  statistical,  or  stochastic.  The  sim¬ 
plest  deterministic  model  is  one  where  the  target  is  represented  as  a  fixed-position,  point-source  radia¬ 
tor  with  constant  RCS  and  no  glint.  Deterministic  scattering  solutions  for  simply  shaped  objects  can 
often  be  obtained  by  using  the  geometric  theory  of  defraction  or  physical  optics  theory.  The  most  com¬ 
plex  deterministic  models  represent  the  target  as  a  collection  of  many,  simply  shaped  objects  for  which 
scattering  solutions  exist.  The  scattered  field  of  the  target  then  becomes  the  superposition  of  the  scat¬ 
tered  fields  of  the  simple  objects,  where  effects  such  as  shadowing,  multiple  reflections,  etc.,  are 
appropriately  accounted  for. 

The  statistical  models  are  those  in  which  only  the  first-order  statistics  of  one  or  more  aspects  of 
the  scattering,  usually  RCS,  are  modeled.  These  models  are  usually  developed  either  from  measure¬ 
ments  of  the  target  or  by  deriving  the  statistics  from  some  broad  assumption  (e.g.,  that  the  target  is 
composed  of  many  independent  scatterers  of  approximately  equal  scattering  areas)  regarding  the  target 
structure.  The  well-known  Swerling  models  of  RCS  are  examples  of  this  type. 

The  stochastic  models  are  those  that  describe  the  scattered  field  as  a  stochastic  process.  Stochastic 
processes  are  useful  because  they  not  only  specify  first-order  statistical  properties  but  they  also  specify 
higher  order  fluctuation  properties  of  the  target  including,  at  least,  the  second-order  correlation  proper¬ 
ties  (e.g.,  the  correlation  function).  Correlation  properties  of  the  radar’s  receive  signal  must  be  known 
if  an  analysis  of  the  tracking  performance  of  the  radar  is  to  be  made.  The  stochastic  models  used  for 
target  modeling  are  usually  second-order  stochastic  models  (i.e.,  they  specify  the  first-  and  second- 
order  properties  of  the  process)  because  to  date,  radars  exploit  only  the  first-  and  second-order  proper¬ 
ties  of  their  receive  signals  and  because  second-order  processes  are  well  developed  mathematically. 

Stochastic  models  usually  are  either  extensions  or  combinations  of  deterministic  and/or  statistical 
models.  An  example  of  the  extension  approach  is  a  model  that  converts  the  many  scatterer  determinis¬ 
tic  model  described  above  to  a  stochastic  model  by  describing  the  target’s  range  and  orientation  relative 
to  the  radar  as  stochastic  processes  and  then  solving  the  scattering  problem  as  a  function  of  time.  An 
example  of  a  combination  model  approach  is  a  model  that  uses  measurements  that  characterize  not  only 
the  first-order  statistics,  but  also  the  correlation  properties  of  a  target’s  scattering  as  a  function  of 
radar-target  relative  motion.  The  stochastic  model  is  then  a  stochastic  process  that  is  defined  to  pro¬ 
duce  the  measured  statistics  as  a  function  of  the  measured  parameters.  We  refer  to  such  a  model  as  an 
empirical  stochastic  model.  We  note  that  when  the  target  being  modeled  is  extremely  large  (e.g.. 
Earth's  surface  illuminated  by  the  radar,  the  moon,  etc.)  the  approach  of  representing  the  surface  as  a 
stochastic  process  is  often  used. 

1. 1. 2  Fundamental  Constraint 

The  fundamental  constraint  imposed  in  the  development  of  the  model  presented  here  is  that  the 
model  muj,  be  useful  in  the  analysis  and  simulation  of  the  pulse-by-pulse  tracking  performance  of 
pulse  radars.  A  major  implication  of  this  constraint  is  that  the  model  must  be  capable  of  efficiently  pro¬ 
ducing  time-series  representations  of  scattered  radar-signals.  Both  empirical  stochastic  models  and 
deterministic  models  using  many  simple  scatterers  can  be  applied  under  this  constraint.  However,  each 
of  these  models  has  serious  drawbacks  for  the  desired  application.  The  stochastic  model,  although  very 
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efficient  for  simulation  purposes,  is  strongly  dependent  on  actual  measured  data  to  provide  the  needed 
statistical  information,  and  so  extrapolation  of  such  a  model  without  supporting  data  is  questionable. 
Further,  this  type  of  model  is  usually  only  weakly  connected  to  the  physical  process  (target-radar  rela¬ 
tive  motion)  that  is  causing  the  scattered-signal  variations.  The  deterministic  model  has  two  major 
drawbacks.  First,  such  a  model  requires  an  enormous  amount  of  detailed  information  about  the  actual 
target  to  be  modeled;  the  structure  of  that  target’s  surface  must  be  precisely  known.  Second,  the 
simulation  of  such  a  model  usually  is  extremly  time  consuming  because  many  exact  scattering  solutions 
must  be  calculated  to  produce  the  required  time  series. 

In  view  of  the  above,  tracking-radar  studies  and  simulations  would  benefit  from  a  model  that  is 
stochastic,  that  can  be  efficiently  simulated,  and  that  can  be  parametrized  by  the  fundamental  physical 
properties  of  target-radar  relative  motion  and  "significant"  target  structure.  Such  a  model  would  not  be 
as  accurate  as  the  deterministic  model  when  the  target  structure  and  motion  are  precisely  known,  and  it 
might  not  be  as  simple  and  efficient  in  simulation  as  the  empirical  stochastic  model  when  extensive 
experimental  data  is  available.  However,  it  would  incorporate  the  fundamental  simulation  efficiency  of 
a  stochastic  model  and  much  of  the  physical  basis  of  the  deterministic  model.  We  will  refer  to  this  type 
of  model  as  a  phenomenological  stochastic  model. 

1.1.3  Statement  of  the  Problem 

The  objective  of  this  report  is  to  develop  a  phenomenological  stochastic  model  for  the  scattering 
of  radar  signals  from  distributed  targets  using  an  jV-source  formulation.  The  fundamental  parameters 
on  which  the  model  is  based  are  target  motion  and  structure.  We  assume  that  the  target’s  motion  can 
be  characterized  as  a  stochastic  process  and  that  the  target’s  approximate  structure  and  dimensions  are 
known. 

We  further  restrict  the  model  by  imposing  several  additional  constraints  on  its  applications.  The 
model  is  intended  for  use  in  analyzing  the  performance  of  continuous  tracking,  monostatic,  pulse  radars 
that  operate  at  microwave  frequencies  and  track  only  in  range,  azimuth,  and  elevation.  The  targets  of 
interest  are  assumed  to  be  large  with  respect  to  the  radar  wavelength,  structurally  complex,  in  the  radar 
antenna’s  far-field,  and  uniformly  (plane  wave)  illuminated  by  the  radar. 

Because  of  the  importance  of  RCS  and  glint  (and  the  existence  of  associated  data  bases),  the 
model  fidelity  is  to  be  inferred  from  the  validity  of  the  resulting  RCS  and  glint  representations.  Valida¬ 
tion  of  any  model  of  scattering  from  complex  targets  must,  of  practical  necessity,  be  statistical.  That  is, 
the  model  must  produce  a  time  series  of  data  that  can  be  tested  for  statistical  accuracy  relative  to  mea¬ 
sured  target  data.  For  the  reasons  given  in  Sec.  1.1.1,  we  require  that  the  model  be  accurate  to 
second-order  statistics;  that  is,  the  model  must  accurately  represent  the  first-order  probability  density 
function,  the  power  spectral  density  function  (and  equivalently,  the  autocorrelation  function),  and  asso¬ 
ciated  parameters  over  intervals  where  the  process  is  piecewise-stationary.  We  do  not  require  that  the 
second-order  probability  density  function  be  specified  because  that  function  is  difficult  to  estimate  prac¬ 
tically  and  though  it  more  completely  specifies  the  process,  it  is  less  useful  in  characterizing  signal 
time-correlation  properties  than  is  the  power  spectral  density. 

The  primary  targets  of  concern  in  this  report  are  surface  ships.  Relative  to  aircraft  modeling,  ship 
modeling  has  received  limited  attention  in  the  literature.  Ships  are  much  more  complex  in  structure 
than  aircraft  and  scattering  from  them  is  complicated  by  sea  surface  multipath.  However,  the  model 
developed  is  of  a  general  nature  and  should  be  useful  in  problems  where  the  assumptions  made  in  its 
development  hold. 

1.2  Previous  Work 

Variations  in  RCS  as  a  function  of  time  (amplitude  scintillation)  were  first  observed  during  World 
War  II  (51.  The  early  models  of  amplitude  scintillation  were  developed  for  application  in  the  detection 
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of  aircraft  targets  beginning  with  the  work  of  Marcum  [6]  and  Swerling  [7,8],  These  models  represent 
the  target  as  a  point  source  whose  fluctuation  statistics  are  chosen  depending  on  the  type  of  radar  used 
and  the  measured  or  assumed  fluctuation  properties  of  the  target.  During  the  1960s,  numerous  papers 
and  books  appeared  describing  various  methods  of  /V-source  modeling  of  RCS.  The  August  1965  issue 
of  the  IEEE  Proceedings  [9]  was  devoted  to  radar  reflectivity  and  is  a  good  reference  for  the  early  RCS 
work. 


The  glint  phenomena  were  not  observed  until  radar  angle  tracking  capabilities  improved  in  the  late 
1940s.  The  first  models  of  glint  appeared  in  Mead  et  al.  [10],  De  Lano  [11],  and  Howard  [12]  during 
the  1950s.  Each  was  based  on  the  /V-sourct  concept. 

We  focus  our  review  on  models  that  are  based  on  the  /V-source  concept  because  that  is  the  con¬ 
cept  on  which  the  analysis  of  this  report  is  based. 

1.2.1  Random  Models 

One  method  of  modeling  complex  targets  has  been  to  represent  the  target  as  an  /V-source  model 
with  scatterers  whose  individual  scattering  properties  are  assumed  to  be  random.  When  simple  random 
properties  are  assumed,  this  type  of  model  allows  an  analytical  solution  to  be  obtained  for  the  target 
scattering.  This  tic thod  was  used  by  Muchmore  [13,14]  to  describe  the  RCS  of  aircraft  and  reasonable 
spectral  estimates,  relative  to  measured  data,  were  obtained.  However,  as  pointed  out  by  Peters  and 
Weimer  [15-17],  this  method  has  serious  drawbacks  when  used  to  analyze  radar  tracking  performance. 
This  is  because  over  short  time  and  aspect  intervals  the  individual  scatterers  do  not  behave  indepen¬ 
dently. 

This  random  method  was  used,  for  arbitrary  complex  targets,  by  De  Lano  [11],  Gubonin  [18], 
and  Mumford  [19]  to  study  glint  statistics;  by  Mohanty  [20],  Gruner  [21],  and  Mitchell  [22,23]  to  study 
RCS  statistics;  and  by  Varshavchik  [24]  and  Borden  [25]  to  study  the  combined  RCS-glint  problem. 
Baras  [26]  used  an  /V-source  representation  of  the  target  to  obtain  the  fundamental  parameters  of  an 
equivalent-point-source  model  of  glint. 

Jakeman  [27,28]  and  Jakeman  and  Pusey  [29-32]  have  investigated  using  AT -distributions  [27]  to 
describe  the  amplitude  statistics  of  scattered  radiation  in  a  variety  of  experiments  involving  scattering 
from  turbulent  media  (e  g.,  sea  clutter  and  optical  scintillation).  This  A -distribution  model  arises  when 
it  is  assumed  not  only  that  the  individual  sources’  amplitudes  and  phases  are  independent  random  vari¬ 
ables  but  also  that  the  number  of  sources  is  a  random  variable. 

1.2.2  Deterministic  Models 

The  most  common  type  of  model  of  complex  targets  based  on  the  /V-source  concept  assumes  that 
the  target  is  composed  of  scatterers  with  simple  geometric  shape  for  which  scattering  solutions  are 
known.  The  field  scattered  by  the  target  is  then  the  sum  of  the  individual  scattered  fields  for  the  given 
radar-target  geometry.  Physical  optics  theory,  the  geometric  theory  of  diffraction,  integral  equation 
methods,  and  combinations  thereof  are  used  to  obtain  scattering  solutions  for  the  component  scatterers. 
Discussions  of  these  methods  can  be  found  in  many  papers  and  books;  examples  are  Crispin  and  Siegel 
[33],  Bechtel  and  Ross  [34],  Bowman  et  al.  [35],  Ruck  et  al.  [36],  Keller  [37],  and  Oshiro  et  al.  [38,39]. 
Many  examples  of  aircraft  and  missile  modeling  are  contained  in  or  referred  to  by  the  previous  refer¬ 
ences.  Examples  of  ship  models  using  this  approach  are  given  by  Toothman  [40]  and  Radza  and  Stenger 
[41],  We  note  that  in  modeling  ships,  thousands  of  scatterers  are  often  required  and  usually  up  to 
twenty  generic  types  of  scatterers  are  used.  Further,  multipath  effects  on  each  scatterer  and  multiple 
scattering  must  be  accounted  for. 
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The  major  studies  of  glint  have  been  performed  without  assigning  specific  properties  to  the  indivi¬ 
dual  scattered  fields.  This  approach  has  led  to  generic  representations  for  glint.  Howard  [12]  used  this 
approach  to  show  that  glint  was  equivalent  to  the  slope  of  the  phase  front  of  the  reflected  field  In  the 
radar  angular  coordinate  of  interest.  Lindsay  [42]  expanded  the  concept  of  phase-front  slope  to  phase- 
front  gradient.  He  also  showed  the  relationship  between  glint  and  Doppler  scintillation.  Dunn  and 
Howard  [43]  showed  that  glint  was  equivalent  to  Poynting  vector  tilt  and,  independent  [5]  of  Lindsay, 
showed  the  relationship  between  glint  and  Doppler  scintillation. 

1.2.3  Unified  Scattering  Models 

Very  few  attempts  have  been  made  to  develop  a  unified  RCS  and  glint  model  for  the  time-varying 
return  from  a  complex  target.  The  geometric  models  described  above  can  be  used  to  generate  time 
series  of  data  by  moving  the  target  randomly  and  solving  for  the  scattered  field  as  a  function  of  time. 
Such  a  procedure  is  of  course  extremely  time  consuming  even  on  a  computer,  especially  for  ships. 
Wright  and  Haddad  [5]  developed  such  a  model  for  an  airborne  drone  where  the  individual  scatterers 
were  represented  as  ellipsoids.  Both  RCS  and  two-dimensional  glint  were  modeled.  The  stochastic 
returns  from  several  target  aspects  were  obtained  and  analyzed,  assuming  appropriate  random  motion  of 
the  drone.  Borden  [25]  developed  a  random  model  for  the  unified  target  return  of  an  aircraft,  assum¬ 
ing  only  one-dimensional  glint.  He  modeled  the  target  as  N  sources  and  then  assigned  statistical  values 
to  the  individual  scatterer  amplitudes,  phases,  and  Doppler  frequencies.  Varshavchik  [24]  investigated 
the  relationship  between  the  amplitude  and  phase  characteristics  of  targets  composed  of  an  array  of 
identical  isotropic  scatterers.  He  studied  returns  resulting  from  small  angle-oscillations  of  the  targe*  in 
a  fixed  plane. 

1.2.4  Multipath  Models 

Low  grazing-angle,  forward-scatter,  over-water  multipath  has  been  studied  extensively  during  the 
past  quarter  century.  The  most  common  representation  of  such  scattering  has  been  made  in  terms  of 
equivalent-point-source  models  where  the  scattered  field  is  viewed  as  being  composed  of  two  com¬ 
ponents:  a  "specular"  (or  "coherent")  component  that  is  deterministic,  and  a  "diffuse"  (or  "incoherent") 
component  that  is  random.  Although  it  is  well  known  that  the  specular  component  appears  to  a  radar 
as  coming  from  a  point-source  that  is  located  at  the  geometric  image  of  the  transmitting  point-source, 
the  apparent  point-source  location  (s)  of  the  diffuse  component  (s)  is  not  well  understood.  This  lack  of 
understanding  has  caused  difficulty  in  analyzing  the  tracking  errors  induced  by  diffuse  multipath 
[44,45], 

In  a  series  of  papers,  Beard  et  al.  [46-48]  developed  a  second-order,  statistical  model  of  the  scat¬ 
tered  field.  The  model  is  based  on  empirical  measurements  and  the  fundamental  parameters  are  diffuse 
and  specular  scattering-coefficients.  Beard’s  model  is  valid  in  the  sense  of  predicting  the  received 
power  at  microwave  frequencies.  It  implicitly  assumes  that  the  diffuse  scattering  can  be  viewed  as  com¬ 
ing  from  a  point-source  located  at  the  specular  point.  Northam  [49]  presents  a  second-order  stochastic 
model,  based  on  this  work,  that  is  useful  for  simulating  multipath  effects  as  a  function  of  time.  Beck¬ 
mann  and  Spizzichino  [50]  theoretically  predicted  that  the  diffuse  scattering  arose  primarily  from  sur¬ 
face  areas  near  the  transmitter  and  near  the  receiver.  Barton  [51]  modified  the  diffuse  coefficient  that 
was  derived  by  Beckmann  and  Spizzichino  with  a  roughness  factor  and  formulated  the  coefficient  in 
terms  of  a  bistatic  scattering  parameter  for  the  surface.  He  showed  that  diffuse  scattering  arose  from 
surface  areas  that  vary  as  a  function  of  geometry  and  surface  structure.  Mrstik  and  Smith  [44]  investi¬ 
gated  the  limitations  caused  by  multipath  on  low-angle  tracking  using  four  different  models  of  the 
poorly  understood  bistatic  scattering  parameter.  Smith  and  Mrstik  [45]  investigated  multipath-induced 
tracking  errors  in  elevation  scanning  and  mcnopulse  radars  by  assuming  that  the  diffuse  tracking  errors 
are  small  so  that  the  diffuse  power  can  be  viewed  as  a  noise-like  interference  superimposed  on  the 
direct-path  signal.  Baras  [52]  has  summarized  the  various  models  from  the  perspective  of  stochastic 
modeling. 
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2.  SUMMARY 

This  report  addresses  the  modeling  of  the  effects  of  electromagnetic  scattering  from  complex,  dis¬ 
tributed  targets  in  a  multipath  environment.  The  basic  objectives  are  to  develop  a  scattering  model  that 
is  useful  in  analyzing  and  simulating  the  correlated  processes  of  RCS  and  glint,  and  to  develop  the 
equations  of  these  processes  appropriate  to  that  model.  The  model  is  specifically  intended  to  be  useful 
for  studying  the  tracking  performance  of  pulse-radars.  Ships  are  the  major  targets  of  interest.  We  are 
primarily  concerned  with  representing  scattering  over  short-time  intervals.  A  digital  simulation  is 
developed,  and  some  implications  of  results  from  that  simulation  are  discussed. 

2.1  Approach 

Our  approach  is  to  use  a  (largely)  deterministic  model  for  the  target  scattering  and  a  stochastic 
model  of  the  target’s  motion.  Using  the  target’s  motion  as  input  we  obtain  a  stochastic  representation 
of  the  received  field.  Over-water  multipath  effects  are  also  modeled  as  stochastic  processes. 

We  model  the  target  as  a  finite  collection  of  point-source  scatterers  that  represent  the  dominant 
scattering  components  of  the  target.  These  "unit-scatterers"  are  defined  to  allow  a  straightforward 
representation  of  the  received  field  at  the  radar.  Because  the  target  is  represented  as  a  collection  of 
point-source  scatterers,  we  draw  on  existing  models  of  RCS  and  glint  that  were  developed  assuming 
point-source  scattering.  A  simulation  of  correlated  RCS  and  glint  is  developed  for  use  in  studying  the 
effects  of  model-parameter  variations. 

2.2  Outline 

The  concept  of  unit-scatterer  is  introduced  in  Section  3,  and  an  analytic  representation  of  the  con¬ 
cept  is  presented.  Methods  of  identifying  and  measuring  the  target’s  component  unit-scatterers  are  pro¬ 
posed  and  anticipated  properties  of  the  associated  amplitudes  and  phases  are  discussed. 

In  Section  4,  the  usual  model  of  the  received  field  from  a  collection  of  point-source  scatterers  is 
presented,  and  for  completeness  a  model  of  a  generic  radar  system  is  developed  in  the  context  of 
point-source  scatterers.  A  point-source  multipath  model  is  then  incorporated  into  the  equation  that 
represents  the  received  field. 

The  motion  equations  for  the  individual  unit-scatterers  of  a  ship  are  derived  in  Section  5.  For 
completeness,  two  formulations  of  these  equations  are  made,  each  assuming  different  inputs:  the  first 
assumes  knowledge  of  the  driving  forces  and  moments  that  are  applied  to  the  ship;  the  second  assumes 
a  frequency-response  model  of  the  ship  motion. 

Section  6  describes  a  ship-motion  simulation  that  we  developed  as  the  stochastic  input  to  the 
scattering  simulation.  It  is  based  on  the  cataloged  ship-motion  data  generated  by  the  David  Taylor 
Naval  Ship  Research  and  Development  Center  (DTNSRDC)  using  frequency-response  models  for  the 

ship  motions. 

In  Section  7,  equations  for  RCS  are  developed  in  terms  of  unit-scatterers  based  on  the  formula¬ 
tions  of  Section  4.  Because  tracking  radars  process  data  over  short  time-intervals,  RCS  variations  over 
such  intervals  are  analyzed,  and  the  nonstationarity  of  the  resulting  processes  is  illustrated. 

In  Section  8,  two  approaches  to  glint  modeling  are  used  to  deve'  p  correlated  azimuth  and  evalua¬ 
tion  glint  models  for  targets  that  are  modeled  by  unit-scatterers.  ’  is  shown  that  one  approach  is  more 
useful  than  the  other,  especially  when  point-source  multipath  effects  are  included. 

Section  9  presents  the  results  of  simulation  studies  The  effects  of  varying  some  of  the  model 
parameters  are  investigated 
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3.  THE  UNIT-SCATTERER  CONCEPT 

As  described  in  Section  1,  one  of  the  most  common  approaches  to  modeling  complex  targets  is  to 
represent  them  as  being  comprised  of  a  finite  number  of  point-source  scatterers.  Using  this  approach, 
various  scattering  models  have  been  developed.  The  differences  in  these  models  are  due  to  the  dif¬ 
ferent  ways  in  which  the  scatterers  are  characterized.  We  use  this  approach  but  attempt  to  minimize 
the  number  of  scatterers  and  the  statistical  assumptions  made  regarding  the  parameters  of  those  scatter¬ 
ers. 

Determining  the  characteristic  parameters  of  a  scatterer  given  the  incident  and  scattered  fields  is 
referred  to  as  an  inverse  scattering  problem.  Recent  years  have  seen  a  rapid  growth  of  interest  in  such 
problems  in  many  areas  of  science,  especially  in  optics  and  electromagnetics;  however,  for  the  EM 
inverse  problem,  "the  prospects  of  expressing  scattering  characteristics  of  complex  shapes  successfully 
and  simply  are  still  inadequate"  [53], 

3.1  Scattering  Matrix 


When  an  object  is  illuminated  by  a  plane  wave  (implied  by  large  radar- to-scatterer  range),  the 
scattered  field,  £*,  can  be  viewed  as  a  linear  transformation  of  the  incident  field,  E‘ .  The  matrix 
representation  of  this  transformation  is  called  a  scattering  matrix  [33].  Such  matrices  have  the  form 


S- 


Su  exp 

Sn  exp(/(J»2i> 


S 12  exp  0^12) 
S22  exp  0>22) 


(3.1) 


(We  note  that  due  to  phase  being  relative,  the  conservation  of  energy  principle,  and  the  reciprocity 
theorem,  only  five  of  the  eight  parameters  of  S  need  to  be  determined  to  completely  specify  S  [54].) 
The  scattered  field  representation  is 

Es  -  SE'.  (3.2) 


The  scattering  matrix  provides  a  significant  characterization  of  radar  targets;  knowledge  of  it  for  a  given 
set  of  radar  parameters  allows  a  complete  characterization  of  the  target’s  RCS  for  those  parameters. 
Further,  the  scattering  matrix  for  any  polarization  pair  can  be  obtained  from  the  matrix  expressed  in 
terms  of  any  two  arbitrary  polarization  vectors  that  are  orthogonal  [33].  As  emphasized  by  Boerner 
[53]  in  his  paper  discussing  the  utilization  of  polarization  in  the  EM  inverse  scattering  problem,  when 
measuring  target  scattering  all  components  of  the  scattering  matrix  should  be  estimated  or  a  loss  of 
information  will  result. 


Given  a  radar  with  receiver  polarization  pR,  the  RCS  of  the  object  can  be  defined  to  be  [54] 

rj  ,  ;  1 7 

E  Pr 


(t  =  lim  47 rr2 


(3.3) 

—  lim  47rr2 

r—  00 

where  r  is  the  target-to-radar  range.  If  the  transmitter  has  the  polarization  pT ,  then  we  can  write  E‘  = 
\E'\pr  and 


SE1  ■  pR 
\E'\ 


<7  -  lim  4i rr2(SpT  ■  pR)2. 

r— 00 


(3.4) 
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In  general,  polarizations  are  defined  in  terms  of  horizontal  and  vertical  unit  vectors,  p(H)  and 
p(  V).  Using  these  polarizations,  we  write  the  scattering  matrix  as 


S 


S(HH)  S(HV) 
S(  VH)  S(VV)  ■ 


(3.5) 


It  follows  from  Eqs.  (3.4)  and  (3.5)  that  there  are  four  fundamental  components  that  characterize  the 
RCS  of  scattering  objects;  they  are 

<r(/y)  =  lim  4nr2\S(IJ)\\  /,  J  =  //,  V.  (3.6) 


Therefore,  the  RCS  of  an  object  can  be  described  by  the  matrix 

r  _  \cr(HH)  o-(HV)  1 
£  ~  I (t{VH)  «r(m|- 


(3.7) 


The  above  concept  of  scattering  can  be  extended  to  targets  composed  of  individual  scatterers  for 
which  scattering  matrices  are  known.  The  procedure  is  to  sum  the  components  of  the  individual 
scattering  matrices  to  yield  a  target  scattering  matrix;  relative  phasing  and  multiple  scattering  between 
the  scatterers  must,  of  course,  be  taken  into  account.  The  target  RCS  can  then  be  determined  from  the 
resulting  matrix.  However,  the  complexity  of  this  procedure,  due  to  the  relative  phasing,  has  led  to  the 
use  of  two  assumptions  to  simplify  the  problem.  Because  an  average  value  of  RCS  is  often  all  that  is 
required  to  be  known  about  the  target,  it  is  often  assumed  that  the  individual  scatterers  are  phase 
independent.  This  allows  a  straightforward  estimate  of  the  RCS  (cross  products  of  independent  phases 
cancel).  Second,  because  most  radars  operate  using  only  one  polarization,  it  is  often  assumed  that  the 
three  scattering  components  associated  with  the  orthogonal  polarization  can  be  ignored.  Both  of  these 
assumptions  have  proved  quite  useful  for  practical  problems. 

3.2  Definition 

It  is  well  known  that  complex  targets  produce  scattering  that  appears  to  result  from  several 
apparent  "sources"  ("dominant  scatterers,"  "bright  spots,”  "hot  spots,”  and  "flare  spots”)  that  are  located 
at  "scattering  centers."  The  number,  amplitude,  and  phase  of  these  sources  vary  with  target  aspect  and 
radar  frequency,  polarization,  and  resolution.  Examples  of  this  phenomenon  for  measured  targets  can 
be  found  in  Refs.  55  through  59.  These  apparent  sources  result  because  scattering  arises  not  from  the 
target’s  entire  surface  but  from  points  of  discontinuity  and  specular  points  of  that  surface;  cancellation 
of  the  field  occurs  along  continuously  varying  surfaces  of  the  target  (341 .  Kell  [60]  refers  to  the  portion 
of  the  target  surface  that  is  near  a  scattering-center  and  that  produces  the  return  associated  with  that 
scattering-center  as  the  "cophased  area."  The  scattering-centers  are  not  necessarily  due  to  specific 
geometric  objects;  they  may  arise  from  phasing  between  adjacent  objects  whose  returns,  when  summed 
incoherently,  would  not  be  significant  (341.  This  phenomenon  may  lead  to  scattering-centers  that  are 
not  located  on  the  target. 


In  light  of  the  above,  we  define  a  unit-scatterer  (US)  to  be  an  apparent  point-source  scatterer  as 
perceived  by  a  radar.  Clearly,  the  definition  is  radar  dependent;  such  objects  depend  directly  on  radar 
frequency  and  pulse  width.  We  characterize  a  US  by  a  complex  amplitude-matrix,  X.  Our  definition  of 
this  matrix  is  analogous  to  that  of  the  scattering  matrix: 


A  ,|  exp  0<£n) 

T- 

A2\  exp  (j<f>2\) 


A 12  exp  C/^12) 
A 2 2  exp  (J<t>22) 


(3.8a) 


where 
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The  amplitudes  and  phases  are  functions  of  US  orientation,  and  the  phase  is  also  a  function  of  radar-US 
range.  The  amplitude  functions,  AtJ,  are  the  square  roots  of  the  associated  component  RCSs.  The 
phase  functions,  i>0,  are  defined  relative  to  some  reference  point  on  the  radar-to- target  sight  line;  this 
allows  coherent  summing  of  the  matrices  that  represent  the  target.  For  the  standard  polarizations,  Eq. 
(3.8a)  becomes 


X 


A hh  e*P  U^hh)  AHv  exp  (J4>hv ) 
A  YH  exp  (j<t>  VH )  A  VV  exp  ( j<t>  yy) 


(3.9) 


The  RCS  of  a  US  that  is  characterized  by  X  becomes 

<r-\pfXpR\2  (3.10a) 

where  Pt  is  the  transmitter  polarization  and  pR  is  the  receiver  polarization.  For  a  target  composed  of  N 
USs, 

(3.10b) 


ZpIXiP* 


The  significance  of  X  is  that  the  electric  field  at  the  receiver  antenna  and  associated  received  signals 
can  be  written  (see  Section  4)  in  terms  of  a  scattering  function  T(\  •)  which  we  define  by 


f  (j 7r>  Pr) 


*Pt  jix,  Pr- 


O.U) 


Although  this  approach  to  modeling  scattering  is  not  unique  (similar  formulations  to  obtain  RCS  and 
glint  are  implicit  in  other  work),  the  definition  of  the  unit-scatterer  is  new  in  that  it  is  based  on  the 
fundamental  scattering  objects  being  large,  complex,  and  not  necessarily  physically  identifiable. 


Because  of  the  complexity  of  USs  comprising  large  targets,  such  as  ships,  it  is  likely  that  experi¬ 
mental  methods  or  detailed  simulation  studies  are  necessary  to  accurately  determine  the  AT,  matrices. 
However,  as  is  discussed  in  Sec.  3.4,  relatively  simple  representations  of  the  USs  may  yield  good 
approximations  to  the  components  of  X,. 


3.3  Identification 


Given  the  above  definition  of  the  US  we  now  discuss  methods  of  identifying  the  USs  on  complex 
targets.  There  are  two  basic  approaches  to  US  identification  (and  characterization):  using  analytic  or 
experimental  methods.  We  first  consider  several  analytic  methods. 

Kell  [60]  investigated  bistatic  RCS  using  the  ^-source  formulation.  In  this  and  the  next  para¬ 
graph,  we  discuss  that  work  as  it  relates  to  the  monostatic  problem.  Using  the  Straton-Chu  integral  in 
the  definition  of  RCS,  he  showed  that  a  target’s  bistatic  RCS  could  be  written  in  the  form 

o-  -  -jj  |/  (3.12) 

where 


k0  is  the  wavenumber  of  the  incident  field, 
0  is  the  bistatic  angle,  and 

A  is  the  radar's  wavelength. 
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Although  Eq.  (3.1)  is  an  exact  expression,  the  difficulty  with  using  it  is  that  /(z)  is  rarely  known  pre¬ 
cisely  (60).  However,  the  analytic  continuity  of  /(z)  allows  it  to  be  divided  into  a  sum  of  integrals  over 
subregions  of  the  target  and  "reasonable  approximations"  to  /(z)  lead  to  contributions  from  these 
integrals  only  near  the  end  points  of  each  integral  [60).  These  contributions  are  then  identified  as  the 
scattering  centers.  He  distinguishes  between  two  types  of  scattering-centers:  "simple  centers,"  which 
result  from  direct  illumination  only,  and  "reflex  centers,"  which  result  from  multiple  reflections.  He 
states  that  it  is  common  for  reflex  scatterers  to  produce  much  larger  RCS  values  than  simple  scatterers 
do. 

Once  the  sources  are  identified,  the  target  RCS  can  be  written  as  a  finite  sum  of  those  sources. 
Kell  shows  that  in  terms  of  N  sources,  the  bistatic  RCS  can  be  expressed  as  (assuming  fixed  polariza¬ 
tions) 

o’  -  £  eJ*' 

i- 1 

with 

<t>,  -  2 k0r,  cos  ^  (3.13b) 

where 

r,  is  the  distance  between  the  /th  and  first  sources’  phase-centers,  projected  on  the  bistatic 
axis,  and 

is  the  "residual  phase  contributions"  of  the  i<h  center. 

Setting  /3  =  0  yields  the  monostatic  RCS.  (We  note  that  to  derive  the  monostatic-bistatic  equivalence 
theorem  in  Ref.  60.  Kell  assumes  that  the  do  not  vary  over  the  range  of  /3  considered;  we  suggest 
that  for  USs  arising  from  structurally  complex  objects,  the  may  vary  significantly  as  a  function  of  the 
monostatic  aspect  angle.) 

When  the  magnetic  field  H  on  the  surface  of  a  perfectly  conducting  object  is  known  and  the  field 
point  is  at  a  large  distance  from  the  body,  the  scattered  magnetic  field,  Hs,  can  be  approximated  by 
[61] 

Hs  ~  f  (n  x  H)  x  re~Jkr  T‘  dS\  (3.14) 

4ita  js 

where  n  is  the  outward  normal  to  the  surface,  k  =  2 n/k,  f  is  the  unit  vector  from  the  origin  to  the 
field  point,  R  is  the  distance  from  the  origin  to  the  field  point,  S'  is  the  scattering  object’s  surface,  and 
7'  is  the  radius  vector  from  the  origin  to  the  integration  point  on  the  surface.  To  solve  Eq.  (3.14)  we 
must  know  H  and  S'^  Although  obtaining  S'  for  a  known  object  is  straightforward,  though  perhaps 
tedious,  determining  H  for  an  arbitrarily  shaped  object  is  in  general  not  feasible.  However,  for  objects 
that  are  relatively  smooth  (large  radii  of  curvature  relative  to  the  radar’s  wavelength),  it  is  well  known 
that  the  KirchhofMphysical  optics)  approximation  (i.e.,  assume  that  the  surface  current  vvnsity  is 
-2  h  x  H'  where  H'  is  the  incident  magnetic  field)  leads  to  an  accurate  description  of  the  scattered 
field.  Assuming  an  incident  field  of  the  form 

£ 

H‘  =  _o  cr  x  e)e-jk?  7'  (3.15) 

V 

where  e  is  the  polarization  direction  for  the  incident  electric  field  E‘ ,  17  is  the  characteristic  impedance 
of  free  space,  and  E0  =  \E‘  I,  the  scattered  electric  field  can  be  approximated  by 

£*  =  JV  —  l«  x  V  *  e)l  *  re~2jk‘r  T  dS'.  (3.16) 

2nR  J S  Tj 


(3.13a) 
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If  the  surface  can  be  separated  into  N  independent  (no  multiple  scattering)  surfaces  for  which  solutions 
to  Eq.  (3.16)  can  be  obtained,  then  the  object’s  scattered  field  can  be  described  by  the  sum  of  N  fields 
(due  to  the  N  "elementary"  scatterers).  In  general,  N  will  be  large  and  multiple  scattering  cannot  be 
ignored  so  that  the  resulting  scattering  model  will  be  very  complicated  (see  e.g..  Ref.  40).  Scattering 

centers  for  complex  objects  will  arise  from  one  of  two  situations:  the  return  from  an  elementary 

scatterer  dominates  (in  magnitude)  the  return  from  most  of  the  other  elementary  scatterers,  or  several 
elementary  scatterers  produce  scattered  fields  that  add  in  phase  yielding  a  total  field  that  dominates 
most  of  the  other  returns.  However,  if  we  cannot  make  the  above  separation  but  we  can  represent  the 
(known)  surface  in  cartesian  coordinates  by 

z-f(x,y)  (3.17) 

so  that 

7'  -  xi  +  yj  +  /(x,  y)k,  (3.18) 

we  can  apply  the  method  of  stationary  phase  to  the  integral  in  Eq.  (3.16)  if  we  can  find  the  points  on 
the  surface  where 

T-  (r  7') (r  •  7’)  «  0.  (3.19) 


The  method  of  stationary  phase  is  an  asymptotic  expansion  technique  for  evaluating  integrals  of 
the  form 

/(*)  -  f*  g(z)e»hh)dz,  (3.20) 

^  a 

where  x  is  a  large  positive  parameter  and  h(z)  is  a  real  function  of  the  real  variable  z  [62].  The  major 
contribution  to  the  integral  arises  from  the  immediate  vicinity  of  the  end  points,  and  at  stationary 
points  of  h  (z)  (points  where  h'(z)  —  0)  and  in  the  first  approximation,  the  contribution  of  the  station¬ 
ary  points  is  more  important  than  the  contribution  of  the  end  points  [62].  If  z,  /  —  1,  ....  n  are  the 
points  of  stationary  phase  (i.e.,  where  h(z)  has  extrema),  then  the  integral  Eq.  (3.16)  can  be  approxi¬ 
mated  by 

/ Or)  ~  T  /  '+<  g(z)eJxl,(z>  dz  (3.21) 

i-i  z‘  * 

where  the  «,  are  small.  The  /'th  element  of  the  summation  is  evaluated  by  expanding,  h(z )  about  z. 
The  method  applies  to  the  integral  Eq.  (3.16)  since  for  wavelengths  of  interest  to  us,  2k  is  large  (on 
the  order  of  400). 

Kodis  [63]  investigated  the  scattering  from  a  random,  perfectly  conducting  irregular  surface  whose 
radii  of  curvature  are  continuous  and  large.  Evaluating  the  scattered-field  integral  by  using  the  method 
of  stationary  phase,  he  showed  that  to  a  first  approximation  the  RCS  of  the  surface  is  proportional  to 
the  average  number  of  specular  points  which  are  illuminated.  Although  he  assumed  that  there  was  no 
multiple  scattering,  implying  that  the  principal  radii  of  the  surface  have  a  lower  bound  that  is  much 
greater  than  the  wavelength  of  the  scattered  radiation  [63],  his  results  also  suggest  that  the  scattered- 
field  integral  can  be  approximated  by  a  finite  sum  and  that  the  method  of  stationary  phase  may  prove 
useful  in  identifying  scattering  centers  of  objects  that  are  relatively  smooth. 

The  method  of  stationary  phase  is  of  interest  because  it  identifies  scattering  centers  on  the  scatter¬ 
ing  object.  When  the  method  can  be  used  we  anticipate  that  the  scattering  centers  identified  by  it  will 
be  USs  of  the  object.  The  method  may  yield  not  only  US  locations  but  also  analytic  representations  for 
those  USs;  by  investigating  the  scattered  field  at  these  points,  we  may  be  able  to  develop  analytic 
representations  for  USs  resulting  from  "smooth  surface"  scattering.  Examining  the  structure  of  ships  we 
see  that  in  general  their  surfaces  are  not  smooth.  In  fact  there  are  many  flat  surfaces  connected  at  right 
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angles  (zero  radius  of  curvature).  Therefore  the  method  of  stationary  phase  will  not  in  general  be  use¬ 
ful  for  identifying  USs  on  targets  such  as  ships.  However,  the  method  may  prove  useful  for  targets, 
such  as  aircraft,  whose  surfaces  are,  in  general,  "smooth." 

Precise  identification  of  the  USs  for  a  specific  target  will  most  likely  require  actual  target  measure¬ 
ments.  Such  measurements  must  be  made  by  some  type  of  high  resolution  radar.  Possibilities  include 
narrow  pulse  radar,  pulse  compression  radar,  synthetic  aperture  radar  (SAR),  and  inverse  synthetic 
aperture  radar  (ISAR).  Such  radars  must  measure,  as  a  function  of  target  aspect,  the  return  signal 
from  the  target  in  cells  that  are  small  relative  to  the  target's  largest  dimension.  Although  our  studies 
will  show  that  even  large  ship  targets  can  be  well  represented  at  a  given  aspect  by  as  few  as  six  USs,  we 
suggest  a  cell  size  at  least  as  small  as  1/10  to  1/20  of  the  target’s  largest  dimension.  The  cell  size  can¬ 
not  be  too  small  because  for  studies  of  radars  with  less  resolution  than  the  measured  data,  the  smaller 
scatterers  can  be  grouped  (summed  coherently  or  incoherently)  to  produce  USs  appropriate  for  the 
desired  resolution.  High  resolution  data  exists  on  many  targets  including  ships  but  analysis  of  that  data 
has  not  been  done  as  part  of  this  research;  the  effort  required  to  obtain  and  analyze  such  data  is 
beyond  the  scope  of  this  report.  High  resolution  radars  identify  scatterers  in  range  and  amplitude  for 
the  associated  radar-to-target  aspect.  Figure  3.1  illustrates  such  a  range  profile  for  an  aircraft  model; 
for  this  aspect  (nose-on)  there  are  five  dominant  scatterers.  By  rotating  the  target  or  moving  the  radar, 
the  amplitudes,  relative  phases,  and  position  of  the  dominant  scatterers  can  be  determined  as  a  function 
of  aspect.  Care  must  be  taken  to  vary  the  aspect  by  increments  small  enough  to  yield  accurate 
(smooth)  representations  of  these  parameters.  Figure  3.2  illustrates  scatterer  data  as  a  function  of 
aspect  for  a  simple  target. 


Fig.  3.1  —  A  high  range-resolution  radar's  output  (from  Ref.  64) 
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Fig.  3.2  —  Amplitude-position  data  for  a  high  range-resolution  radar  as  a 
function  of  aspect  (from  Ref.  65) 


Another  method  of  identifying  the  USs  is  to  use  existing  geometric-model  simulations.  (Various 
organizations,  including  the  Naval  Research  Laboratory  (NRL),  Georgia  Institute  of  Technology,  Tech¬ 
nology  Services  Corporation,  and  Rockwell  International,  have  developed  such  simulations.)  This 
method  involves  running  such  a  simulation  for  the  target  of  interest  and  for  the  desired  radar-target 
geometry.  The  resulting  scattered  field  would  then  be  examined  to  identify  the  major  scatterers  (say, 
those  that  contribute  95%  of  the  total  return)  and  then  grouped  as  appropriate  into  USs.  Similarly, 
SAR  or  ISAR  simulations  could  also  be  used  to  determine  the  USs.  Figure  3.3  is  an  illustration  of  the 
resulting  SAR  image  from  a  simulation  of  a  KC-135  aircraft  for  various  SAR  resolutions.  The  "blobs" 
in  these  images  are  candidate  USs. 

A  less  accurate  but  very  simple  method  of  identifying  some  USs  is  to  examine  geometric  images 
of  the  target  (photographs,  line  drawings,  blueprints,  etc.).  This  method  is  useful  because  scattering 
from  certain  generic  objects  is  known.  For  example,  Table  3.1  lists  the  22  objects  used  by  Toothman 
[40)  in  a  geometric  type  of  simulation  of  ship  scattering.  An  examination  of  the  structure  of  ships 
shows  that  there  often  are  obvious  generic  scatterers  whose  dimensions  are  such  that  they  will  dominate 
the  scattering  from  the  ship  at  specific  aspect  angles.  For  example,  the  broadside  and  stern  aspects  of 
the  DD963  (Fig.  3.4)  illustrate  the  presence  of  flat-plate  scatterers  in  the  superstructure.  The  aft- 
quarter  aspects  indicate  the  presence  of  large  corner  reflectors  and  the  two  exhaust  stacks  should  have 
reflections  that  are  analogous  to  those  from  cylinders.  Using  these  observations  we  can  obtain  a  set  of 
USs  that,  to  first-order,  approximate  the  actual  USs  for  a  given  aspect.  This  is  the  method  that  we  will 
use  to  obtain  USs  for  the  simulation  (Section  9).  In  the  next  two  sections  we  discuss  in  some  detail  the 
application  of  this  method  to  obtain  US  amplitude  and  phase  characterizations. 

3.4  Amplitude 

The  measurement  and  simulation  procedures  of  the  previous  section  yield  amplitude  information 
and  can  be  used  to  determine  the  amplitude  functions.  The  visual  method,  however,  does  not  allow 
such  an  accurate  determination  except  for  obvious  geometric  objects.  For  analysis  and  simulation  pur¬ 
poses  we  would  like  to  obtain  closed  form  expressions  for  the  amplitude  functions.  Such  expressions. 
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Table  3.1  —  Generic  Scattering  Types 
Used  by  Toothman  [40] 


Concave  Edge 

Rectangular  Plate 

Convex  Edge 

Paraboloid 

Edge  Caustic 

Ogive 

Elliptic  Disc 

Point 

Ellipsoid 

Concave  Dihedral  (2  reflections) 

Hyperboloid 

Concave  Trihedral  (3  reflections) 

Cylinder 

Straight  Edge  (convex  dihedral) 

Elliptic  Cone 

Convex  Trihedral  (3  plane  tip) 

Inner  Torus 

General  Curved  Surface  Specular 

Outer  Torus 

Cavities 

Elliptic  Tip 

Antennas 

if  accurate,  would  allow  the  US  model  to  be  more  useful  in  general;  the  model  would  not  be  depen¬ 
dent  on  detailed  information  about  the  target.  Therefore,  our  goal  in  specifying  the  amplitudes  is  to 
obtain  simple,  visually  identifiable,  closed  form  expressions  that  are  valid  over  short  time-intervals.  To 
achieve  this  we  rely  on  optical  scattering  theories  and  RCS  measurements  of  distributed  targets  to  guide 
our  choice  of  amplitude  functions. 

It  is  well  known  that  the  spectrum  of  the  RCS  of  a  distributed  target  is  band-limited  and  it  is 
often  characterized  as  being  low-pass.  This  implies  that  a  finite  (though  perhaps  large)  number  of 
scatterers  can  be  used  to  accurately  approximate  the  RCS.  It  is  also  known  that  for  most  targets,  due  to 
radar-target  relative  motion,  the  target  RCS  is  not  a  stationary  process  because  the  mean,  variance,  den¬ 
sity  type,  and  spectral  properties  all  vary  with  aspect  as  indicated  by  Figs.  3.5,  3.6,  and  3.7.  The  most 
obvious  implication  of  this  is  that  some  of  the  target’s  component  scatterers  have  amplitude  patterns 
that  vary  significantly  with  aspect.  Visual  examination  of  ship  structures,  for  example,  finds  obvious 
flat  plates,  dihedrals,  trihedrals,  etc.  The  RCSs  of  such  generic  objects  have  been  calculated  theoreti¬ 
cally  and  are  well  known.  The  flat  plate  RCS,  for  example,  is  strongly  dependent  on  plate  orientation; 
it  varies  dramatically  near  vertical  incidence  where  a  very  large  RCS  relative  to  its  physical  size  is  gen¬ 
erated  [70], 

Mitchell  [23]  has  suggested  that  there  are  only  four  basic  scattering  mechanisms  within  a  radar 
resolution  cell.  To  describe  these  mechanisms,  he  defines  four  basic  scattering  elements;  he  claims 
that  for  simulation  purposes,  they  describe  all  scattering  effects.  The  basic  elements  are: 

1.  rough  surface  (diffuse  scattering), 

2.  point  source  (sphere  and  dihedral  corner), 

3.  flat  plate  (specular  scattering),  and 

4.  line  source  (edge,  dihedral  corner,  thin  cylinder,  and  wire). 
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In  view  of  this,  we  propose  using  a  small  set  of  generic  functions  to  model  the  amplitude  functions  of 
the  USs.  This  set  may  vary  depending  on  the  type  of  target  being  modeled.  For  example,  ship  models 
would  tend  to  include  flat  sided  scatterer  types  (plates,  dihedrals,  dihedrals,  etc.)  whereas  aircraft 
models  would  tend  to  include  curved  surface  scatterer  types  (ogive,  parabaloid,  etc.). 

In  the  simulations  (Section  9),  we  model  the  US  amplitude  functions  using  only  three  generic 
scattering  types:  sphere,  flat  plate,  and  corner  reflector.  These  types  were  chosen  because  they 
represent,  to  first  order,  the  most  obvious  scattering  effects  and  they  are  consistent  with  the  fundamen¬ 
tal  set  proposed  by  Mitchell.  The  fidelity  of  the  model  by  using  only  these  types  of  scatterers  must  be 
determined  relative  to  measured  data  which  we  do  not  have  access  to  for  this  report. 


The  sphere  type  was  chosen  to  represent  scattering  that  is  specular  but  is  also  significant  (in 
amplitude)  over  a  wide  angle  of  illumination.  We  model  the  amplitude  of  such  a  scatterer  as  being  pro¬ 
portional  to  the  projected  area  of  the  object  as  viewed  by  the  radar,  as  it  is  (in  the  optical  limit)  for  a 
sphere  [33]. 


To  represent  scattering  that  is  specular  and  highly  directional,  in  analogy  with  the  equation  for  the 
RCS  of  a  flat  plate  near  specular  incidence  we  suggest  the  representation 


A  = 


sin‘(27ra  sin  0) 
2na  sin  0 


sin  (2 nb  sin  <ft) 
2nb  sin  <f> 


1/2 


(3.22) 


where  a  and  b  are  the  plate  dimensions  and  0  and  tf>  are  azimuth  and  elevation  angles  between  the 
plate  and  the  radar-to-US  sight  line  and  a  line  perpendicular  to  the  flat  surface.  We  anticipate  that  this 
type  of  US  will  be  realized  by  actual  flat  plates  which  for  a  ship  appear  primarily  at  the  bow,  stern,  and 
broadside  aspects.  Unlike  the  sphere  type  of  US,  this  type  is  frequency  dependent. 


The  corner-reflector  type  was  chosen  to  represent  frequency-dependent  scattering  that  is  nonspec- 
ular.  This  type  of  scattering  results  from  multiple  reflections  [61].  We  model  the  associated  amplitude 
by  the  ratio  of  the  US’s  projected  area  to  the  radar  wavelength.  This  is  analogous  to  the  equation  for 
the  maximum  RCS  of  a  corner  reflector. 


_  4na 2 

°"max  =  2 

A 


(3.23) 


where  a  is  the  area  of  the  reflector’s  aperture  [61]. 


3.5  Phase 

For  a  point-source  scatterer  located  at  some  position  R  relative  to  the  radar,  the  phase  of  the  scat¬ 
tered  field  at  the  radar  results  from  the  properties  of  and  the  range  to  the  scatterer.  For  a  sphere  the 
phase  varies  as  the  range  to  the  sphere’s  center  varies  but  is  constant  with  sphere  orientation.  How¬ 
ever,  for  nonspherical  objects,  the  phase  depends  not  only  on  some  location  point  for  the  object  (say, 
its  geometrical  center)  but  also  on  the  object’s  orientation.  This  is  because  as  the  object  rotates, 
scattering  from  components  of  the  object  (edges,  protrusions,  etc.)  vary  as  a  function  of  aspect.  This 
apparent  point-source  location  of  the  object  is  called  the  object’s  scattering  center.  Aggregating  several 
physical  objects  may  create  a  US  of  substantial  physical  size  relative  to  the  radar  wave  length.  As  such 
a  scatterer  rotates,  the  scattering-center  location  may  vary  rapidly  for  the  wave  lengths  of  interest  in  this 
report  (on  the  order  of  a  few  centimeters).  Therefore,  the  US  phases  (the  <6y  in  Eq.  (3.8a))  are  defined 
as  the  apparent  location  of  the  US’s  scattering  center  relative  to  some  fixed  point  on  the  US,  as  a  func¬ 
tion  of  orientation. 

The  scattering  center  will  also  depend  on  the  type  of  ranging  technique  employed  by  the  radar. 
For  example,  narrow  pulse,  leading-edge  tracking  will  yield  different  apparent  ranges  for  scattering 
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objects  than  will  peak-signal  techniques.  However,  in  our  development  we  do  not  attempt  to  explicitly 
include  ranging-technique  effects  primarily  because  of  the  lack  of  existing  data  describing  US  phase- 
variations  as  a  function  of  orientation. 

We  suggest  that  there  are  two  dominant  physical  causes  (excluding  multipath  effects)  of  variations 
in  the  US  phase  for  small  changes  in  aspect  angle.  First,  we  expect  the  most  significant  phase  varia¬ 
tions  to  be  caused  by  the  motion  of  the  surfaces  that  yield  the  USs;  that  is,  the  relative  motion  of  the 
geometric  center.  Second,  we  expect  that  US  rotation  will  also  introduce  phase  variations.  We  expect 
two  types  of  effects  due  to  US  rotation:  low  frequency  and  high  frequency.  The  low-frequency  effects 
should  result  from  the  objects  that  generate  the  US  moving  closer  or  farther  away  from  the  radar,  thus 
producing  a  change  of  phase  (due  to  earlier  or  later  pulse  arrivals).  The  high-frequency  variations 
should  result  for  USs  which  arise  from  the  relative  phasing  of  various  objects.  As  the  objects  rotate, 
their  relative  phasing  may  change  producing  glint  [Section  8],  which  is  a  high-frequency  phenomenon 
[Section  9].  However,  we  expect  this  effect  to  be  small  because  glint  is  inversely  proportional  to  RCS 
and  the  existence  of  the  US  implies  a  large  RCS. 


We  will  model  the  US  phase  by 

2|J?|  +  0(0)  +  f(co,  l) 

where 


<£(/?,  0,  ai,  t)  = 


2n 

X 


(3.24) 


R  is  the  location  of  the  geometric  center  relative  to  the  radar, 

0  is  the  orientation  of  the  US  relative  to  the  radar, 

X  is  the  radar’s  wavelength, 

0()  represents  the  low-frequency  orientation  effects,  and 
f  (o>,  t)  represents  the  high-frequency  effects. 

Because  of  the  causal  relationship  between  US  motion  and  target  motion,  the  £(&>,  t)  process  should  be 
correlated,  to  some  extent,  with  the  target  motion.  The  variance  of  £(w,  t)  most  likely  will  be  small 
relative  to  2tt,  otherwise  f  (w,  t)  could  be  approximated  by  a  phase  process  that  is  uniformly  distributed 
over  (0,  2tt ) .  This  would  yield  a  target  RCS  that  is  Rayleigh  distributed  (see  Section  7)  which  is  in 
general  not  true.  This  observation  is  consistent  with  our  comments  in  the  previous  paragraph. 

Because  of  the  lack  of  experimental  data  and  the  inlractibility  of  US-phase  analysis,  we  will  not 
attempt  to  model  or  simulate  the  £(w,  i)  process  by  other  than  uniformly  distributed,  uncorrelated 
processes.  The  interval  over  which  the  processes  are  defined  are  specified,  in  a  simulation  input,  as  a 
percentage  of  the  radar  wavelength.  Simulation  results  have  indicated  that  if  the  random  phase  varia¬ 
tions  on  a  pulse-by-pulse  basis  exceed  about  5%  of  the  radar  wavelength,  then  they  become  a  signifi¬ 
cant  cause  of  variations  in  the  scattered  field,  which  intuition  suggests  should  not  happen. 

Finally,  we  expect  that  for  small  changes  in  aspect  angle,  scattered-field  variations  are  due  pri¬ 
marily  to  phase  variations  rather  than  amplitude  variations. 


4.  MODEL  OF  THE  RECEIVED  SIGNAL 


Although  a  radar  measures  the  EM  field  at  its  antenna  aperture,  evaluation  of  a  radar's  perfor¬ 
mance  can  be  made  knowing  only  a  few  projections  of  that  field.  The  most  useful  projections  are 
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received  power  (RCS)  and  the  azimuth  and  elevation  error  signals  (glint).  The  validity  of  any  model  of 
the  scattered  field  that  is  to  have  practical  applications  must  be  measured  relative  to  such  projections. 
In  this  chapter  we  formulate  the  received  signal  resulting  from  scattering  by  an  N- source  target  that  is 
modeled  by  unit-scatterers.  We  develop  representations  for  the  received  field  and  the  associated 
azimuth  and  elevation  error  signals.  The  received  field  representation  is  used  in  later  chapters. 

4.1  Scattered  Field 

Because  we  are  modeling  the  target  as  a  finite  number  of  USs,  the  received  electric  field  at  the 
radar  is  the  linear  superposition  of  the  fields  from  each  US.  We  are  interested  in  representing  only  the 
pulse-to-pulse  variations  in  the  received  signal,  not  intrapulse  properties,  so  we  assume  that  the 
received  pulses  are  ideal  in  the  sense  that  they  are  completely  defined  by  an  amplitude,  phase,  and 
pulse-width.  Further,  we  assume  that  the  target  extent  is  small  relative  to  the  radar-to-target  range  and 
that  the  radar  pulse  width  is  large  relative  to  target  extent.  This  allows  the  assumption  that  the  return 
pulses  are  time  coincident.  This  time-coincident  assumption  would  not  be  acceptable  if  we  were  analyz¬ 
ing  range  tracking  errors  or  the  structure  of  the  received  pulse’s  leading  or  trailing  edge. 

Assuming  uniform  target-illumination  (due  to  large  radar-to-target  range)  and  identical 
transmitter  and  receiver  polarizations,  the  received  field  at  the  radar  due  to  scattering  by  a  point-source 
target  is  usually  represented  by  an  equation  of  the  form 

£(£,,?„  f)  -  F2(R, )g(R,) a, (9, ) (4.1) 


where 


w  is  the  radar  frequency, 

Rt  defines  the  location  of  the  scatterer  relative  to  the  radar, 

9,  defines  the  orientation  of  the  scatterer  relative  to  the  radar, 

g()  describes  the  antenna  and  propagation  effects  on  the  signal, 

F(-)  describes  the  effects  of  multipath, 

a,  is  the  amplitude  of  the  scattered  field,  and 

</>,  (•,  •)  is  the  phase  angle  of  the  scattered  field. 

Three  transmit  path  effects  are  accounted  for  here:  phases  due  to  radar-to-scatterer  path  lengths, 
attenuation  due  to  path  losses,  and  variations  in  the  illuminating  signal  due  to  multipath.  The  phases 
are  accounted  for  in  t^((\  ■)  and  the  loss  effects  in  g().  The  multipath  effects  are  discussed  in  detail  in 
Sec.  4.4,  so  for  the  remainder  of  this  section  we  assume  that  F2(  )  =  1  (i.e.,  no  multipath  effects). 

Summing  the  individual  fields  from  an  Af-source  target,  the  received  field  becomes 

£(/)-  £  g  ( /?, ) a,  (?, )eJ*'{*' ' * eJu ' .  (4.2) 

i-i 

The  assumption  of  long  range  implies  that 

g(Rt)  =  g(Rj ),  for  all  /,  j  (4.3) 
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so  that  we  can  approximate  Eit)  by 

EU)  -  g(R)eJul  £  at( 


(4.4) 


i-i 


where  R  locates  a  fixed  point  on  the  target;  for  example,  the  target’s  center-of-gravity.  Assuming  that 
the  N  sources  of  the  target  are  USs,  we  see  from  Eq.  (3.11)  that 


E(t)  =  g(R )<?>'  T ipT,  pR )  -  g(R)e*‘pI 


I  X,(R„  9,) 

i-i 


Pr 


(4.5) 


where  pT  and  pR  define  the  polarization  that  is  implicit  in  Eq.  (4.1).  If  we  make  the  usual  assumption 
that  pR  —  pR  =  (1,  0) r,  then  Eq.  (4.5)  becomes 


EU)  -  g(R)e'“'  '£AUi(0l)e  ' 


j<tll(Rrel) 


(4.6) 


<- 1 


We  next  discuss  the  effect  of  the  radar  system  on  the  scattered  field. 


4.2  Radar  System  Model 

Although  our  primary  concern  is  target-scattering  modeling  and  not  radar-system  modeling,  prac¬ 
tical  validation  of  any  scattering  model  will  of  necessity  involve  actual  target  measurements  made  by 
specific  radars.  Further,  we  are  concerned  with  scattering  models  from  the  point  of  view  of  analyzing 
the  effects  of  target  scattering  on  tracking  radars  rather  than  that  of  predicting  scattered  electric  fields. 
Therefore  for  completeness,  we  include  a  representation  of  tracking  radars  that  can  be  used  with  the 
US-model  formulation. 


We  view  tracking  radars  as  being  composed  of  three  basic  functional  units:  a  "front-end,"  error 
detectors,  and  feedback  loops.  Figure  4.1  is  a  block  diagram  of  a  generic  amplitude-comparison,  mono¬ 
pulse  radar  with  each  functional  unit  indicated.  Such  a  radar  is  commonly  used  as  a  trackin'*  at  ar  at 
microwave  frequencies.  Because  this  type  of  radar  is  widely  used  in  practice,  we  limit  our  model  of 
radars  to  this  type.  This  limitation  does  not  significantly  restrict  the  results  to  be  obtained  because  the 
primary  influence  on  tracking  performance  is  the  fluctuation  of  the  field  scattered  from  the  target. 
Further,  phase-comparison  monopulse  has  been  shown  to  be  functionally  equivalent  to  amplitude- 
comparison  monopulse  [71],  and  other  forms  of  tracking  radars  (e.g.,  sequential  lobing,  conical  scan, 
etc.)  differ  primarily  in  the  form  of  the  front-end  transformation  and  error  detector  implementations. 
In  any  event,  the  US  concept  is  independent  of  specific  radar  implementations. 

4.2.1  Front-End 


We  first  consider  the  effects  of  the  radar  front-end:  the  components  that  transform  the  input  EM 
field  (at  the  antenna’s  aperture)  to  a  complex  signal  at  some  intermediate  frequency  (IF).  This  IF  sig¬ 
nal  is  then  processed  to  obtain  tracking  information.  Our  approach  is  to  model  the  components  that 
convert  the  input  field  to  an  IF  signal  as  a  simple,  memoryless  transformation,  which  we  represent  by 

VU)-H[EU)  1  (4.7) 

where  V(t)  is  a  voltage.  In  general,  the  front-end  transformation,  //(•),  is  very  difficult  to  specify  pre¬ 
cisely  because  it  includes  the  effects  of  the  antenna,  waveguide,  hybrids,  mixers,  and  IF  amplifiers. 
Clearly,  Hi)  depends  strongly  on  the  specific  radar  implementation,  particularly  the  antenna  used.  For 
the  applications  of  interest  here  (where  target-radar  motion  is  causing  the  major  fluctuations  in  the 
received  signal),  these  devices  can  be  accurately  characterized  by  simple  models:  hybrids  with  input 
signals  A  and  B  as  producing  outputs  A  +  B  and/or  A-B,  mixers  as  transforming  signals  with  spectra 
Si co)  im  >  0)  to  signals  with  spectra  Sim  -  m„)  im  -  mm  >  0)  where  m„  is  the  mixer  frequency,  and 
IF  amplifiers  as  constant  gain  (or  logarithmic)  devices  with  feedback  gain  control.  An  overall  power 
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Fig.  4.1  —  A  two-coordinate  (azimuth  and  elevation),  amplitude-comparison, 
monopulse  tracking  radar  (adapted  from  Ref.  68) 


loss,  due  mainly  to  the  hybrids,  and  an  overall  phase  shift,  due  mainly  to  physical  lengths  and  the 
mixers,  complete  the  description  of  the  guided  signals.  The  video  signal  is  produced  by  detectors 
(amplitude,  phase,  or  sum-and-difference  [711)  that  transform  the  guided  fields  to  voltages.  The  most 
complex  aspect  of  //(•)  is  the  transformation  of  an  EM  field  across  the  antenna  aperture  to  an  EM  field 
within  an  associated  waveguide.  This  transformation  is  typically  represented  by  the  antenna  amplitude- 
gain  pattern  (to  be  referred  to  as  the  antenna  pattern). 

Since  we  assume  that  the  scatterers  are  point-source  radiators  which  are  at  long  range,  it  will  be 
sufficient  (as  has  been  the  case  in  many  other  studies;  for  example,  see  Refs.  72  and  73)  to  model  the 
front-end  as  the  antenna  pattern  coupled  with  a  complex  gain  factor  and  a  frequency  shift.  Specifically, 
if  the  antenna  amplitude-gain  (complex  in  general)  is  given  by  G(0O,  /3f)  where  /9a  and  /3,  represent 
azimuth  and  elevation  angles  relative  to  the  radar  boresight,  then  for  a  point-source  target  whose  angu¬ 
lar  coordinates  relative  to  the  radar  boresight  are  given  by  ifia,  f}t),  we  model  the  front-end  transfor¬ 
mation  as 

V(t)  -  H[E(t)]  -  KG(fia,  P')E0(fla,  £„  t)  (4.8) 

where  K  is  complex  and  represents  the  nonantenna  front-end  gains  and  phase  shift,  and  £0(/3o,  0r,  t) 
is  the  value  of  the  scattered  electric  field  at  the  aperture  location  defined  by  /3a  and  /9,.  The  carrier  fre¬ 
quency  of  V(t)  is  the  frequency  of  E0(fia,  0, ,  t)  lowered  by  the  mixer  frequency. 

Because  we  are  assuming  tracking,  we  are  only  interested  in  antenna  characteristics  near  boresight, 
say  ±3  dB  about  the  main  beam  maximum.  Therefore  it  may  be  useful  to  assume  an  ideal  approxima¬ 
tion  for  the  antenna  pattern  that  has  a  simple  analytical  form  yet  accurately  represents  typical  antenna 
amplitude-gain  characteristics  in  this  restricted  region: 

G(9 ,  0)  -  G0  -  k(9  -  90 )2  -  fr(0  -  0O)2  (4.9) 

where  0O  and  0O  define  the  pattern  maximum  relative  to  a  given  radar  axis  (say,  boresight)  and  G0  and 
k  are  constants.  However,  we  note  that  the  long  range  assumption  will  often  allow  the  antenna  pattern 
to  be  modeled  as  a  constant  gain  because  of  the  resulting  negligible  separations  of  the  unit-scatterers. 
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4.2.2  Error  Detectors 

The  error  detectors  transform  guided  EM  waves  to  video  signals  useful  for  determining  antenna 
pointing  (angle)  errors.  Rhodes  171]  has  shown  that  in  monopulse  radars  there  are  only  three  distinct 
kinds  of  angle  detection:  amplitude,  phase,  and  sum-and-difference.  Because  any  monopulse  error  law 
can  be  described  as  a  function  of  the  difference-to-sum  ratio  (A/Z)  [74],  we  model  the  detectors  as 
producing  sum  and  difference  signals  which  yield  as  output  A/Z.  Assuming  instantaneous  automatic 
gain  control  (I AGC),  the  outputs  of  the  phase-sensitive  detectors  of  Fig.  4.1  are  the  real  part  of  A/Z 
due  to  the  sum  signal  being  fed  back  to  the  various  IF  amplifiers  [71], 

4.2.3  Feedback 

There  are  two  feedback  loops  of  interest:  IF  amplifier  feedback  and  antenna  servomechanism 
(servo)  feedback.  The  amplifier  feedback  is  simply  a  feedback  of  the  sum  signal  to  control  the  IF 
amplifier  gains. 

Angle  error  signals  that  are  output  from  the  detectors  are  used  to  control  the  pointing  direction  of 
the  antenna.  This  is  done  by  feedback  of  these  signals  to  servos  that  control  the  position  of  the 
antenna.  Type  II  servos  are  usually  employed  for  this  purpose  [68].  We  note  that  the  long  range 
assumption  (implying  small  angular  separations  of  the  scatterers)  will  often  allow  the  effects  of  the  ser¬ 
vos  to  be  ignored. 

4.3  Received  Signals 

We  now  describe  the  received  signal  in  terms  of  the  scattered  field.  First  we  define,  in  earth-  or 
radar-coordinates,  the  radar  boresight  direction  pBS  =  (/3 a,Pe),  and  the  radar-to-/th  scatterer  sight 
line, 

P,  ~  Pe)-  The  antenna  gain  associated  with  the  ith  scatterer  is  then  given  by  Gifi,  -  (iBS )  where 

G0i  ~  Pbs )  =  Gifiat  -  /3a ,  (3e.  -  (ie).  (4.10) 

With  this  definition  each  IF  signal  has  the  form  (from  Eqs.  (4.8)  and  (4.5)) 

VU)  -  £(£)e7(“"<‘’",'  £  KG0,  -  p BS)pfXl(Ki,9i)pK.  (4.11) 

;-l 

More  specifically,  there  are  three  signals  of  interest:  one  sum  and  two  difference  signals.  Each 
signal  has  an  associated  sum  and  difference  pattern  (Gi(/3),  GAo(j 3),  and  GAf(/3))  and  complex,  front- 
end  gains  (Kz,  Kita,  and  K 4  ).  Suppressing  the  first  two  terms  of  Eq.  (4.11),  the  sum,  azimuth- 
difference,  and  elevation-difference  signals  become 

Vi (t)  -  Kz  £  Gr( P,  -  pBS)pIX,(Rn  9,)pR ,  (4.12a) 

/-i 

FA(1(f)  =  A4a  J  —  PBs)f>T  XjiRj,  9,)pB ,  and  (4.12b) 

/-I 

yA,(t)  “  Kif  £  GAe0,  -  ?«)#*,(*„  ®r)Ar-  (4.12c) 

i- 1 

The  three  antenna  patterns  are  derived  from  Eq.  (4,9)  where  we  have  assumed  that  each  of  the 
four  radiators  that  comprise  the  antenna  have  identical  patterns  and  that  each  is  offset  relative  to  the 
radar  boresight.  The  sum  pattern  is  represented  by 
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Ci0,>  -  Ci(/S.,.  0.)  -£GW.,  -  0.,.  0,,  -  0.1  -  f  G<0,  -  0j)  (4.13a) 

J- 1  j- 1 

where  /8y  —  fie)  represents  the  associated  radiator’s  offset  relative  to  the  radar  boresight.  The 
difference  patterns  are  represented  by 

Gia(i 3,)  -G0-  £,)  +  G0-  p})  -  G0-  p2)  ~  G0-  fi4)  (4.13b) 

and 

G±t0i)  -  G0-px)  +  G0-p2)  -  G0-pi)  -  G0-pA).  (4.13c) 


Equation  (4.11)  is  our  fundamental  representation  of  the  monopulse-radar  received-signal  due  to 
N  point-sources.  Representations  of  this  for  radar-system  have  been  used  for  analyses  by  others;  the 
most  common  use  assumes  a  mix  of  independent  processes  and  deterministic  functions  for  one  or 
several  of  the  target  parameters  (e.g..  Refs.  11,  19,  20,  23,  24,  25,  75,  and  76). 

4.4  Sea-Surface  Multipath  Effects 

If  the  target  is  a  surface  ship,  then  for  small  and  moderate  radar-to-ship  grazing  angles,  sea- 
surface  multipath  affects  the  received  signals.  These  multipath  effects  alter  both  the  signal  illuminating 
the  ship  and  the  scattered  field. 

As  discussed  in  the  introduction,  the  effects  of  multipath  on  the  received  power  of  a  signal 
transmitted  by  a  point  source  is  well  understood,  but  an  apparent-point-source  representation  of  the 
scattered  field  has  only  been  determined  for  the  specular  component.  However,  because  the  diffuse 
scattering  arises  from  scattering  over  a  large  area,  there  are  at  least  several  scatterers  comprising  the 
ship  at  any  aspect,  and  because  we  are  interested  in  only  a  representation  of  "average"  seas  (determined 
only  by  a  sea-state  parameter),  we  will  model  the  diffuse  vector’s  location  point  as  the  associated 
target’s  image-point.  Although  we  will  develop  glint  equations  using  this  model,  clearly  it  is  not  restric¬ 
tive  and  the  method  applies  for  other  diffuse  models  as  well. 

To  be  consistent  with  the  multipath  models  employed,  we  assume  that  the  radar’s  antenna  has  a 
beamwidth  wide  enough  to  fully  illuminate  the  first  several  Fresnel  zones  about  the  specular  point 
(defined  in  Fig.  4.2).  We  note  that  this  assumption  will  hold  even  for  relatively  narrow  beams  if  the 
radar-to-target  range  is  large.  Because  of  our  assumption  of  low  grazing  angles  and  the  reasons  stated 
above,  we  also  ignore  time  delays  associated  with  the  multipath  signals,  noting  that  when  analyzing 
range-tracking  errors  such  delays  should  be  accounted  for. 


Fig.  4.2  —  Multipath  geometry 
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Beard  et  al.  [46]  show  that  if  the  received  field  is  D  without  multipath  effects,  then  with  mul¬ 
tipath  effects  the  received  field  becomes 

T  -  D  +  R  (4.14) 

where  the  total  received  field,  T,  is  the  sum  of  the  direct  field,  D,  and  a  reflected  field,  R.  The 
reflected  field  appears  to  the  receiver  to  be  transmitted  from  a  point  located  at  the  transmitter  image 
reflected  about  the  plane  of  the  sea  surface  (Fig.  4.2).  The  reflected  field  is  modeled  as  the  sum  of  two 
components,  one  that  is  deterministic  and  coherent,  relative  to  the  direct  field,  and  one  that  is  stochas¬ 
tic  and  incoherent, 

R-C  +  7.  (4.15) 


The  coherent  and  incoherent  fields  are  also  referred  to  as  the  specular  and  diffuse  fields.  For  conveni¬ 
ence  Eq.  (4.14)  is  usually  normalized  relative  to  the  direct  field  D  by 


FD, 

(4.16a) 

where  F  is  a  multipath  coefficient  defined  by 

F  “  1  +  Pc  +  Pi 

(4.16b) 

where 

Pc-Pseja, 

(4.16c) 

and 

Pi  ”  Pp  +  jpQ • 

(4.16d) 

The  phase  angle  a  is  due  to  path  length  differences  and  sea-surface  reflection  characteristics.  The  in- 
phase  and  quadrature  terms,  pP  and  pQ,  are  zero-mean  Gaussian  processes.  Curves  defining  the  vari¬ 
ance  of  the  p/  process  (assuming  that  pP  and  pe  have  equal  variances)  are  presented  in  Ref.  49.  To 
account  for  the  fact  that  the  scattering  coefficients  pc  and  p,  depend  on  geometry  we  use  the  notation 

F  -  F(.R,  0),  (4. 1 6d) 

and  for  the  ith  US  we  abbreviate  Eq.  (4.16d)  by 

Fj  -  F(Riy  ?,).  (4.16e) 

Because  the  sea  surface,  radar,  and  target  effectively  do  not  move  during  the  time  of  pulse 
transmission  and  reception  at  the  radar  (for  the  problems  of  interest  here),  the  transmit  and  receive 
modes  are  reciprocal.  Therefore  the  transmit  and  receive  multipath  effects  are  identical  except  for  the 
effects  of  the  transmit  and  receive  radiation  patterns  of  the  antenna  and  target. 

Because  the  ship  is  close  to  the  scattering  surface,  the  multipath  image  locations  are  close  to  their 
associated  scatterer  locations.  Note  that  if  the  antenna  pattern  and  geometry  are  such  that  G0,  -  PBS) 
=  G(J3j  -  fJB$)  f°r  '  and  Ji  then  the  received  signal  can  be  approximated  by 

VU)  -  g(R)eJ(a~a‘m)'KG0  -  0BS)  £  F2{Ri,'9i)plXi{RlJi)pR.  (4.17) 

1 

However,  if  the  antenna  effects  on  the  direct  vs  the  reflected  signals  are  to  be  accounted  for,  we  must 
view  the  reflected  signal  as  arising  from  a  source  located  below  the  surface  of  the  sea.  As  is  often  done 
[49]  we  model  the  source-location  of  the  reflected  signal  as  the  apparent  location  of  the  specular  reflec¬ 
tion  (Fig.  4.2). 

To  modify  the  equivalent-point-source  multipath  coefficient  Eq.  (4.16a)  to  become  an  jV-source 
model  (i V  —  2),  we  view  the  problem  as  follows  (where  for  clarity  we  assume  that  the  propagation  and 
scatterer  reflection  effects  are  accounted  for  in  g(  )  and  X(-,-)).  Assume  that  the  incident  field  at  the 


26 


NRL  REPORT  8887 


scatterer  has  the  form  £  —  FE0  where  £  is  the  direct  Held  (no  multipath).  Then  if  we  view  the 
scatterer  as  transmitting  the  signal  ££0 ,  then  the  signal  received  at  the  radar  will  be  £2£0.  From  the 
point  of  view  of  the  radar,  the  received  signal  appears  to  have  a  direct  component  D  =  FE0  and  a 
reflected  component  R  **  (nc  +  p,)££0.  Therefore  we  separate  the  two-way  multipath  coefficient  F 2 
for  the  ith  US  into  direct  and  reflected  coefficients: 

£,  -  F2E0  -  (1  +  pc  +  Pi)F,E0 


-  F,E0  +  iF,  -  1  )£(£0-  A  +  R,. 

Finally,  we  need  the  direction  of  propagation  for  D,  and  Rt. 

(4.18) 

If  in  a  given  coordinate  system  /?,  locates  the  /th  scatterer  and  Rr  locates  the  receiver,  then  the 
direction  of  propagation  for  the  /th  direct  signal  is  given  by 

Rd,  “  Rr  ~  Ri- 

In  Earth-axis,  cartesian  coordinates, 

(4.19) 

Rd,  —  (RXr  -  Rx>  RyR  —  Ry.,  hR  —  h,) 

(4.20) 

where  hR  and  hi  are  the  heights  of  the  radar  and  the  ith  source  above  Earth, 
of  propagation  for  the  direct  and  reflected  field  components  are 

Rd, 

Therefore  the  direction 

*  1 
'°r  i*0,i 

and 

Rr, 

(4.21a) 

r*‘~  1^1 

where 

(4.21b) 

Rr,  A  (RXr  -  Rxr  RYr  -  Ry,,  hR  +  h,).  (4.22) 


With  the  separation  of  Eqs.  (4.18)  and  (4.21),  we  see  that  by  viewing  the  reflected  fields  as  an  addi- 
lional  set  of  ,Y  scalterers,  Eq.  (4.11)  becomes 

V(i)  =»  K(R)eJ'u-uJ'K  £  [G0i  -  Pbs)F,  +  G0,+N  -  )£,(£,  -  1  )\p}X,pl  (4  23) 

i-i 


Though  we  use  this  point-source  model  to  account  for  the  effects  of  multipath,  USs  result  from 
surfaces  that  are  distributed  in  space.  When  there  is  no  multipath  we  can  represent  the  scattered  field 
from  a  US  by 

Es  =  f  —  In  x  (}  x  e)]  x  ?e~2jk'r  T’dS'.  (3.16) 

2 rrR  J  5  -q 

Using  the  point-source  multipath  model  to  represent  the  incident  field,  we  represent  the  scattered  field 
in  the  presence  of  multipath  by 


£' 

*^mp 


EqF2(7  ') 
V 


In  x  (?  x  e)]  x  re~2'k}  7  dS' 


(4.24) 
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where  F(r')  represents  the  complex  multipath  factor  at  the  point  r  If  over  S', 

F(T)  =  F0,  (4.25) 

then  we  have  the  point-source  representation  described  above: 

El,  =  FlE\  (4.26) 

Although  data  collected  by  NRL  and  analyzed  by  the  author  indicate  that  FiT)  can  often  be  approxi¬ 
mated  by  an  azimuth-independent  process  over  relatively  large  variations  in  azimuth,  other  data^have 
shown  [77]  that  F(7)  rapidly  decorrelates  in  elevation.  Therefore  for  most  USs,  representing  El,  by 
Eq.  (4.26),  with  Fq  —  F( 7?)  where  7]  is  the  /th  US  location,  will  not  be  valid.  The  difficulty  raised 
here  is,  to  our  knowledge,  an  open  question:  can  the  target  effects  and  the  multipath  effects  be 
represented  in  Eq.  (4.24)  in  a  way  that  yields  an  equivalent-point-source  model  for  the  scattering? 


As  stated  above,  the  point-source  multipath  model  is  valid  when  the  transmitting  and  receiving 
amplitude  patterns  (gain  pattern  for  an  antenna)  are  not  narrow.  This  constraint  will  likely  hold  for  the 
sphere  type  and  corner-reflector  type  of  US  amplitude  functions.  For  the  highly  specular  fiat-plate  type 
a  modification  to  the  multipath  effects  may  be  required  analogous  to  the  antenna  gain-pattern  correc¬ 
tion  suggested  in  Ref.  49.  Such  a  modification  was  not  made  in  the  simulation. 

Because  multipath  is  but  one  parameter  (though  an  important  one)  affecting  the  ship  scattering 
and  because  the  point-source  model  presented  here  is  the  only  validated  stochastic  model  of  over-water 
multipath  that  we  are  aware  of,  we  have  used  it  in  our  simulation  studies.  In  Section  9  we  discuss 
results  obtained  using  this  modeling  approach  in  the  simulation. 

5.  UNIT-SCATTERER-MOTION  MODELS 

In  this  section  we  develop  equations  of  motion  for  an  arbitrary  US  location  because  this  motion 
determines  the  amplitude  and  phase  variations  in  the  received  field.  Our  purpose  is  not  ship-motion 
modeling  per  se  but  US-motion  modeling  for  use  in  analyzing  and  simulating  scattered  fields.  How¬ 
ever,  for  completeness  we  investigate  two  basic  approaches  to  ship-motion  modeling.  In  Sec.  5.2  we 
formulate  a  motion  model  that  assumes  knowledge  of  the  sea  forces  that  cause  ship  motion.  Then  in 
Sec.  5.3,  we  discuss  a  frequency-response  model  that  has  been  quite  useful  in  ship-motion  modeling 
given  knowledge  of  the  sea  spectrum.  The  DTNSRDC  has  used  the  model  to  generate  libraries  of 
center-of-gravity-motion  (cg-motion)  spectra  for  various  ships  (e  g..  Refs.  78  and  79).  In  later  chapters 
we  begin  our  analyses  and  simulations  by  assuming  knowledge  of  these  cg-motion  spectra.  Because 
these  motions  are  zero-mean  Gaussian,  the  associated  spectra  are  all  that  is  required  to  completely 
specify  them  as  stochastic  processes. 

Our  reliance  on  the  DTNSRDC  spectra  is  not  without  its  costs.  Because  of  a  strong  desire  to 
develop  models  that  are  intuitively  appealing  and  easy  to  work  with  (assuming  a  minimal  background  in 
mathematics),  the  sea-surface  model  that  has  been  developed  by  the  oceanographic  community  is  in  a 
form  that  is  nonstandard  relative  to  the  system-theoretic  formulation.  The  differences  are  slight  but 
confusing  and  can  lead  to  errors  if  not  correctly  accounted  for:  the  spectrum  is  defined  for  positive  fre¬ 
quencies  only  and  so  that  the  integral  of  the  spectrum  is  the  energy  in  the  sea  surface  (at  a  point  on  the 
surface).  To  minimize  the  confusion  that  would  result  when  the  tabulated  spectra  are  used  in  a  simula¬ 
tion,  we  have  chosen  to  follow  the  nonstandard,  oceanographic  spectrum-formulation.  This  allows  us 
to  directly  use  the  tabulated  (nonstandard)  spectra. 

5.1  Approaches  to  Linearization 

We  use  linear  equations  because  they  are  widely  used  in  the  study  of  both  ship  and  aircraft 
motion  and  are  easy  to  work  with.  We  are  interested  primarily  in  small  deviations  in  the  target’s 
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motion  from  a  nominal  trajectory,  and  also  in  target  motion  modeling  only  over  relatively  short  periods 
of  time  (usually  on  the  order  of  seconds  to  no  more  than  several  minutes).  There  are  two  basic 
approaches  used  in  developing  linear  equations  of  ship  motion.  The  first  approach  is  to  transform 
Newton’s  six  nonlinear  force  and  moment  equations  from  the  target’s  eg  axis-system  to  an  axis  system 
that  is  determined  by  target  symmetry.  The  equations  are  then  linearized  in  the  symmetry  axis-system. 
We  call  this  the  transform-first  approach.  The  second  approach  is  to  assume  linear,  second-order  dif¬ 
ferential  equations  in  the  eg  axis-system  to  describe  the  six  forces  and  moments.  The  various  coeffi¬ 
cients  are  then  estimated  via  calculations  or  measurements.  We  call  this  the  linearize-first  approach. 
The  choice  of  approach  usually  depends  on  the  ease  of  obtaining  values  for  the  equation  coefficients 
and  the  feasibility  of  measuring  the  modeled  forces  and  moments  for  the  target  of  interest. 

The  transform-first  approach  to  linearization  is  strongly  dependent  on  target  structure,  especially 
symmetry.  Using  this  approach,  various  sets  of  linear  equations  have  been  developed  for  ships  and  air¬ 
craft.  The  form  of  these  equation  sets  varies  depending  on  intended  applications. 

The  linearize-first  approach  begins  with  a  set  of  linear,  second-order  differential  equations  for  the 
ship’s  cg-motion.  Therefore  assuming  rigid  body  motion,  the  equations  of  motion  for  an  individual  US 
have  the  same  form  regardless  of  target  structure.  We  will  use  the  linearize-first  approach  in  develop¬ 
ing  ship  motion  models  that  are  useful  in  ship  scattering  studies. 

5.2  Motion  Equations  Given  the  Driving  Forces  and  Moments 

We  begin  by  deriving  a  set  of  motion  equations  for  an  arbitrary  US  location  by  assuming 
knowledge  of  linearized  equations  of  motion  for  the  target’s  eg  and  the  driving  forces  and  moments. 
We  also  assume  that  the  target  undergoes  only  rigid-body  motion  and  that  the  six  motions  of  the  target 
are  uncoupled  (i.e.,  u  «■  x,  v  «  y,  . . . ,  r  —  d»). 

5.2.1  Definition  of  the  Axis  Systems 

We  use  two  axis-systems  in  describing  US  motion:  an  Earth-axis  system  and  a  target-body  axis 
system.  Each  system  is  right  handed  with  z-axes  that  are  positive  upward  (relative  to  Earth). 

The  Earth-axis  system  is  fixed  relative  to  Earth’s  surface  (assumed  to  be  flat)  with  its  origin 
located  at  some  fixed  but  arbitrary  point  relative  to  Earth’s  surface.  The  x-y  axes- plane  is 
located  parallel  to  the  plane  of  Earth’s  surface  and  the  z-axis  is  oriented  positive  upward.  Vec¬ 
tors  represented  in  Earth  coordinates  have  no  superscript. 

The  target-body  axis-system  is  fixed  relative  to  the  target  with  its  origin  at  the  target’s  eg.  The 
x-axis  is  parallel  to  a  longitudinal  lihe-of-symrhetry  for  the  target,  positive  forward.  The  y-axis 
is  parallel  to  a  transverse  line-of-symmetry  for  the  target.  The  z-axis  is  positive  upward.  Vec¬ 
tors  represented  in  target-body  coordinates  are  superscripted  with  B. 

Figure  5.1  illustrates  these  axis  systems. 

5.2.2  Motion  Equations  for  the  Target’s  Center  of  Gravity 

Following  the  notation  in  Ref.  80,  the  six  equations  of  eg  motion  have  the  form 
(fix i  dxi 

M,  +  N,  +  R,x,  -  F,(w,  t ),  i  -  1,  ...  ,  6. 
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In  Eq.  (5.1)  we  use  the  notation  /"/(to,  t)  to  describe  each  input  force  and  moment, 
the  F,(ot,  t )  are  stochastic  processes.  For  convenience  we  rewrite  Eq.  (5.1)  as 


+  2ht 


dXi 

Hi  + 


fiioj,  t),  i 


1, 


We  assume  that 


(5.2) 


where  2 h,  represents  a  damping  factor  and  i>(  represents  a  natural  undamped  frequency  of  oscillation  of 
the  target.  Writing  Eq.  (5.2)  in  state-variable  format  we  have 


where 


and 


%  -  A, It,  +  Bj, 


X,  “  (x,,  X,)T 


f,  -  (0,  f,)T. 


(5.3a) 


(5.3b) 


(5.3c) 


(5.3d) 


■•2 


We  define  the  state  of  the  target-motion  system  as  the  vector  of  position  and  velocity  variables  of 
the  six  motion  equations: 

x  =  Or,  x,  y,  y ,  z,  z,  <*>,  0,  0,  0,  <i»)r  (5.4) 

where  x,  y ,  and  z  are  the  linear  coordinates,  and  d>,  0,  and  are  the  angular  coordinates  of  the  target's 
eg.  Because  the  motions  are  assumed  to  be  uncoupled,  the  state  vector  becomes 

x  =  (x,  u,  y,  i>,  z,  h\  </>,  p ,  0,  q ,  </,,  r)T  (5.5) 

where  «,  i>,  w,  /?,  (7,  and  r  are  the  usual  variables  of  a  six  degrees-of-freedom  model;  that  is,  m,  v , 
and  w  are  the  linear  rates  of  change  and  p,  <7,  and  r  are  the  angular  rates  of  change.  The  system  dif¬ 
ferential  equations  can  now  be  written  in  the  form 

x  *■  Ax  +  Bf  (5.6) 
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where  A  is  a  12  x  12  matrix  with  the  form 


A  - 


Az 


'•  A* 


with  the  Aj  defined  by  Eq.  (5.3b).  The  B  matrix  is  the  12  x  6  matrix  (from  Eq.  (5.3c)) 


(5.7) 
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1 
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(5.8) 


The  input  vector  is 


/  ”  (/**  /jm  fv>  /*»  /*>  /*)r- 


(5.9) 


For  given  h,  and  v,  coefficients  Eq.  (5.6)  can  now  be  solved.  The  solution  is 

xU)  —  <Mf)x(0)  +  J*Q  4>(f  -  t)B  f(r)dr 
where  is  the  12x12  block  diagonal  matrix 


<Mf)  - 


<M') 


<M*) 


V') 


4>g(f) 


0*(f) 


(5.10) 


(5.11) 


Each  4>,(f)  in  Eq.  (5.11)  is  the  2  x  2  transition  matrix  associated  with  the  appropriate  A,  from  Eq. 
(5.3a). 

5.2.3  Motion  Equations  for  the  Unit-Scatterer  Locations 

The  location  of  the  *th  US  (rth  point)  on  the  target  (in  Earth  coordinates)  is  given  by 
R,(t)  -  Rt( 0)  +  f  '  Vt{r)dr 


-  (RXi,  R2)t 


(5.12) 
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where  V,(t)  is  the  velocity  of  the  /th  point,  and  /?,  (0)  is  the  initial  location  of  the  /th  point.  Now 

Vj(t)  =  Vcgit)  +  wit)  x  R*  (5.13) 

where  P<*(f)  is  the  linear  velocity  of  the  target’s  eg  iVcgit)  =  iu,  v,  w)r),  wit)  is  the  angular  ve¬ 
locity  of  the  target  (S(r)  =  ip,  q,  r)T),  and  R?  =  (a,  r,)r  is  the  location  of  the  /th  point  relative 
to  the  target  eg  in  target-body  coordinates  (which,  for  rigid  bodies,  is  constant)  so  that 

R,it)  =  R,i 0)  +  Jo'  [pt(r)  +  Sir)  x  j?*(T)]rfr.  (5.14) 

Expanding  the  cross-product  term  we  have 

w  x  R*=  iztq  -  y,r)i  +  ix.r  -  z,p)j  +  iy,p  -  X/q)k.  (5.15) 


Equation  (5.14)  can  now  be  written  in  the  form 


Rx  (0) 

i 

u  +  iz,q  -  y,r) 

Ryit) 

= 

Ry.iO) 

+  /0 

v  +  ix,r  -  z,p) 

dr 

R2  it) 

RAO) 

w  +  iy.P  ~  x,q ) 

Rx  (0) 

r  t 

1  0  0  0  2, 

-y, 

= 

Ryi  0) 

+  /o 

0  1  0  -2,  0 

X, 

R‘i0) 

0  0  1  y,  —Xj 

0 

dr . 


By  the  decoupling  assumption  (Eqs.  (5.4)  and  (5.5)),  Eq.  (5.16)  becomes 


U 

V 

W 

P 

Q 

r 

/?,(/)  -  R,i0)  +  /' 


1  0  0  0  z,  —y, 

0  1  0  -r,  0  x ) 

0  0  1  >>,  — Xj  0 


dr 


(5.16) 


(5.17a) 


-  $«»  +  f0'  Tjdr 
=  Rti0)  +  7;  [7(f)  -  7(0)] 


(5.17b) 

(5.17c) 


where  7  is  defined  by  Eq.  (5.4)  and 
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(5.18) 

0  0 

0 

0 

1 

0 

y, 

0  —Xj 

0  0 

0 

Substituting  Eq.  (5.10)  into  Eq.  (5.17c)  we  have  the  position  of  the  /th  US  as 
R,it)  -  £,(0)  +  T,  |<l>(f)7(0)  +  /o'<t>(f  -  r)Bfir)dr\  -  7)7(0) 

-  R,i0)  -I-  T,  |(<t>(f)  -  <t>(0))  +/0'<b(/  -  r)Bfir)dr\. 


(5.19) 
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Next  we  wish  to  obtain  the  components  of  the  impulse  response  matrix  for  the  motion  of  the  ith 
scatterer  (in  Earth  coordinates).  Relabeling  (see  Eq.  (5.11))  the  transition  matrix  as 


where 


we  have 


1>(r) 


*/,21  <**(,22  ' 


(5.20) 


(5.21) 
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(5.22) 


so  that  the  impulse  response  matrix  becomes 
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5.3  Motion  Equations  Given  the  Sea  Spectrum 

Although  ship-motion  studies  date  back  to  Froude  [81]  and  Krylov  [82],  it  has  only  been  within 
the  last  30  years  that  ship  motions  in  realistic  seas  have  become  predictable  with  some  confidence  [83], 
The  sea  surface  (waves)  causes  most  of  the  nonmaneuvering  ship  motion  and  the  major  difficulty  in 
modeling  ship  motion  was  the  lack  of  an  accurate  model  of  this  surface.  The  early  sea-surface  models 
were  deterministic  and  therefore,  because  of  the  complexity  of  the  surface,  quite  limited.  In  1952,  a 
stochastic  model  for  the  sea  surface  was  introduced  by  Pierson  [84],  Then  in  1953,  St.  Denis  and  Pier¬ 
son  [80]  introduced  a  ship  motion  model  that  coupled  the  stochastic  sea-surface  model  with  linear 
ship-motion  equations  to  produce  a  method  of  describing  ship  motion  that  is  still  used  today.  The 
linearize-first  model  of  this  section  is  based  on  this  method. 

5.3.1  Assumptions 

There  are  four  fundamental  assumptions  made  by  St.  Denis  and  Pierson  that  allow  the  lineariza¬ 
tion  of  the  ship  motion  equations.  They  are  [80]: 

(1)  The  theory  is  restricted  to  vessels  of  the  displacement  type  (e.g.,  this  excludes  planing  and 
hydrofoil  craft). 

(2)  The  Froude-Krylov  hypothesis  holds:  the  waves  act  on  the  ship  but  the  ship  does  not  act  on 
the  waves. 

(3)  The  theory  is  restricted  to  uncoupled  motions,  and  the  motions  are  assumed  to  be  independent 
(e  g.,  this  restricts  the  theory  to  vessels  where  the  water-planes  are  quasi-symmetric  fore  and 
aft). 

(4)  The  motion  response  is  assumed  to  be  a  linear  function  of  the  exciting,  restoring,  damping,  and 
inertia  forces  (e.g.,  this  restricts  the  theory  to  vessels  that  are  wall  sided). 

Given  these  assumptions,  each  of  the  six  (decoupled)  motions  can  be  described  by  a  linear  equation  of 
the  form 

+  2h  ^  +  tt/(/).  (5.25) 

dr  at  M 

Of  the  six  possible  motions  (referred  to  as  roll,  pitch,  yaw,  surge,  sway,  and  heave),  three  are 
dominant  contributors  to  the  ship  motion:  heave,  pitch,  and  roll.  However,  for  generality  we  include 
all  six  motions  in  the  solutions. 

5.3.2  Sea  Surface  Model 

As  stated  previously,  the  Pierson  sea-surface  model  is  the  basis  for  describing  ship  motion.  This 
model  describes  the  height  of  the  surface  (wave  height)  relative  to  its  mean-level  as  a  function  of  posi¬ 
tion  and  time.  The  model  is  a  second-order,  Gaussian  stochastic  process  with  a  spectrum  that  is  direc¬ 
tional.  The  spectrum  is  the  most  important  parameter  of  this  model,  and  various  models  for  it  have 
been  developed  including  those  by  Neumann  [85],  Bretschneider  [86],  and  Pierson  and  Moskowitz 
[87],  In  general,  experimentally  measured  wave  height  spectra  have  the  following  major  characteristics 
(to  >  0): 

•  The  spectra  are  smooth. 

•  There  is  a  dominant  frequency  (single  mode). 

•  The  spectral  shape  is  "bandpass." 

•  There  is  asymmetry  to  the  spectra;  there  is  more  energy  in  the  higher  frequencies. 
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Figure  5.2  is  an  illustration  of  a  spectrum  exemplifying  these  characteristics.  Next  we  present  a  descrip¬ 
tion  of  the  sea  surface  using  the  nonstandard,  oceanographic  approach. 


TWO-PARAMETER  BRETSCHNEiDER  WAVE  SPECTRA 
SIGNIFICANT  WAVE  HEIGHT  -  10  ft 


WAVE  FREQUENCY,  u>  { RADIANS/ S) 

Fig.  5.2  —  Sample  Bretschneider  wave  spectrum:  T0  is  the  modal  wave 
period  (from  Ref.  78) 

Given  a  process  t j(t)  that  is  stationary,  an  autocovariance  R  (r)  and  a  spectrum  F(w)  for  the  pro¬ 
cess  can  be  defined  by  the  relations 


R(t)  - J”  eJaTF(u)dw 
F(w)  -  v-  f~  e~JmrR  (r)dr. 

7ir  oo 


(5.26a) 


(5.26b) 


Note  that  associating  the  27t  with  the  forward  transform  allows  the  variance  to  be  defined  as  simply  "the 
integral  of  the  spectrum."  The  "spectrum"  of  the  sea  surface,  S( a>),  is  defined  for  positive  frequencies 
only  by 

S(w)  -  2 FM.  (5.27) 

Because  wave  height  is  a  real  process,  Eqs.  (5  26)  are  usually  expressed  as  the  real  cosine  transforms 

R  (t)  —  J  cos  b>TS{(o)da>  (5.28a) 

Jo 


1  c 00 

SUti)  ™  —  I  cos  cur/?  (t)</t. 

IT  Jo 


(5.28b) 


The  Gaussian  sea-surface  is  represented  by  the  stochastic  integral 

r)(t)  —  J*Q  cos  (a >t  +  ^(a>))  V2S(ai)</o>.  (5.29a) 

The  integral  in  Eq.  (5.29)  is  defined  as  a  quadratic-mean  limit  of  a  sequence  of  random  partial  sums. 


■»}(/)—  l.i.m.  T  COs[<U2/+if  +  ^(<U2/+i)]  >/ S’ (co2/+|)  (o>2;-*-2  “  «2/)> 
“2  AT*"*  “ 

*2V+2"“2W“° 


(5.29b) 
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where  <£(<o,)  is  a  sequence  of  independent  random  variables  that  are  uniformly  distributed  over 
(0,  2n).  This  integral  was  formally  introduced  by  Levy  [88]  and  was  used  in  sea-surface  modeling  by 
Pierson  [89]  (and  noise-current  modeling  by  Rice  [90]).  It  allows  the  intuitive  interpretation  that  17(f) 
is  made  up  of  an  infinite  number  of  randomly  phased  sinusoids  with  amplitudes  that  are  determined  by 
the  sea  spectrum.  Most  importantly,  Eq.  (5.28a)  is  useful  for  simulating  samples  of  the  process  77(f)  by 
making  the  approximation 

N  - - 

7) (f )  ==  Tja(f)  —  £cos  [a>,f  +  4>(to,)]  ■\/2S(to/)Acu,.  (5.30) 

1-1 

Also,  samples  of  the  solution  of 

$ +  u  f  +  « -  ’(') 

can  be  simulated  by 

x(t)  =  £x,(f) 

1-1 

where  x(t )  is  the  solution  of  Eq.  (5.31)  for 

77(f)  =  77fl(f)  =  cos  [<u,f  +  0 Co> ( ) ]  -\/2S’(<«>()Aci>,.  (5.33) 

This  is  the  equation  that  we  use  in  Section  6  to  simulate  the  various  ship-motion  processes. 

Although  Eq.  (5.29a)  describes  the  vertical  motion  (force  magnitude)  of  the  sea  at  some  location 
on  the  surface,  it  contains  no  information  describing  the  "direction  of  propagation"  (force  direction)  of 
the  sinusoidal  components.  Though  we  do  not  require  a  two-dimensional  representation  in  later 
chapters,  for  completeness  we  present  a  brief  description  of  the  one  commonly  used.  This  model  also 
follows  the  nonstandard  formulation. 

The  directional  effects  are  accounted  for  by  a  simple  heuristic  modification  of  the  sea-surface 
spectrum.  (Note  that  because  of  the  nonstandard  formulation,  the  directional  sea-spectrum  is  not  a 
spectral  density  in  the  sense  of  two-dimensional  random  fields  [91]  and  that  it  must  also  be  used  with 
caution.)  The  directional  sea-surface  spectrum  S(a>,  9)  where  9  represents  an  angle  measured  relative 
to  a  fixed  direction,  has  been  described  [92],  based  on  empirical  measurements,  by 


(5.31) 

(5.32) 


S(«,  9) 


—  [1  +  (0.50  +  0.82  e~l/2*  w)  cos  29 

7 r 


+(0.32  e->/2*'4“4“4)  cos  40], 


tog  ^  to  ^  °° 

2  2 
otherwise. 


(5.34) 
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(The  gravitational  constant  is  g,  and  the  mean  wind  speed  is  u.)  Another  model,  used  by  Pierson,  Neu¬ 
mann,  and  James  [93],  is 


SU,  0) 


cos2* 

0 


ft»o  ^  w  ^ 

T  ^  0  ^  2 
otherwise 


(5.35) 


where  g  ™  980  cm/s2,  and  C  =  1.53  x  104  cm2s-5.  Not  only  is  Eq.  (5.35)  simpler  than  Eq.  (5.34), 
but  it  has  been  found  to  be  a  very  good  approximation  for  actual  seas  [92,94]  where  high-frequency 
components  (frequencies  greater  than  about  1  Hz)  are  not  significant.  The  "spectral  representation"  of 
the  sea  surface  at  some  point  (x,  y)  as  a  function  of  S(o>,  0)  is  [89] 


i)(x,  y ,  t)  - 


cos 


—  (x  cos  *  +  y  sin  0)  -  tot  +  <t>(<o,  0) 

g 


V25(w,  9)dtod0. 


(5.36) 


The  usefulness  of  Eq.  (5.36)  is  that  it  represents  r/O )  as  a  summation  of  sinusoids  of  frequency  to  and 
direction  0,  traveling  in  the  mean  wind-direction  (0  “  0),  distributed  in  energy  and  direction  as 
described  by  Eq.  (5.35).  Figure  5.3  is  a  qualitative  illustration  of  this  interpretation. 


Fig.  5.3  —  Qualitative  description  of  the  two-dimensional  wave 
spectrum  (from  Ref.  94) 


5.3.3  Response  Amplitude  Operators 

As  we  saw  in  the  previous  subsection,  the  sea  surface  can  be  modeled  as  a  Gaussian  process  so 
that  different  seas  are  described  by  specifying  their  associated  "spectra."  Building  on  this  fundamental 
role  for  the  spectrum,  St.  Denis  and  Pierson  [80]  developed  a  frequency  response  model  of  ship 
motions.  The  system  response  functions  describing  the  interactions  between  the  ship  and  sea  are  called 
response  amplitude  operators  (RAOs).  The  RAOs  are  the  frequency  response  of  the  various  degrees  of 
freedom  of  the  ship  to  the  wave  motion.  Given  the  wave-height  directional  spectrum,  S{a> ,  0),  and  the 
ith  degree-of-freedom  (dof)  RAO,  /?,(&>,  0),  the  associated  response  spectrum  is 

/w/2  p*/2 

_w/J  S,(to,  0)d0  -  J  w/2  |rt,(a>,  0)|2S(o>,  0)d0,  i  -  1 . 6.  (5.37) 
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In  terms  of  system  theory,  one  can  think  of  /?,(&»,  0)  as  the  frequency  response  of  a  system  whose 
input  is  the  stochastic  wave  motion  along  direction  0  and  whose  output  is  the  /'th  dof  ship  motion.  The 
total  motion  due  to  the  sea  is  then  computed  by  incoherent,  linear  superposition  (Eq.  (5.37))  of  the 
contributions  from  each  wave  direction.  Figure  5.4  illustrates  this  concept  assuming  a  finite  number  N 
of  wave  directions. 


S  (uj, 

S  (a),  02* 

S  (cu,  0jg) 


S,  (u>) 


directional  sea,  Eq.  (S.38) 

When  a  ship  moves  in  a  directional  sea,  the  response  of  the  ship  to  the  sea  is  a  function  of  both 
ship  speed  and  heading  relative  to  the  sea-wave  direction.  The  RAOs  are  derived  with  reference  to  a 
system  of  coordinates  fixed  in  space  and  oriented  so  that  the  x-axis  is  positive  in  the  direction  of  the 
ship’s  heading.  To  account  for  the  speed  and  direction  effects,  a  transformation  from  the  fixed  coordi¬ 
nate  system  to  the  ship’s  coordinate  system  is  introduced.  This  is  simply  a  frequency  transformation 
that  maps  w  into  the  "frequency  of  encounter"  a>e  by  [80) 

<ue  —  at  —  —  V  cos  x  (5.38) 

g 


where  V  is  the  ship’s  speed  (F  >  0),  *  is  the  ship’s  heading  (relative  to  the  sea),  and  g  is  the  gravita¬ 
tional  constant.  Equation  (5.37)  now  becomes 

rn/1  /•»/  2 

S'(<°)  “  J_w,2  0)d0  =  J  !/?,(<*>«.,  0)l2S(o>,  0)d0 ,  i  -  1,  ...  ,  6.  (5.39) 


6.  UNIT-SCATTERER-MOTION  SIMULATION 

Our  goal  in  this  chapter  is  to  develop  a  simulation  of  the  6  dof  motions  for  the  US  geometric- 
centers  which  are  caused  by  sea  surface  motion  and  ship  maneuvering.  It  is  well  known  that  sea- 
induced  motions  can  be  accurately  described  as  zero-mean  Gaussian  processes.  The  DTNSRDC  has 
used  the  St.  Denis-Pierson  model  of  ship  motion  (described  in  Section  5)  to  generate,  via  simulations, 
libraries  of  ship-motion  "spectra,"  which  complete  the  specification  of  the  processes.  These  "spectra" 
have  been  generated  for  various  ships  assuming  various  sea  conditions  and  ship  speeds.  The  method  of 
RAOs  described  in  Sec.  5.3.3  is  used  to  generate  the  ship-motion  "spectra"  and  time-domain  realizations 
for  the  6  dof  motions  of  various  locations  on  each  ship.  Some  results  of  this  work  are  presented  in 
Refs.  78  and  79.  Because  we  are  interested  in  studying  the  ship  motions  per  se  and  not  in  studying  the 
physical  processes  by  which  these  motions  are  generated,  we  base  our  motion  model  on  the  DTNSRDC 
"spectra." 


6.1  Motion  Model 


Because  we  are  assuming  rigid  body  motion  and  that  the  sea-induced  motions  are  small  and 
uncoupled  (assumptions  consistent  with  the  St.  Denis-Pierson  model),  we  only  require  knowledge  of 
the  eg  motions  and  of  the  scatterer  locations  relative  to  the  ship’s  eg.  Given  these,  we  use  Eq.  (5.17) 
to  yield  the  position  of  the  /th  scatterer  (relative  to  an  Earth-coordinate  system)  as  a  function  of  time. 
We  modify  those  equations  here  to  be 

/?,(/)  -  *,(0)  +  r,[xf,(t)  -  xf,(0)]  (6.1) 
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and 


T, 


1  0  0  0  z,  -y, 

0  1  0  -z,  0  X/ 

0  0  1  y,  - x i  0 


(6.2) 


where  ?„(/)  is  the  position  vector  (reduced  state  vector)  of  the  ship’s  eg,  and  R,(t)  is  the  position  vec¬ 
tor  of  the  <th  scatterer,  both  in  Earth  coordinates. 


Equation  (6.1)  was  derived  assuming  small  uncoupled  motions.  For  our  studies  we  are  interested 
in  scattering  from  various  aspect  angles.  These  aspect  angles  result  from  apparent  translations  and/or 
rotations  of  the  ship  relative  to  the  radar.  For  example,  the  ship  may  be  located  at  some  position  R0 
and  aspect  angle  0o  relative  to  the  radar  boresight  axis  (Fig.  6.1).  These  large  aspect  variations  are  due 
to  our  choice  of  radar-ship  orientation  and  to  ship  maneuvering.  Our  approach  to  accounting  for  these 
effects  is  to  specify  the  sea-induced  motions  in  a  local  coordinate  system  and  then  transform  that  solu¬ 
tion  to  an  Earth-coordinate  system  located  at  the  radar.  This  allows  us  to  retain  the  linear  uncoupled 
solution  Eq.  (6.1),  rather  than  the  general  nonlinear  solution  (from  Eq.  (5.16)), 


R,(t) 


*,«>+/; 


(6.3) 


We  start  building  the  simulation  by  assuming  knowledge  of  the  sea-induced  ship-cg  motion, 
(r).  This  motion  is  given  in  local  (ship)  Earth-axis  coordinates. _This  local  system  is  defined  by  the 
ship-axis  and  the  Earth-axis  systems  being  coincident  for  xcg(t)  —  0.  We  then  locate  the  ship-axis  sys¬ 
tem  in  space  according  to  ship  maneuvering  and  positioning  by  the  translation 

R,  (r)  -  R,  (0)  +  Jo'  Vs  (r )dr  (6.4a) 

and  the  rotation 


T,(t) 


COS  0  cos# 
sin  4*  cos  9 
’-sin  9 


cos  0  sin  9  sin  0  ~  sin  0  cos  0 
sin  0  sin  9  sin  0  +  cos  0  cos  0 
cos  9  sin  0 


cos  0  sin  9  cos  0  +  sin  0  sin  0 
sin  0  sin  9  cos  0  -  cos  0  sin  0 
cos  9  cos  0 


(6.4b) 


where  F,(r)  is  the  linear  velocity  of  the  ship’s  eg  in  Earth-coordinates,  and  0,  9 ,  and  0  here  denote 
ship-axis  system  orientation  relative  to  the  radar-site  Earth-axis  system,  not  components  of  a  motion 
vector.  The  position  of  the  /th  scatterer  now  becomes 
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(6.5a) 

(6.5b) 


Rf(t)  =  r,(/){j^(0)  +  T,[x^(t)  -  4(0)]}  +  Rf(t) 

=  Ts(t)Rf-(t )  +  R?(t) 

where  the  superscript  L  denotes  local  Earth-coordinates  and  the  superscript  E  denotes  radar-site  Earth- 
coordinates.  We  note  that  the  sea-induced  motion  xcg(t)  is,  in  general,  a  function  of  ship  maneuver¬ 
ing.  Ship  rotations  due  to  maneuvering  (included  in  Ts(t))  are  not  accounted  for  in  the  DTNSRDC 
data,  and  it  is  well  known  that  in  hard  turns,  ships’  responses  to  the  sea  are  not  the  same  as  in 
straight-line  motion.  Therefore  we  restrict  our  studies  to  maneuvers  characterized  by  straight-line 
motions  with  constant  forward  speeds  and/or  turning  maneuvers  in  which  we  assume  that  the  turns  are 
small  enough  that  their  resulting  motions  are  independent  of  the  sea-induced  motions.  This  allows  us 
to  use  the  solutions  Acg(r)  that  are  tabulated  by  DTNSRDC  as  a  function  of  ship  speed.  Once  we 
choose  a  ship  speed  (i.e.,  Vs(t))  we  then  obtain  the  appropriate  xcg(t)  processes. 

Summarizing  to  this  point:  we  begin  our  scatterer-motion  model  by  assuming  that  the  sea- 
induced  and  maneuvering  motions  are  independent  except  that  Pj(f)  determines  the  choice  of 
cataloged  motions  xcg(t).  We  then  solve  Eq.  (6.1)  for  sea-induced  scatterer  motions  and  then  Eq. 
(6.5a)  for  maneuvering  effects.  Note  that  the  sea-induced  motions  can  be  viewed  as  Gaussian  pertur¬ 
bations  of  the  maneuvering  motion. 

6.2  Model  of  the  Spectrum 

Equations  (6.5)  are  our  fundamental  equations  for  calculating  phase  due  to  scatterer  motion. 
They  require  three  parameters  for  solution:  ship  position  and  orientation  due  to  maneuvering,  scatterer 
locations  relative  to  the  ship’s  eg,  and  ship’s  eg  motion.  Maneuvering  is  application-dependent  and  the 
scatterer  locations  are  determined  by  target  structure.  The  third  parameter,  xcg(t),  is  always  a  zero- 
mean  Gaussian  process  with  spectra  that  are  functions  of  ship  type,  ship  speed,  and  sea  conditions. 
This  section  discusses  the  simulation  of  xcg(t). 

We  simulate  xcg(t)  by  using  the  method  discussed  in  Sec.  5.3.  There  it  is  shown  that  we  can 
approximate  a  Gaussian  process  with  arbitrary  energy  ’’spectrum"  S(u>)  by  Eq.  (5.30), 

N  - - 

v(t)  =  £V2SWAo>,  cos  (<Oii  +  <£(&>,))  (6.6) 

/-i 

whe-e  the  </>(«,)  are  independent  random  variables  uniformly  distributed  on  (0,2ir).  The  processes  are 
completely  determined  by  specifying  the  "spectrum"  S(ot)  for  each  component  process.  The  "spectrum" 
can  be  defined  by  an  analytical  expression  or  by  tabulated  data. 

Clearly,  the  motion  spectra  are  our  fundamental  model  inputs.  DTNSRDC  has  calculated  and 
tabulated,  on  microfiche,  ’’spectra"  for  various  ships  and  operating  conditions.  An  example  of  this  data 
is  shown  in  Table  6.1.  Reference  78  describes  the  DTNSRDC  method  and  illustrates  it  with  specific 
examples:  eg  motion  for  the  DD963  moving  at  10  knots  in  two  sea  conditions.  Because  detailed  data 
for  these  examples  are  presented  in  Ref.  78,  we  use  one  of  them  to  obtain  the  motion  processes  for  the 
studies  of  this  report.  In  particular,  we  use  the  processes  resulting  from  the  DD963  moving  at  10  knots 
in  a  fully  developed,  shortcrested  (irregular)  Bretschneider-spectrum  (Table  6.2)  sea  and  moving  at  a 
30°-heading  relative  to  the  dominant  sea-wave  direction.  The  associated  motion  "spectra"  are  shown  in 
the  third  column  of  Fig.  6.2.  Figure  6.3  illustrates  the  cosine-squared  spectral  weighting  used  to  pro¬ 
duce  the  shortcrested  wave  spectrum.  Clearly,  the  spectral  shapes  vary  considerably  as  a  function  of 
ship  type,  speed,  heading,  and  sea  conditions.  Although  the  simulation  can  use  the  DTNSRDC  data 
directly  in  specifying  S( <o)  in  Eq.  (6.6),  we  prefer  to  use  an  analytical  expression.  This  allows  analytical 
formulations  in  analyses  involving  the  motion  processes  and  more  flexibility  in  the  simulation.  How¬ 
ever,  in  doing  so,  we  must  approximate  the  spectral  forms.  This  is  not  a  significant  problem  because 
we  are  concerned  with  very  short  time  periods  (on  the  order  of  1  s)  for  the  scattered  signal  to  evolve 
relative  to  the  time  required  for  the  motion  processes  to  evolve  significantly  (on  the  order  of  tens  of 
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Table  6.1  —  An  example  of  the  Tabulated  DTNSRDC  Motion-Spectra  (from  Ref.  78) 


DDH3 

SHORTCRESTED 

RMS  VER  DISP  IN  FEET/ENCOUNTERED  MODAL  PERIOD,  T 0£,  IN  SECONDS 
CENTER  OF  GRAVITY  -  251.7  FT  FORWARD  OF  AP  AND  21.9  FT  FROM  DL 
SHIP  HEADING  ANGLE  IN  DEGREES 

VT 


0 

0 

15 

30 

45 

60 

75 

90 

105 

120 

135 

ISO 

165 

100 

7 

.034/9.5 

.040/0.3 

.074/  7.9 

.101/  7.5 

.124/  7.0 

.141/  7.0 

.149/  7.0 

.140/  7.0 

.130/  6.1 

.120/6.7 

.090/  6.7 

.075/  6.7 

.064/  6.0 

9 

.071/12.6 

.001/12.1 

.105/11.6 

.132/10.1 

.155/9.0 

.172/  9.5 

.101/  9.2 

.179/  9.0 

.169/  9.0 

.151/  9.0 

.120/  9.0 

.107/  9.0 

.090/  9.2 

11 

.114/14.0 

.121/14.0 

.130/13.4 

.160/13.1 

.100/12.6 

.195/12.1 

.202/11.6 

.200/11.2 

.190/11.2 

.176/11.2 

.155/11.2 

.139/11.2 

.133/11.2 

13 

.150/15.7 

.155/15.3 

.167/15.3 

.104/14.6 

.200/14.3 

.211/13.7 

.217/13.6 

.215/13.1 

.207/12.0 

.194/12.6 

.100/12.6 

.160/12.6 

.163/12.6 

IS 

.1.77/17.5 

.101/17.0 

.190/17.0 

.202/16.5 

.214/61.1 

.223.15.7 

.220/15.3 

.226/15.0 

.220/14.6 

.210/14.3 

.199/14.3 

.190/14.3 

.107/14.3 

17 

.195/19.0 

.190/19.0 

.205/19.0 

.214/10.5 

.224/15.0 

.231/17.5 

.236/17.0 

.233/17.0 

.224/16.5 

.206/16.5 

.220/16.5 

.206/16.5 

.203/16.5 

V9 

.200/19.6 

.210/19.6 

.216/19.6 

.223/19.6 

.230/19.6 

.236/19.0 

.230/19.0 

.230/11.5 

.234/10.0 

.220/10.0 

.222/10.0 

.217/17.5 

.215/17.5 

21 

.219/22.4 

.221/22.4 

.225/22.4 

.230/21.7 

.236/21.7 

.240/20.9 

.242/20.9 

.242/20.3 

.239/20.3 

.234/12.6 

.230/19.6 

.216/19.6 

.225/19.6 

7 

.032/13.4 

.046/10.1 

.073/  7.9 

.101/  7.5 

.126/  7.0 

.143/  7.0 

.153/  7.0 

.153/  6.0 

.144/6.0 

.127/  6.0 

.104/  7  .0 

.001/  7.1 

.071/  7.3 

9 

.069/15.0 

079/IS.0 

.103/12.1 

.131/10.3 

.157/  9.0 

.176/  9.2 

.107/  9.0 

.100/  9.0 

.100/  0.7 

.164/  0.5 

.143/  0.5 

.124/  0.5 

.116/0.5 

11 

.111/16.1 

.110/16.1 

.136/13.7 

.159/15.3 

.100/13.4 

.197/12.1 

.207/11.6 

.200/11.2 

.201/11.2 

.107/10.0 

.171/10.0 

.157/10.4 

.151/10.0 

13 

.146/17.5 

.151/17.5 

.164/17.0 

.102/16.5 

.199/15.7 

.212/14.6 

.220/13.7 

.221/13.4 

.215/12.0 

.204/12.6 

.192/12.6 

.101/12.6 

.177/12.6 

15 

.173/19.0 

.177/19.0 

.107/10.5 

.200/10.0 

.213/17.5 

.223/16.5 

.229/15.7 

.230/15.0 

.226/14.4 

.217/14.6 

.200/14.3 

.200/14.3 

.197/14.5 

17 

.192/20.9 

.194/20  3 

.202/20.3 

.212/19.6 

.222/19.0 

.230/10.0 

.235/17.5 

.235/17.0 

.232/16.5 

.226/  1.5 

.219/16.5 

.213/  1.6 

.211/16.5 

19 

.205/21.7 

207/21.7 

.213/21.7 

.221/21.7 

.229/20.9 

.235.20.3 

.239/19.6 

.239/19.0 

.237/10.5 

.232/10.0 

.227/10.0 

.222/10.0 

.221/17.5 

21 

216/24.2 

.210/24.2 

.222/24.2 

.220/23.3 

.234/22.4 

.239/21.7 

.242/20.9 

.243/20.3 

.241/20.3 

.237/19.6 

.233/19.6 

.230/19.6 

.229/19.6 

7 

.031/14.6 

045/  0.7 

.072/  0.3 

101/  7.0 

.126/  7.0 

.145/  7.0 

.155/  7.0 

.156/  7.0 

.147/  7.0 

.131/  7.1 

.100/  7.3 

.006/  7.5 

.075/  7.5 

9 

.067/19.6 

077/196 

.101/12.5 

.131/12.1 

.159/  9.2 

.100/  9.0 

.196/  0.7 

.190/  0.5 

.193/  0.3 

.179/  0.3 

.161/  0.3 

.144/  0.3 

.137/  0.3 

11 

.100/203 

.115/203 

133/20.3 

.157/20.3 

.101/12.0 

.200/12.6 

.213/11.6 

.210/11.2 

.214/10.0 

.205/10.5 

.192/10.1 

.100/10.1 

.175/10.1 

13 

.143/20.3 

140/203 

161/20.3 

.100/20.3 

.190/15.7 

.214/15.3 

.224/14.0 

.220/13.4 

.226/12.0 

.219/12.6 

.209/12.6 

.201/12.6 

.197/12.6 

IS 

.170/21.7 

.173/21.7 

.103/21.2 

.197/20.9 

.212/20.3 

.224/17.5 

.235/15.7 

.236/15.3 

.236/15.0 

.220/14.6 

.221/14.3 

.215/14.3 

.213/14.3 

17 

.100/22  4 

.191/224 

191/224 

209/224 

.221/209 

.230/20.3 

.237/14.0 

.240/17.5 

.230/17.0 

.234/16.3 

.229/16.5 

224/16.5 

.222/16.5 

19 

202/24  2 

.200/24.2 

.210/24.2 

.210/23.3 

.227/22.4 

.235/20.9 

.240/20.3 

.242/19.6 

.241/19.0 

.230/10.0 

.234/10.0 

.231/17.5 

.229/17.5 

21 

213/262 

.215/26.2 

219/25.1 

.226/25.1 

.233/24.2 

.239/23.3 

.243/20.9 

.245/20.* 

.244/20.3 

.242/19.6 

.239/19.6 

.236/19.6 

.235/19.6 

7 

030/137 

.044/13.7 

071/137 

100/  7.0 

.126/  7.0 

.145/  7.0 

.153/  70 

.156/  7.0 

140/  7  .1 

.132/  7.1 

.109/  7.7 

.017/  7.7 

.076/  7.7 

9 

064/190 

.075/19.0 

100/190 

.130/190 

.160/  9.5 

.104/  9.0 

.200/  0.5 

.206/  0.3 

.204/  0.3 

.192/  0.5 

.177/  0.5 

.162/  §  5 

.155/  8.5 

11 

104/23  3 

111/23.3 

131/233 

.156/23.3 

.102/19.0 

.204/11.6 

.220/11.2 

.220/10.0 

.229/10.1 

.223/  9.0 

.213/  9.5 

.204/  9.2 

.201/  9.2 

IS 

130/26  2 

144/26.2 

150/26  2 

150/26.2 

.170/23.3 

.190/23.3 

.217/15.0 

.230/14.0 

.230/13.7 

.239/12.6 

.235/12.4 

.229/12.1 

.223/11.6 

IS 

166/26  2 

.169/26.2 

.100/26.2 

.195/262 

.211/23.3 

.226/190 

.237/15.7 

.243/15.3 

.244/15.0 

.242/14.6 

.237/16.3 

.233/14.3 

.232/14.3 

17 

104/27  3 

107/27.3 

.195/27.3 

.207/26.2 

.220/26.2 

.231/19.6 

.240/19.0 

.245/17.5 

.246/170 

.244/16.5 

.241/16.5 

.238/16.5 

237/16.5 

19 

199/27  3 

201/27  3 

207/27.3 

.216/27.3 

.226/26.2 

.235/23.3 

.242/20.3 

.246/19.6 

.247/19.6 

.246/10.0 

.244/10.0 

.242/17.5 

.241/17.5 

21 

211/20  6 

212/273 

.217/27.3 

.224/27.3 

.232/27.3 

.239/26.2 

.245/23.3 

.240/20.9 

.249/20.3 

.240/19.6 

.247/19.4 

.245/19.6 

.244/19.6 

7 

030/165 

044/165 

071/16  5 

.100/16.5 

.123/16.5 

.144/  7.0 

.154/  7.0 

.155/  7.1 

.146/  7.5 

.130/  7.5 

.107/  7.7 

.084/  7.9 

.071/  S  I 

9 

062/165 

073/16.5 

099/16  5 

.130/16.5 

.161/16.5 

.106/16.5 

.204/16.5 

.212/  0.5 

.211/0.5 

.201/  0.7 

.106/  0.7 

.172/  8.7 

.166/  0.7 

II 

101/33  I 

109/331 

129/165 

.155/16.5 

.103/16.5 

.200/16.5 

.226/11.2 

.230/10.5 

.241/  9.0 

.230/  9.2 

.231/  9.2 

.224/  9.2 

.221/  9.2 

13 

136/33  1 

141/33.1 

.156/33.1 

.176/24.2 

.199/24.2 

.219/16.5 

.236/13.4 

.247/11.1 

.252/12.6 

.252/12.1 

.240/11.2 

.245/11  2 

.243/11.2 

IS 

163/33  1 

167/33  1 

.170/33.1 

.194/33.1 

.211/24.2 

.220/16.5 

.241/16.5 

.250/16.5 

.255/14.6 

.256/14.1 

.256/14.3 

.252/14.3 

.251/143 

17 

182/33.1 

105/31 1 

194/11  1 

.206/33.1 

.219/242 

.233/24.2 

.241/10.0 

.251/17.5 

.255/17  0 

.256/16.5 

.255/165 

.253/ 16 5 

.253/16.5 

19 

197/33  1 

199/31 1 

.206/31.1 

.213/33.1 

.226/33.1 

.236/24  2 

.245/20.1 

.251/19.6 

.256/19.6 

.255/10.0 

.255/100 

.254/ 17.5 

.253/17.5 

21 

209/33  1 

211/13  1 

.216/33.1 

.223/33.1 

.232/33.1 

240/24  2 

.247/24.2 

.252/20.9 

,255/20.3 

.256/19.6 

.255/19.6 

.255/19.6 

.254/196 

NOTE:  V  i«  ship  speed  in  knou  and  7g  is  modal  wave  period  in  aeconda. 
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Table  6.2  —  Bretschneider  Spectrum  and  Associated  Statistics  (from  Ref.  78) 


Single  Amplitude  Statistics 

<r 

Bretschneider  Spectrum  S£(&>) 

Root  mean  square  amplitude,  rms 

1.00 

S{(  w) 

- 

Aa>  5  exp  [— /3/a*4]  in  ft2/s 

Average  amplitude 

1.25 

A 

- 

483.5  CUyTS,  ftV4 

Average  of  highest  1/3  amplitudes,  significant 

2.00 

B 

- 

1 944.5/  T$,  s'4 

Highest  expected  amplitude  in  10  successive 
amplitudes 

2.15 

- 

Average  of  highest  1/3  wave 
heights 

Average  of  highest  1/10  amplitudes 

Highest  expected  amplitude  in  30  successive 

2.55 

To 

- 

Modal  period  of  spectrum,  i.e. 
period  corresponding  to  peak 

2.61 

of  spectrum 

amplitudes 

Highest  expected  amplitude  in  50  successive 
amplitudes 

2780 

Highest  expected  amplitude  in  100  successive  • 

amplitudes 

3.03 

Highest  expected  amplitude  irt  20(7 successive' 

•l  * 

amplitudes 

3.25 

Highest  expected  amplitude  in  1000  successive 
amplitudes 

3.72 

. 

N  *»  Number  of  successive  amplitudes 

CONSTANT  —  yfl  (In  A0I/2,  where  CONSTANT  relates  a  te  the  highest  expected  amplitude  in  N 
successive  amplitudes. 

>  >.  ■' 

The  highest  expected  amplitude  in  N  amplitudes  is  the  nost  probable  extreme  value  in  N  amplitudes.  This  value  may  be 
exceeded  63%  of  the  time.  *  '  ' 

To  obtain  wave  height  or  double  amplitude  statistics  from  rms  values,  multiply  single  amplitude  constants  by  2.0. 
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180“ 

(0.106) 


NOTE:  NUMBERS  IN  PARENTHESES  ARE 
SHORTCRESTING  WEIGHTS. 

Fig.  6.3  —  Shortcresting  scheme  for  the  example  wave  spectrum  (from  Ref.  78) 


seconds).  In  other  words,  a  ship  does  not  move  very  far  in  1  s.  Therefore,  in  generating  the  scattering 
processes  we  generate  very  short  "pieces"  (realizations)  of  the  motion  processes.  We  quantify  this 
observation  in  later  paragraphs. 

6.3  Spectrum  Algorithm  in  the  Simulation 

Visual  examination  of  the  spectra  yields  the  following  observations:  they  typically  are  unimodal, 
are  band-pass,  are  asymetrical  about  the  mode  with  more  energy  in  the  higher  frequencies,  and  have 
negligible  energy  in  very  low  frequencies.  These  characteristics  and  the  desire  for  simplicity  led  us  to 
propose  using  a  spectrum  of  the  form  (following  the  notation  of  Eqs.  (5.28a)  and  (5.28b)) 

F„( o>)-—7 — -r- - rr  +  —7 — -r - (6.7) 

a2  +  (o>  —  g>o)2  a2  +  (at  +  (uq)2 

to  approximate  the  motion  spectra.  For  the  appropriate  choice  of  a»o  and  a,  the  spectral  approximation 
is,  for  w  >  0,  characterized  as  being  unimodal,  band-pass,  symmetric  about  the  mode,  of  simple  ana¬ 
lytic  form,  and  F(0)  «  F(<dq)  (Fig.  6.4).  In  the  simulation  we  alter  the  upper  and  lower  spectral  tails 
to  introduce  assymetry.  (The  author’s  experience  has  been  that  Eq.  (6.7)  often  well  characterizes  the 
spectrum  of  measured  data  and  may  be  a  better  approximation  to  the  spectra  of  physical  processes  than 
the  commonly  used  low-pass  form 
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Fig.  6.4  —  Approximation  spectrum  for  m  >  0  and  S0  (0)  «  S,  (u0) 

especially  when  <u0  35  «  so  that  F(<u)  has  the  form  shown  in  Fig.  6.5;  note  that  Eq.  (6.7)  becomes  Eq. 
(6.8)  when  <u0  —  0.)  The  autocovariance  function  of  the  associated  energy  spectrum 
( Sa(u> )  «■  2F(u>),  <u  >  0)  is 

R„(r)  —  Jo  2 Fa(w)  cos  (btr)</cu 

™  2 ire~a^  cosoigT.  (6.9) 


Fi(.  6.S  —  Approximation  spectrum  for  u  >  0  and  a  =»  «0 


In  fitting  Eq.  (6.7)  to  the  forms  of  Fig.  6.2,  we  specify  five  parameters:  a>0,  a,  o-2  (the  process 
variance),  /,,  and  /„  (lower  and  upper  cutoff  frequencies).  Clearly,  w0  corresponds  to  the  center  fre¬ 
quency  of  the  process.  We  choose  a  based  on  the  concept  of  bandwidth.  If  Sa(0)  «  Sa(«0),  we 
define  the  half-power  bandwidth  BW  as  the  distance  (in  radians)  between  the  frequencies  where 

Sa  («)  —  y  Sa(o> o).  Since 

Sa  («o)  ~  — , 

cc 

we  solve 

s>„.)  -  £  -  o!+(“-Uo)! 
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yielding 


wl/2  -  w0  ! 


b>o  ±  a. 


(6.10) 


The  process  half-power  bandwidth  is  defined  by 


so  that 


(6.11) 


We  modify  Eq.  (6.7)  for  the  simulation  to  introduce  an  asymmetry  effect  by  specifying  asymmetric 
upper  and  lower  cutoff  frequencies.  For  simplicity,  we  choose  the  lower  frequency  to  be  /,  ** 
o>0  -  BW  and  the  upper  frequency  to  be  fu  —  co0  +  2 BW.  The  resulting  spectral  shape  is  illustrated  in 
Fig.  6.6.  Finally,  we  must  specify  the  variance  of  the  process.  Reference  71  gives  the  spectral  forms 
and  the  associated  rms  values  (standard  deviations)  for  the  component  processes.  To  generate  a  reali¬ 
zation  with  the  same  rms  value,  we  first  generate  a  realization  with  a  variance  of  one  and  then  multiply 
by  the  desired  rms  value.  The  variance  a2  of  the  truncated  spectrum  process  is 


!|tan  1  4 


tan  1  2  +  tan  4 +  4  -  tan'1  4-^  -  2|j. 


(6.12) 


For  our  example  (Appendix  A),  u)q/BW  varies  from  about  2.5  to  about  3.5.  For  simplicity,  and  with 
minimal  loss  of  accuracy,  we  set  gj q/ BW  —  3  so  that  <r2  =  4.940.  Therefore  we  normalize  the  simu¬ 
lated  process  by 

I  .  11/2 


f=|t94o)  rm*Mm 

-  0.45  rms  {t>(/)} 

where  rms  {tj (/)}  is  the  desired  rms  value  for  the  simulated  t j(t)  process. 


(6.13) 


rjxjj] 


FREQUENCY  (Hz) 

Fig.  6.6  —  Truncated  spectral  form  used  in  the  simulation 
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7.  RADAR  CROSS  SECTION 


In  this  section  we  develop  equations  for  the  RCS  of  an  N-source  target  in  terms  of  unit-scatterers. 
Multipath  effects  are  accounted  for  by  using  the  model  described  in  Section  4.  The  transmit  and 
receive  polarizations  are  assumed  to  be  the  same;  no  cross-polarization  effects  are  considered.  We  will 
focus  our  attention  on  the  correlation  properties  of  RCS  over  short  time-intervals. 

7.1  Fundamental  Representation 


In  Section  3  we  showed  that  the  RCS  of  a  target  comprised  of  N  USs  is 


■(/) 


“  2^  Pt X/(Ri,  9t)pR 


(7.1) 


where,  for  the  ith  US,  pT  is  the  transmitter  polarization  vector,  pR  is  the  receiver  polarization  vector,  9, 
defines  the  US  orientation  (as  a  function  of  time),  and  R,  defines  the  US-to-radar  range  (as  a  function 
of  time).  For  simplicity  we  assume  that  pT  -  pK  —  (1,  0)r  so  that  (from  Eq.  (3.8a)) 


r(t)~ 


J*\ i, <*/•»/> 


(7.2) 


For  practical  reasons,  measurements  of  ship  RCS  at  low-grazing  angles  must  include  the  effects  of 
sea-surface  multipath.  Therefore  we  include  a  model  of  multijjath  effects  in  our  representation  of  RCS 
by  modifying  Eq.  (7.2)  by  the  complex  multipath  factor  F(R ,  9)  (described  in  Sec.  4.4)  yielding 


N  _  _  _  _  J4u 

tm,«Mi, (*,.*,)*  ' 


k-i 


(7.3) 


As  discussed  in  Section  4,  F2(R,,  ~9,)  was  obtained  assuming  that  the  US  reradiates  equally  toward  the 
radar  and  toward  the  sea-surface  specular  point  (Fig.  4.2).  This  is  a  good  approximation  for  small 
radar-grazing-angles  and  when  the  US  amplitudes  are  not  strongly  dependent  on  elevation.  Otherwise, 
the  elevation  radiation  pattern  must  be  viewed  as  an  antenna  pattern.  This  pattern  is  then  accounted 
for  by  using  the  same  method  as  was  used  to  account  for  the  radar’s  pattern  in  Sec.  4.4. 


7.2  Useful  Equations 


Equation  (7.3)  is  the  fundamental  equation  that  we  use  to  model  RCS.  In  this  section  we  derive 
from  Eq.  (7.3)  several  other^  equations  for  RCS.  These  equations  will  prove  useful  for  analysis  and 
simulation.  Suppressing  the  R ,  and  0,  terms  of  Eq.  (7.3), 


<r(t) 


yu», 


(7.4) 


and  associating  l/^l2  with  AU/, 
cr(r)  - 


!  I* =u„, 


N 

IS 


A,  e 
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Expanding  Eq.  (7.5), 


and  hence 


and 


N 

2 

N 

<r(t)  — 

£  A,  cos  'fri 

+ 

£  A,  sin  V, 

/-i 

i-i 

N  N  N  N 

-  £  A,  cos  £  Aj  cos  'fry  +  £  Aj  sin  T  Aj  sin  'fry 

i-i  j- 1  /-i  Jt i 

N  N 

'll  4,/4y[cos  'I'j  cos  'Py  +  sin  sin  'fry] 

<-i  j- 1 

N  N 

cr(t)  -  £  I  A,Aj  cos  ('fr,  -  'fry) 
i-l  j- 1 


N  N 

cr(t)  =  2  ]£  £  AjAj  cos  (¥,  -  'fry)  + 
<-i ;- 1 
i*j 


(7.6) 


(7.7a) 


(7.7b) 


For  the  /th  US,  let  /3,(  )  account  for  the  orientation  effects  on  phase  and  y,()  the  multipath 
effects  so  that 

*,(£„?,)  “ 

and 

'fr/  -  'fr;  “  ^  (l«/l  -  \Rj\)  +  ~  PjGj)  +  y,{M„  9,)  -  yj{Rj,9j).  (7.9) 

Equations  (7.7a)  and  (7.9)  are  our  fundamental  equations  for  simulating  RCS.  In  the  simulation,  the 
A,(-)  and  fi(-)  functions  are  user  controlled,  the  R ,  terms  are  generated  by  the  ship  motion  (Section  5), 
and  the  >,  (•)  terms  are  generated  by  using  the  model  in  Section  4. 

As  suggested  by  Eqs.  (7.5)  and  (7.2),  we  define  a  complex  scattering  function  s(t)  to  simplify 
notation: 

s(t)  -  %A,  eJ*'.  (7.10) 

/-I 

Note  that  <r(f)  =  s(t)s*(t). 

7.3  Statistics  for  an  Af-Source  Target 

In  this  section  we  investigate  the  first  and  second  moments  of  the  RCS  process.  We  show  that 
over  short  time-intervals  the  RCS  process  can  be  nonstationary  and  analytically  complex  even  when 
simplifying  assumptions  are  made  regarding  the  scatterers.  We  also  investigate  in  what  sense  the  well- 
known  random-phase  model  (resulting  in  Rayleigh-amplitude  statistics)  is  a  limiting  case  of  the  general 
formulation.  We  use  the  notation  of  o-(t)  for  RCS  and  cr,(r)  for  standard  deviation  of  the  ith  phase 
process;  no  confusion  will  result  given  both  the  context  and  the  subscripting  (standard  deviations  only). 

7.3.1  Mean  and  Correlation  Functions 


2  IT 
\ 


2|R,i  +/3,tf()  +  yi(Ri,0i) 


(7.8) 


Because  the  amplitude  functions  are  much  less  sensitive  to  target  motion  than  the  phase  functions 
are,  over  short  time  intervals  they  can  be  well  approximated  by  deterministic  functions  (typically  as 
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constants,  except  for  specular  scatterers).  However,  even  ignoring  the  US  orientation  and  multipath 
effects,  the  phase  functions  are  rapidly  varying,  nonlinear  functions  of  target  motion.  Therefore  we 
view  the  phase  terms  as  stochastic  processes. 


Given  an  iV-source  target  with  deterministic  amplitudes  A  At)  and  stochastic  phases  0,(r),  the  tar¬ 
get  RCS  is 


N  N 

cr(t)  =  £  £  Ai(t)Ak(t)  e 

i-l  *-i 

and  the  associated  mean  and  autocorrelation  functions  are 

£{<r(/)}=  £  £  .4, (/M* (/)£{*• 

i-l  *-i 


JlBjU)  -  «k(t)] 


-akU)] 


}, 


and 


£{«■(/,) cr(f2)}=  £  £  £  £  Ai(ti)Ak(tl)A,(t2)Am(t2)E[eJl(0i~9kH,')  +  (9'~BmHl2)i}. 

i—  1  /c—1  /—I  m— 1 


(7.11) 


(7.12) 


(7.13) 


The  expectations  in  Eqs.  (7.12)  and  (7.13)  suggest  characteristic  functions.  Such  functions  are 
useful  when  the  joint  densities  of  the  phase  processes  are  known.  Writing  the  characteristic  functions 
in  the  form 


<M<»„  tn<»k,  tk)  = 

Qiklm^h  tn  «J*,  4,  w,,  th  (Om,  tm)  -  £{e'"“'#,(',)  +  ■ 


we  have 


(where  we  use  the  notation  (0,  -  0A )  (r)  =  0,(r)  -  0k(t))  and 

=  ?1,  ~  1,  f|,  1,  t2,  -1,  tj). 

Using  this  formulation  in  Eqs.  (7.12)  and  (7.13)  we  obtain 


.v  v 


£{crU))  =  £  A,2(t)  +  £  £  A,(l)Ak(t)<t>ik(  1,  t ,  -1,  t) 

/*1  /-I  *-l 

k 

=  £/l,2(r)  +  2  £  £  A,(t)Ak(t)<t>lk(l,  r,  -1,  /), 

/-l  /- 1  *-/+! 

i*k 

£{<r(/,)a(/2))  =  £  Ai2(tl)A/2(t2)  +  £  £  A^A/H^) 


/-I  /-i 
i*i 


+  XII  I  (Ti)^,(/2)>4^(r2)«l»fAc,m  (1,  r,  -1,  l,  r2,  -i,  ,3), 


I  A  •> 1  /« l  nt—  I 
except  /  =  A- /=m 
and  /-A  *1—  m 


md 


-V  /V 


,v  ,v 


41  r-(t)}  =  £  A*(t)  +  £  £  A,2(t)A,2(t)  +  £  £  rf,2(/M*20) 


/-l  /-I 
i*l 


i- I  A- I 
t*k 


+  X  X  I  X  /l,(/)/l*(r)/l,(r)/tm(/)«I>,*/m(l,  r,  -1,  t,,  1,  r,  -1,  r). 

/*!  A  — ]  /-l  m— 1 
except  /  =A--/—m 


(7.14a) 

(7.14b) 

(7.15a) 

(7.15b) 


(7.16a) 


(7.16b) 


(7.17) 


and  /-A  ^  /•« 
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In  general  the  phase  processes  are  complicated  functions  with  probability  laws  that  at  best  can  only 
be  approximated.  Over  long  periods  of  time  (relative  to  target-motion  time  constants)  they  most  likely 
have  probability  density  functions  that  are  unimodal  and  symmetric  when  they  are  not  viewed  as 
modulo  2 it  (i.e.,  distributed  over  (-<»,  oo)).  This  suggests  a  Gaussian  approximation.  If  the  phases 
can  be  represented  by  joint  Gaussian  processes,  then  the  form  of  the  associated  characteristic  functions 
are  known;  they  are  [95] 

j\a,lThlll)+a,lcT,k(tk)\-~leTlUl)al2+2rlk(ll,tk>allut+<rlUk)'»j‘\ 

lj ,  wk,  tk)  =  e  (7.18a) 

and 


^  iklm  (to  j  ytj ,  (t)  k ,  <0 1 ,  // ,  to  m,  tm)  —  e  p 

(7.18b) 

where 

v,(t,)  -  £{#<(/, •)), 

(7.18c) 

cr}(t ,)  =  £{[0 ,•(//)  ~  i7,(r()J2}, 

(7.1 8d) 

,  ,  £{[0,('/)  '  ((/)][«*(/*)  -  (/*)!) 
r‘k  i,tk  <rM)<Tk{tk) 

(7. 1 8e) 

ft  =  (o>  j  ,  (Ofay  (0  /  ,  W  )  * 

(7.180 

and  h  is  the  covariance  matrix  for  0,,  6k,  0,,  and  9m  with 

Ppq(lp'  lq 

)  =  E{[0p(tp)  -  1? p(tp)][9q(tq)  ~  i7»(^)]}  for  p,  q  -  /,  k,  /,  n. 

(7.1 8g) 

We  note  that  from  Eq.  (7.18a) 

<J>*(1  t  -1  t)  =  e7|’>/<,)-1’*('),~1/2,<r.2(','':r*<,')+<7*(,)1 

(7.19a) 

and  from  Eq.  (7.18b), 

<*W1.  tu  1,  (2.-1,  r2)  =  ejFrU1Fl 

(7.19b) 

where 

F\  =  V/(0  -  Vk(tk)  +  V/(0  -  Vm(U 

(7.19c) 

and 

F2  =  nuU\,  t\)  -  nik(ti,  1 1)  4  1 2)  -  1.  h) 

-  1 1)  +  flkk(t],  1 1)  -  nki(t  1,  t2 )  +  Hkm(t  1,  t2) 

+  hi, (r2,  1 1)  -  Hik(t 2,  t\)  +  h)  ~  h) 

-  //,(/;,  t\)  +/Um*(r2,  /,)  -  t2)  +  n„„(t2,  r2).  (7.1 9d ) 


7.3.2  Random-Phase  Model 

We  define  a  random  phase  model  for  RCS  to  be  an  /V-source  model  in  which  the  phases  are 
represented  by  independent,  uniformly  distributed  (over  (0,  2ir))  processes.  Such  models  are  widely 
used  in  analysing  the  RCS  of  complex  targets.  They  are  useful  because  they  greatly  simplify  the  RCS 
analysis.  They  are  valid  when  the  time  period  of  interest  (for  analysis)  is  large  relative  to  the  correla¬ 
tion  times  of  the  phase  processes. 
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Because  we  are  interested  in  RCS  analysis  and  simulation  over  short  time-intervals,  it  is  important 
to  consider  in  what  sense  the  general  formulation  can  be  replaced  by  a  random-phase  model.  Therefore 
in  the  following  we  examine  the  behavior  of  the  first  and  second  moments  of  RCS  as  the  random-phase 
assumptions  are  approached. 

We  now  consider  the  simplest  random-phase  model:  assume  that  the  scatterer  amplitudes  are 
equal  and  constant  and  that  the  individual  phases  are  white-noise  processes.  To  obtain  the  first  two 
moments  of  this  model  we  write 


s(f)  -  £  A,eJ9‘  -  A,  cos  0,  +  j  jr  A,  sin  0, 


(7.20) 


where  we  recall  that  <r(t)  —  sit)  s*it).  By  the  central  limit  theorem,  the  real  and  imaginary  sums  in 
Eq.  (7.20)  approach  being  Gaussian  distributed  for  large  N.  Further,  they  are  uncorrelated  because 


(7.21a) 


£|  £  A,  cos  9 1  £  Ak  sin  0*J  -  A1  £  £{cos  9 ,  sin  9k ) 

-  A2  J?  £{cos  9 ,  sin  0,}  -  A2  £  -  0. 

Clearly,  the  mean  of  each  sum  is  zero.  Now  consider  the  variances: 
fjj  21^' cos  9 1  IM  £  cos  0,  cos  0*1 


-  A2  £{cos2  0/} 


>S2  9i  _  NA2 


and  similarly. 


£  52i4sin0,  -  A2  £{sin2  0,)  -  — — 


(7.21b) 


(7.21c) 


Because  of  the  white-noise  assumption,  RCS  is  also  a  white-noise  process.  Finally,  the  RCS  process. 


<r(t)  «■  |s(f)|2  —  IS  A  cos  0/j  +  jjj  A  sin  0/J  , 


(7.22) 


approaches  being  exponentially  distributed  because  the  summations  are  independent  Gaussian  processes 
(approximately)  with  zero  means  and  equal  variances  [95].  From  Eqs.  (7.21b)  and  (7.21c), 

£{<x (f))  -  NA2  (7.23a) 

and  because  a- it)  is  exponentially  distributed, 

£{(T2(t)}  -  IN2 A*.  (7.23b) 


Next  we  consider  Eqs.  (7.16)  in  the  context  of  the  random-phase  assumption.  We  assume  that 
the  phases  can  be  represented  by  Gaussian  processes  and,  for  simplicity,  we  assume  that  the  means  of 
the  individual  phase  processes  are  zero.  Therefore  Eqs.  (7.16)  become 

£{<r(f)}  -  T  A,2U)  +  2  £  £  Ai{t)Akit)e~mW}U)~2r*U',HalW  (7.24) 

1- \  *-/+! 
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£Ml,)«r(f2)}  -  £^(/,)^(/2)  + 

i-i  /-I  /-i 

/*/ 

+  ££££  Aj  {t{)Ak  (t\)A)  (12)  Am  (t2)e 

i-i  k-\  /-] 
except  I— k-t-m 
and 


where 


1*1 m  (<i»  h)  ”  +  tfjKfi)  +  <r,2(r2)  +  <r^,(f2) 

-  2rik(tu  ti)<Tj(ti)(rk(ti )  -  2r/m  (t2,  r2)o-,(f2)crm(r2) 
+  2ri,(t\>  r2)o-l(/1)o-,(r2)  +  2/*m(f,,  f^cr^f^o-*,^) 

-  2rim(/1,  /^(r^o-*^)  -  2rki(ti,  r2)o-*(rI)o-/(r2). 
Now  assume  that  the  0,s  are  not  cross  correlated.  In  the  limit  of  large  variances, 

lim  £{<r(r)}  -  £A,Ht) 


(7.25a) 


(7.25b) 


(7.26) 


lim  EMr.VM  -  £  A?Ui)A?{t2)  +  £  £  A^tl)A^t2) 

1-1  t-l  /-I 


+  £  ^AiUi)Ai(t1)Al(tx)A,U1)e~Xnim'',l) 


i-i  /-1 

i*i 


(7.27a) 


where  the  function  in  the  exponent  is 

(*1,  h)  -  W('i)  +  <r}u2)  +  <r}Ui)  +  <r,2(r2)  -  2 /?,(*,,  f2)  -  2/?,(f,,  /2)J  (7.27b) 

(R-tUu  t2)  is  an  autocovariance  function).  (The  third  sum-term  in  Eq.  (7.27a)  remains  because  even 
for  <rf  large,  as  (r2  -  tx)  approaches  zero  Eq.  (7.27b)  approaches  zero;  in  other  words  as  t2  approaches 
?i,  even  for  a2  large,  (0,(ft)  -  0m(/2)l  and  [#*(fi)  ~  0/(t2)]  approach  zero  when  /'  —  m  and  k  —  /.) 
For  ft  —  r2  —  r,  Eq.  (7.27a)  becomes 


lim  £{o-2(r)}  -  £^(/>  +  2  £  £  A,Ht)Ak2(t). 

i- 1  l-l  k-\ 


(7.28) 


Equations  (7.26)  and  (7.27)  imply  that  even  if  the  variances  are  time-varying,  as  they  become  large  the 
RCS  process  becomes  wide-sense  stationary,  assuming  the  amplitudes  are  constant.  Equation  (7.27b) 
shows  that  the  RCS  process  can  remain  correlated  even  for  large  <r,()  if  the  individual  phase  processes 
remain  strongly  correlated. 


Next  we  assume  that  A,{t )  =  Aj(t)  —  A.  Then  Eqs.  (7.26)  and  (7.28)  become 

lim  £{<r(r)}  -  NA2 


(7.29a) 


lim  £{<r2(/)}  -  [N  +  2 N(N  -  1) ]A*  -  2N2 


2  N  ' 


(7.29b) 
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and  Eq.  (7.27a)  becomes 


V  i 


Assuming  equal  and  constant  variance  functions  in  Eq.  (7.30),  ( R ,  —  RA ,  at  —  <rA ,  i 

lim  £{<r(f,)<7(f2)}  -  N2A 4  +  £  e~2l^-^(l^ 

*}—<*>  i-\  /-I 


-  n2a  4 


1  +  N(N-  1)  -2(<ri-^(»,./2)l 

N2 


(7.30) 


1, 


AO, 


(7.31) 


We  see  from  Eq.  (7.31)  that  the  RCS  process  remains  correlated  if  the  individual  phase  processes 
remain  correlated  (not  white  noise). 

For  real  targets  the  phase  variances  remain  finite  and  the  phase  processes  are  not  white-noise 
processes.  It  is  likely  to  be  the  case  that  for  target  motions  that  produce  large  phase  variances  (over 
(-oo,  oo)),  the  phase  process  (due  primarily  to  geometric-center  motion)  will  be  strongly  correlated. 
Clearly,  for  high  sampling  rates  (short  time  between  samples),  the  correlation  of  the  phase  processes 
cannot  be  ignored. 

7. 3. 3  Range-  Variation  Effects 

Next  we  consider  a  specific  model  of  the  geometric-center  motion  to  further  illustrate  the  diffi¬ 
culty  of  modeling  phase  variations  over  short  time-intervals.  Ignoring  US  orientation  and  multipath 
effects  we  have 

o-O)  -  £  £  A,Aj  cos  ~  (|tf>,  r)l  -  |/?,(co,  r)l).  (7.32) 

i-i  j- i  A 

The  range  process,  /?,(&>,  r),  is  (Eq.  (6.5)), 

R*(<o,  t )  -  rs(f)^/-(a>,  t)  +  Rf(t).  (7.33) 

The  stochastic  component  of  Eq.  (7.33)  is  R,Hw,  r);  it  is  a  linear  combination  of  Gaussian  processes 
and  therefore  is  Gaussian.  The  ship  maneuvering  terms,  Ts(t)  and  R£(t),  are  deterministic,  and  over 
small  time-intervals  can  be  approximated  by  linear  functions  of  time.  If  Ts(t)  and  Rf(t)  were  con¬ 
stant,  then  RfU u,  r)  would  be  a  three-dimensional  Gaussian  process.  If  in  addition,_each  component  of 
/?,£(&i,  t)  had  a  mean  of  zero  and  equal  variances  (which  they  do  not),  then  |/?£(a>,  r)|  would  be 
Maxwell  distributed  [95). 

Clearly,  the  |/?,(«,  t)l  process  is  in  general  nonstationary  and  analytically  complex.  Even  assum¬ 
ing  that  the  ship  does  not  maneuver,  all  that  can  be  said  in  general  about  lR,£(w,  t)l  is  that  it  is  a  non¬ 
linear  transformation  of  large-time-constant  Gaussian  processes  (relative  to  a  radar’s  PRO,  is  oscillatory 
and  positive,  and  probably  has  a  density  function  that  is  symmetric  about  a  mode.  Relative  to  a  radar’s 
PRI,  long  time  periods  may  be  required  to  produce  realizations  of  |R,(«o,  r)|  that  exhibit  stationary 
statistics. 

As  a  specific  example,  we  investigate  range-induced  variations  in  the  RCS  phase  for  our  example 
ship.  _We  assume  a  (moderately  high)  radar  frequency  of  10  GHz  (X  -  3  cm).  Over  short  intervals 
the  l£,(w,  f)l  processes  vary  slowly  (a  numerical  estimate  is  made  in  Appendix  B)  and  therefore  we 
will  model  them  as  linear  functions  of  the  form 

I r)|  -  R?(m)  +  K,((0,  t)t.  (7.34) 
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Normalizing  by  4tt/\,  we  obtain  the  phase  angle  process 

<M<o,  t )  —  tli?(co)  +  Cj(o>,  t)t  (7.35) 

where  ifi?(ui)  is  a  uniformly  distributed  random  variable  (over  (0,  27r)).  For  nonmaneuvering  motions, 
c,(w,  t)  is  a  slowly  varying  process  whose  density  function  is  determined  by  the  ship’s  motion.  Because 
this  motion  is  oscillatory,  the  density  function  of  t)  will  likely  be  symmetric  with  mean  equal  to 
zero.  For  short  periods  we  approximate  V,  (to ,  t)  by  a  random  variable  and  the  RCS  process  becomes 


<r(f)  =  £  £  At(t)AkU) 


(7.36) 


The  RCS  which  is  generated  by  the  target  over  the  short  periods  are  "pieces"  of  realizations  of  the 
RCS  process.  Because  the  radar  must  process  only  these  pieces,  we  examine  the  "apparent  variance"  of 
them.  For  an  interval  of  length  T,  the  apparent  mean  of  the  phase  process  is 


i  A'1 

Vt  -  -j.  J_T/2  (<l>?  +  Cit)dt  =  ifi? 


and  the  apparent  variance  is 


aT“  ~f  J_r/2  (<,')2<*=—  • 


(7.37a) 


(7.37b) 


Equation  (7.37b)  confirms  that  which  is  intuitively  obvious:  either  or  both  c,  or  T  must  be  large  for 
the  uniformly  distributed-phase  approximation  to  be  valid.  If  the  standard  deviation  oft ^  is  required  to 
be  greater  than  2nir  for  the  approximation  to  hold,  then  over  a  time  interval  T0  we  require  that 

“TT “  ^  (2nn)2  (7.38a) 


Since  c,  —  4 -it  VJ\, 


.  ^  V48  TT 

C,  >  — = —  n. 
1  o 


1/  V.  -v/3  \ 
vi  &  — =  n 

1  o 


(7.38b) 


(7.39a) 


V3  A 


(7.39b) 


To  illustrate  the  limitations  of  the  uniformly  distributed-phase  model,  consider  a  radar  tracking 
the  stern  of  our  example  ship  at  a  low-grazing  angle.  To  estimate  the  value  of  V,  we  use 


v  _  £{A7*}1/2 


(7.40) 


where  r  is  the  radar’s  PRI.  In  this  scenario,  pitching  motion  has  the  dominant  effect  on  the  variations 
in  the  scatterer  ranges:  from  Appendix  B, 


£{AR2)  =  z2a*£{A02) 

“  2zLx<rj[l  -  e~a'M  coswqtI. 


(7.41) 


54 


NRL  REPORT  8887 


For  t  -  0.002  s,  <o0#  -  0.409  rad/s  and  <t,  -  0.0089  rad  (Appendix  A),  and  -5m  (Appendix 
B),  £{AR2),/2  -  0.00127  m  and  V,  -  0.64  m/s.  If  the  pitch  process,  9,  is  approximately  Gaussian, 
then  n  must  be  no  less  than  about  10.  For  a  3  cm  radar,  Eq.  (7.39a)  says  that  for  the  random-phase 
model  to  apply,  V,  must  be  greater  than  0.51  m  for  T  -  1  s.  We  see  that  in  our  example  this  is  only 
marginally  true  and  that  for  lower  frequency  radars,  it  will  not  be  true. 

8.  GLINT 

The  second  (next  to  RCS)  most  common  projection  of  the  scattered  field  is  glint  (angle  noise, 
angle  scintillation).  This  parameter  represents  the  apparent  angular  location  of  the  target  as  seen  by  the 
radar.  It  is  correlated  with  the  target’s  RCS  and  is  a  function  of  the  radar’s  frequency,  polarization,  and 
antenna  pattern  and  of  the  target’s  shape,  size,  and  composition.  We  make  the  same  assumptions  for 
glint  that  we  made  for  RCS:  the  radar  is  horizontally  polarized  with  fixed  frequency  and  the  target  is  a 
perfect  conductor  whose  polarization  effects  are  accounted  for  in  the  unit-scatterer  amplitudes. 

As  described  in  Section  1,  two  concepts  have  been  used  to  understand  and  analyze  the 
phenomenon  of  glint:  phase-front  distortion  and  Poynting-vector  direction.  In  the  development  of 
glint  models  based  on  these  concepts,  it  was  assumed  that  the  target  was  composed  of  a  finite  number 
of  point-source  or  dipole  scatterers.  Therefore,  these  models  are  directly  applicable  in  our  modeling 
approach.  In  this  chapter  we  use  both  concepts  to  develop  glint  equations  for  a  target  described  in 
terms  of  unit-scatterers. 

Various  authors  have  applied  the  glint  concepts  to  specific  problems,  but  to  our  knowledge  no 
deterministic  analytic-formulation  of  two-dimensional  glint  has  been  made  for  the  /V-source  target  in 
the  presence  of  multipath.  We  develop  this  formulation  for  both  concepts.  The  resulting  equations  will 
be  analytically  complex  but  ideal  for  digital  simulation. 

8.1  Definitions 

Most  existing  models  of  glint  represent  the  target  and  its  motion  in  only  one  plane:  the  target  is 
represented  as  a  line  in  the  azimuth  plane  of  the  radar.  (This  approximation  is  made  because  the  long¬ 
est  dimension  of  most  targets  is  in  that  plane.)  Therefore,  these  glint  models  are  one  dimensional.  We 
refer  to  such  glint  as  azimuth  glint.  The  less-often-used  but  more  general  glint  models  are  two- 
dimensional.  Because  of  our  three-dimensional  representation  of  the  ship  and  of  multipath  effects,  we 
will  develop  a  two-dimensional  glint  model.  The  second  glint  dimension  is  elevation. 

8.1.1  Phase-Front  Gradient 

The  phase-front  definition  was  introduced  by  Howard  (121.  He  showed  that  azimuth  glint  was 
equal  to  the  azimuth  slope  of  the  phase  front  of  the  scattered  field.  Figure  8.1  illustrates  this  concept. 
For  the  one-dimensional  problem,  the  phase  front  is  a  curve  that  passes  through  the  location  point  of 
the  radar.  This  curve  is  defined  by  mapping  the  locations  of  points  where  the  phase  of  the  scattered 
electric  field  is  equal  to  the  value  of  phase  at  the  radar  point.  The  analytical  expression  for  azimuth 
glint  is  derived  by  expressing  the  phase  of  the  scattered  electric  field,  d>,  as  a  function  of  the  angle 
about  the  sight-line  from  the  radar  to  a  fixed  point  on  the  target,  +  (see  Fig.  8.1).  Each  scatterer  is 
assumed  to  yield  identically  polarized  plane  waves  at  the  radar.  The  azimuth  glint  then  becomes,  in 
linear  units. 


1  d<t> 
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RADAR  POINT 


Fig.  8.1  —  Azimuth  glint  visualized  in  terms  of  phase  front  slope 


where  k  is  the  wave  number,  2ir/X.  It  was  shown  by  Howard  [12]  that  for  a  target  composed  of  N 
point  scatterers  on  a  line  (Fig.  8.2),  azimuth  glint,  in  linear  units,  is 

N  V 

£  £  aiajL<  cos  [2*  (I,  -  Lj) cos  0] 

Ga  (*)  -  -  sin  0  •  -  (8.2) 

£  £  fl/fly  COS  [2*U,  -  Lj)  COS  \fl] 
l-i  J-\ 

where  L,  and  0  are  defined  in  Fig.  8.2  and  a,  is  the  amplitude  of  the  ith  scatterer.  The  minus  sign 
arises  because  we  define  i/t  in  the  opposite  sense  than  that  of  Ref.  12. 
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The  obvious  extension  of  the  phase-front-slope  concept  (Eq.  (8.1))  is  the  phase-front  gradient 

/ 

m  8 


rr/»>  -  8  .  /  8  , 

V<p  *  n  —  t  —  —  + 

dr  r  69  r  sin  0  80 


<J>. 


(8.3) 


From  Eq.  (8.3)  and  Fig.  8.2,  we  see  that  in  linear  units  azimuth  glint  is 


if,  i 


8<fr/80 


Vd>  •  n 

Similarly,  the  elevation  glint  is,  in  linear  units, 


sin  8  8d>/8  r' 
-8d>/88 


V<t>  •  h  sin  8  84>/dr 
For  9  —  n/2,  these  equations  reduce  to  the  one-dimensional  problem. 


(8.4a) 


(8.4b) 


8.1.2  Poynting-  Vector  Direction 


The  phase-front  model  indicated  that  the  flow  of  energy  was  not  necessarily  radial  at  the  radar. 
This  led  Dunn  and  Howard  f43]  to  show  that  the  Poynting  vector  from  an  N-source  target  was  normal 
to  the  phase  front.  Therefore  the  Poynting  vector  concept  was  shown  to  be  equivalent  to  the  phase- 
front-slope  concept  in  predicting  glint.  Assuming  that  the  Poynting  vector  of  the  scattered  Held  at  the 
radar  is 


P  -  P,  +  P+  +  P9  -  Prh  +  P+m  +  P9j  (8.5) 

where  Pr ,  P+,  and  P9  are  the  radial,  azimuth-tangential,  and  elevation-tangential  components,  we  see 
from  Fig.  8.3  that  the  azimuth  glint,  in  linear  units,  is 

Ga~\Rs\-^-.  (8.6a) 

*r 

Similarly,  elevation  glint,  in  linear  units,  is 

Ge  -  \R,\  ■  (8.6b) 

‘r 


Fig.  8.3  —  Azimuth  glint  visualized  in  terms  of  the  Poynting  vector 


D.Y.  NORTHAM 


8.1.3  Axis  Systems 

To  obtain  equations  that  are  useful  for  simulating  the  glint  from  an  N-source  target,  we  use  N  + 

2  axis-systems:  an  Earth-axis  system,  a  target-cg-axis  system,  and  N,  US-location-axis  systems.  All  of 
the  axis  systems  are  right-handed  with  the  r-axis  positive  upward  relative  to  Earth’s  surface: 

The  origin  of  the  Earth-axis  system  is  at  the  radar  point.  This  system  is  fixed  relative  to  Earth 
with  the  x  —  y  plane  defined  as  being  parallel  to  Earth’s  surface.  Subscripting  a  vector  with  the 
letter  e  implies  that  the  vector  is  expressed  in  Earth  coordinates.  The  unit  vectors  »,  j,  and  k 
are  the  unit  vectors  of  this  system. 

The  origin  of  the  target-cg-axis  system  is  fixed  at  the  radar  point.  The  x-axis  is  colinear  with 
and  directed  along  the  sight  line  from  the  radar  point  to  the  target’s  eg.  Subscripting  a  vector 
with  the  letter  s  implies  that  the  vector  is  expressed  in  target-cg  coordinates. 

The  origins  of  the  N  US-location-axis  systems  are  fixed  at  the  radar  point.  The  x-axes  are  co¬ 
linear  with  and  directed  along  the  sight  lines  from  the  radar  point  to  the  US  locations.  Sub¬ 
scripting  a  vector  with  an  index  implies  that  the  vector  is  expressed  in  coordinates  of  the  associ¬ 
ated  US-axis  system. 

Coordinates  and  vectors  that  locate  points  on  the  target  relative  to  the  radar  point  are  indicated  by  sub¬ 
scripting  with  the  letter  R :  location  vectors  are  subscripted  with  a  target-point  indicator,  location-vector 
coordinates  are  first  subscripted  with  a  coordinate  indicator  (x,  y,  or  2)  and  then  with  a  location-point 
indicator  (index  or  5).  Whenever  the  letter  R  is  used  as  a  vector  coordinate,  it  is  an  Earth-system,  rec¬ 
tangular  coordinate. 

For  both  the  gradient  and  Poynting-vector  approaches,  the  most  natural  formulation  is  in  spherical 
coordinates.  Figure  8.4  illustrates  the  spherical-system  unit  vectors  that  we  use.  We  define  the  ith 
Poynting-vector  plane  (ith  plane)  as  the  plane  passing  through  the  radar  point  that  is  normal  to  the  ith 
Poynting  vector.  The  spherical  unit  vectors  are  defined  relative  to  this  plane: 

h  is  perpendicular  to  the  ith  plane  and  positive  away  from  the  target, 

m  is  the  intersection  of  the  ith  plane  and  the  x  -  y  plane  of  the  Earth-axis  system,  and 

/  is  in  the  ith  plane  such  that  /  -  m  x  h  and  is  positive  downward. 

From  these  definitions  we  see  that  m  indicates  the  azimuth  direction  and  /  indicates  the  (negative) 
elevation  direction.  We  note  that  h  is  also  aligned  with  the  target’s  Poynting  vector  and,  equivalently, 
with  the  phase-front  gradient. 

8.2  Polarization  Effects 

The  electric  field  polarization  is  an  important  parameter  in  determining  glint  but  it  has  not  been 
explicitly  addressed  in  previous  work.  The  one-dimensional  glint  problem  has  always  been  analyzed 
with  the  implicit  assumption  that  the  target  elements  radiate  with  the  same  polarization  as  the  radar. 
This  is  a  re^onable  assumption  in  many  practical  problems.  However,  if  this  assumption  is  not  valid, 
then  current  glint  models  can  yield  erroneous  results.  For  example,  consider  the  classic  problem  of  a 
target  composed  of  two  scatterers  (Fig.  8.2  with  N  “  2).  If  one  of  the  scatterers’  reradiation  polariza¬ 
tion  is  horizontal  and  the  other’s  is  vertical,  then  an  ideal  horizontally  or  vertically  polarized  radar  will 
not  sense  one  of  the  scatterers  and  there  will  be  no  nontrivial  glint.  The  physical  components  of  a  US 
will  often  exhibit  multiple  scattering  and  such  scattering  causes  cross  polarization  [531. 
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Fij.  8.4  —  Spherical-coordinate  unit  vectors 


Implicit  in  the  glint  concepts  are  different  assumptions  regarding  the  received-field  polarizations. 
The  gradient  concept  requires  knowledge  of  the  phase  of  the  received  field.  This  phase  can  be  obtained 
accurately  by  assuming  that  the  received  fields  from  the  USs  are  plane  waves  which  are  identically 
polarized.  The  Poynting-vector  concept  requires  knowledge  of  the  Poynting-vector  direction  and  there¬ 
fore  only  assumes  that  the  received  fields  are  planar;  the  received-field  polarizations  and  directions  of 
propagation  must  be  accounted  for  in  this  method  because  they  contain  the  directional  information. 
For  long  ranges,  the  two  methods  produce  results  that  are  equivalent  for  practical  purposes.  However, 
we  prefer  the  Poynting-vector  method  because  it  accounts  for  polarization  effects,  it  is  more  closely 
related  to  physical  concepts  (power  flow),  and  it  does  not  require  knowledge  of  rate-of-change  proper¬ 
ties  of  the  US  or  multipath  functions  (Sec.  8.3). 

8.3  Phase-Front  Gradient 

From  Eq.  (4.5),  the  magnitude  of  the  received  field  has  the  form 

£  -  g(R)T(pr,  P*)  -  g(RVpi  |l  A)J P„-  (8-7) 

Equation  (8.7)  represents  only  that  portion  of  the  received  field  that  is  sensed  by  the  radar.  Assuming 
plane  wave  propagation,  the  vector  representations  of  the  fields  sensed  by  the  radar  are  simply 

i ?  *  EpK  (8.8a) 

and 

H  —  -  —  plR  (8.8b) 

V 

where  17  (—  377  fl)  is  the  characteristic  impedance  of  free  space  and  p[R  is  orthogonal  to  pR  such  that 
when  the  received  polarization  is  horizontal,  in  the  target-cg  axis-system  p  —  m  and  p  —  /.  Given  the 
large-range  assumption,  this  formulation  yields  accurate  estimates  of  the  magnitude  and  phase  of  the 
received  field. 
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and 


b<t> 

br 


N  N 

II 

J-LH 


A,Aj  -jf-  cos  ‘4,  -  tj)  -  A,-jj-  sin  (0,  -  0y) 


N  N 
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(8.14) 


cos  (0,  —  <f>j) 


The  a,  functions  C4,  —  I/1,2!  a,)  depend  on  US  orientation,  not  on  range.  Therefore  bajbr  —  0  for  all 
i  and 


N  N 


■■  -  bd>,  N  n  a|/r| 

X  X  A/Aj  .  cos  (0,  —  4>j)  -  £  i  2a^  —  sin  (0/  —  0y) 

9d>  _  i-i  >-i _ ^ _ /-i  j- 1 _ _ 

dr  ”  *  iv 

I  I  A,Aj  cos  (0,  -  <f>j) 
i- 1  ;-i 


(8.15) 


The  resulting  glint  equations  are 


N  N 


C 


90;  bAj 

A,Aj  -g—  cos  (0/  -  0,)  -  i4,  sin  (0,  -  0y) 


AT  N 


and 


G, 


90  9|/t2| 

I  I  [4  4  -g7  cos  foi  -  -  aJ  Ai  —frT  sin  (0/  -+j) 


n  n  [  9 <*.  9,4, 

II  «»  “  Ai  m  sin  (<*>i  ~  <t>j 3 


We  recall  that 


*  Af  90,  9|f/2| 

II  a‘aj  cos  4  "TT"  sin  <♦/  - 

/-i  y-i  or  or 


0/(*„9,)  -  2fc|/?,|  +y(R„9l)  ,(9,) 


(8.16a) 


(8.16b) 


(8.17) 


where  y(7?M  9,)  is  the  phase  due  to  multipath  and  /M0/)  is  the  phase  due  to  US  orientation.  Clearly, 
specific  knowledge  of  the  US  functions  is  required  to  use  the  glint  equations  in  the  form  of  Eqs.  (8.16). 
In  particular,  the  derivatives  of  the  US  and  multipath  coefficient  amplitudes  and  phases  (with  respect  to 
azimuth  and  elevation)  are  required. 


For  the  classic  problem  of  N  sources  on  a  rigid  rod  (Fig.  8.2),  various  assumptions  are  made: 
ranges  are  large,  the  scatterer  amplitudes  and  phases  are  constant  with  respect  to  rotations,  and  there  is 
no  multipath.  These  assumptions  yield 


90/  9  lit/ 1 

90  90  ’ 

90/  9  Iff,  I 

99  99  ’ 


(8.18a) 

(8.18b) 


(8.18c) 


|F/|-  1, 


(8. 1 8d) 


and 


bAj  bAj  9  If)  I 
90  99  br 


(8.18e) 
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So  that 

I  T"7  COS  (</>,  -  4>j ) 

r  _  <-l  J- 1  qij< _ 

/V  V 

£  £  cos  to,  -  <£,) 

/-I  j- 1 

We  see  from  Fig.  8.2  that  the  large-range  assumption  implies  that 

<(>,  -  <ftj  =■  2k(Lj  -  Lj)  cos  t/i 


(8.19) 


(8.20a) 


so  that 


|/?^|  =  Rs  +  Lj  cos  ill 

8l*,l  ,  .  , 

~it~  =  -L‘s  ,n* 


(8.20b) 


(8.21) 


where  the  L,  are  the  signed  distances  of  the  scatterers  from  the  target  reference  point.  (We  note  that 
Howard  [12]  incorrectly  uses  the  unsigned  distances  in  his  Eq.  (5).)  Using  Eqs.  (8.20a),  (8.20b),  and 
(8.21)  in  Eq.  (8.16a)  we  obtain  the  well-known  Eq.  (8.2). 


8.4  Poynting-Vector  Direction 

Dunn  and  Howard  [43]  showed  that  glint  was  equivalent  to  the  direction  of  the  time-average 
Poynting  vector.  The  time  averaging  is  performed  over  the  period  T,  where  o>  —  In  I T  is  the  radar’s 
frequency.  The  averaging  is  valid  because  the  phase  terms  of  the  received  field  (Eq.  (8.7))  are  approxi¬ 
mately  constant  over  the  period  T.  The  time-average  Poynting  vector  is  defined  by 


2 


For  an  /V-source  target  and  no  multipath. 


g{R)pj  £  X(R/,  0/)  Pr  Pr- 
!-\ 


(8.22) 


(8.23) 


Given  that  the  polarization  of  both  the  transmitter  and  the  receiver  are  horizontal,  in  the  target-cg 
axis-system  the  electric  field  and  the  Poynting  vector  can  be  represented  in  the  forms 


E  -  g(R)  £a,ey  '  ms 
/-] 


(8.24a) 


N  N 

II  cos  to,  -  <t»k) \ns 

/-i  *-i 


(8.24b) 


where  the  a,  and  terms  are  obtained  by  taking  into  account  the  associated  US-radar  geometries. 
(Note  that  the  pT  and  pK  vectors  must  be  expressed  in  each  of  the  individual  US  axis-systems.)  Equa¬ 
tion  (8.24b)  is  not  in  a  useful  form  because  «,  is  not  known  directly.  Therefore,  we  next  derive  an 
equation  for  P  in  terms  of  the  N,  known  Poynting  vectors. 

Equations  for  the  locations  of  the  USs  were  derived  in  Sections  5  and  6;  given  the  radar  point  and 
the  US  locations,  the  Poynting  vector  for  each  US  can  be  determined.  The  total  fields  due  to  the  USs 


t  -  g(R)  £  a,eJ*<m, 
/-! 


(8.25a) 
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and 


H  =  -  £  a,^*% 


and  the  associated  Poynting  vector  is 

g2(R)  n  N 

p - - -  £  £  a, a*  cos  (0,  -  <t>k)m,  x  4. 

ZTJ  1-1  *-1 

We  define  the  position  of  the  /th  scatterer  relative  to  the  radar  (in  Earth  coordinates)  as 

R,  =  RXli  +  Ry)  +  Rzfi 

where  for  clarity,  the  R,  time-dependence  is  suppressed.  In  the  /th  US-system 


n,  “  - 


1 


Uty  +  V  +  V>- 

Because  m(  is  in  the  jc->>  plane  of  the  Earth-axis  system  and  is  perpendicular  to  Rh 

fhi  =  sin  «/»,  /  —  cos  <fij 

where  t/t,  —  tan-1  (Ryi/Rx)  and  Rx^  =  /?Xj/  +  /?y.y  (Fig.  8.5).  Finally, 


«  Rt, 

l,~  m,x  n,  -  j=zj  cos  «/»,/  +  j=j  sin  «/»,;  - 


sin  /?x  cos  1)1  i 


\R< 


I/?, I 


k. 


(8.25b) 

(8.26) 

(8.27) 

(8.28a) 

(8.28b) 

(8.28c) 


Fig.  8.5  —  The  /th  azimuth  and  elevation  angles 
in  the  Earth-axis  system 


By  using  Eqs.  (8.28),  the  cross  product  term  of  Eq.  (8.26)  is 

mt  x  4  —  ntj  x  ( mk  x  hk). 
Expanding  Eq.  (8.29a)  in  terms  of  known  quantities, 


(8.29a) 


Ry  sin  0*  +  Rx  cos  0* 
m,  x  4  - - r=jn - cos  <M 


1**1 

K,*  sin  0*  +  /?X|k  cos 


sin  tM  +  -j-=-r  cos  (t pk  -  t \>i)k. 


(8.29b) 


Equations  (8.28)  express  the  /th-system  unit  vectors  in  (rectangular)  Earth  coordinates.  Referring  to 
Fig.  8.5  we  see  that  the  transformation  from  Earth  coordinates  to  target-cg  coordinates,  Tet , 's 


Tet- 


\R,\ 


Rt t  cos  \f/s  Rz ,  sin  ips  -Ry>  sin  i/ts  -  RXj  cos  0, 
sin  </>j  -cost (i,  0 


R r 


(8.30) 


.’.liUi 
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I 


I 


where  the  subscript  s  indicates  the  target’s  eg  location  and  *jis  —  tan-1  (RyJRx).  The  Poynting  vector 
is  now 


P  ”  ~  I  I  a‘ak  cos  ~  <f>k )  TET(nti  x  \)t 

zrl  I  /-i  *-i 


g2(R)  T 

“  2tj 


A  JV 

II  cos  (<f>t  —  x  lk)e 

t- 1  *-i 


(8.31a) 

(8.31b) 


where  the  ( m ,  x  4)  vector  is  expressed  in  Earth  coordinates.  The  azimuth  and  elevation  glint  can  now 
be  expressed  as 


I R,\ 


N  N 

JII  a,ak  cos  (<t>j  -  (f>k)  TET(m,  x  lk)e I  •  m, 
|/-i  *-1 

N  N 

JII  a,ak  cos  (<f>,  -  </>*)  TET(mt  x  /fc)e| 

/-t  *-1 


(8.32a) 


G,=  \R,\ 


N  N 

II  atak  cos  (0,  -  <(>k)  TET(m,  x  lk)t 

/- 1  k—  1 _ 

A  A 

II  a,ak  cos  (^,  -  <(>k)TET(mi  x  4), 
/-1  *-1 


I-4J 


(8.32b) 


We  note  that  for  small  US-separations  and  long  ranges,  m(  x  4  =  n,  so  that  the  denominator  of  Eqs. 
(8.32)  is  approximately  the  target’s  RCS  (Eq.  (7.7a)). 


8.5  Multipath  Effects 

For  low-grazing-angle  illumination  of  ships,  over-water  multipath  affects  the  received  signal.  In 
Section  4,  we  described  a  multipath  model  appropriate  for  point-source  scatterers.  In  this  section,  we 
discuss  the  modeling  of  the  effects  of  multipath  on  glint. 

A  different  method  of  multipath  modeling  is  required  for  each  glint  concept  because  they  are 
developed  under  different  assumptions.  The  phase-front  gradient  method  requires  explicit  knowledge 
of  only  the  received  amplitudes  and  phases  whereas  the  Poynting-vector  method  also  requires  explicit 
direction-of-propagation  information. 


8.5.1  Phase-Front  Gradient 


For  the  phase-front  gradient  model,  we  view  multipath  as  an  alteration  of  the  received  amplitude 
and  phase  from  the  US.  Referring  to  Eq.  (4.17),  we  rewrite  Eq.  (8.9)  to  explicitly  show  the  multipath 
effects: 


4> 


tan  1 


j^\F,2\  a,  sin  <t>, 
1- 1 _ 

N 

£  If,2 |  a,  cos  d>i 

I- 1 


(8.33) 


where  <£,  is  defined  by  Eq.  (8.17).  From  Eqs.  (8.16)  we  see  that  the  terms  of  interest  are  the  partial 
derivatives  of  \F2\a,  and  <t>,  with  respect  to  r,  1 (1,  and  9  (in  the  target-cg  axis  system).  Because  there 
are  in  general  no  physical  interpretations  for  these  rate-of-change  processes,  they  would,  at  best,  be  dif¬ 
ficult  to  estimate.  They  are  also  analytically  complex  so  we  will  develop  the  general  multipath  formula¬ 
tion  only  for  the  Poynting-vector  method,  which  yields  physically-meaningful  results.  However,  we  will 
analyze  the  special  case  of  very  calm  (smooth)  seas. 
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where,  from  Fig.  8.6,  Rh  -  8,  -  8,.  Now 

( &Z,  +  1 

a,  =  288,  — p=- — 


so  that 


=  288,  -r~  +  cos  9, 


(8.44) 


(8.45) 


f  8,  8 1/?, I  88/) 

do  ~~  2kh*  f  ~W  ~ sm  e‘  ai"J-  (8  45) 

Because  sin  0,  =  1,  8 9t/d9  =  1,  8,/|/i,|  «  1,  and  d|/?,|/80  «  |/?,|,  we  have  that  for  a  smooth 
sea-surface, 

9a, 

*  _  A  i  i  /«  a  « v 


VU  / 

IT""2***' 

Hence,  the  multipath  components  of  the  phase-front-gradient  formulation  are 

8 IF)2 1  _  4*8,8,  [288,8, 

8r  *  |^|2  C0S|  |^|  ’ 

8r,  _  88,8, 

•'  ■“  iv* 

so  ikh*  005  |  ■ 


(8.46) 


(8.47a) 

(8.47b) 


(8.47c) 


(8.47d) 


(8.47e) 


hi  88 

dg  - 


(8.470 


8.5.2  Poynting-Vector  Direction 

To  introduce  multipath  into  the  Poynting-vector  formulation,  we  modify  the  received  fields  by  the 
multipath  factor,  F,2,  described  in  Sec.  4.4  (Eqs.  (4.18)).  This  factor  modifies  the  direct-return  scat¬ 
tered  fields  and  induces  N  apparent  (image)  scatterers.  This  model  allows  the  obvious  extension  of  Ea. 
(8.31b)  to 

-  e2(R)  \2N  2N 

Pm,~  *  2ri  Tet  2}  L  M,Mka,ak  cos  ( <t> ,  -  $*)(m,  x  /*),  (8.48) 
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where 

j  Fj  ,  1  <  i  <  N 

M‘  “  -  1)  ,  N +  l<  l  <2N, 

/  <  N,  m  -  1  +  N, 

m,  -  mm 

and,  from  Fig.  8.6, 


R,  -  R;  -  RX1  +  Ry  j  -  (hR  +  hj)ic,  N  +  1  <  /  <  IN. 


Fig.  8.6  -  Multipath  geometry 


(8.49a) 

(8.49b) 

(8.49c) 


Equation  (8.49c)  implies  that  Eq.  (8.28c)  becomes 

.  v  .  R>:  .  , 
'■-W\cos*"‘  +  lZ\sm*lJ 


Ry,  sin  ifii  +  RX(  cos  0, 

W\ 


(8.49d) 


for  N  +  1  <  i  <  27V.  Therefore  when  accounting  for  multipath,  the  glint  equations  become 


Ga  -  I*,  I 
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(8.50b) 


2IS  2IS  * 

II  MiMkdjdk  cos  (0,  -  <t>kf^ET^mi  x  k>e\  '  ns 


As  in  Eqs.  (8.32),  for  small  US-separations  and  long  ranges  the  denominator  of  Eqs.  (8.50)  approxi¬ 
mates  the  target’s  RCS. 

It  is  clear  that  Eqs.  (8.50)  are  more  natural  and  accurate  for  use  in  modeling  the  glint  of  A-source 
targets  than  are  Eqs.  (8.16).  Equations  (8.50)  are  used  in  the  simulations  (Appendix  C). 
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9.  SIMULATION  RESULTS 

In  this  section  we  briefly  summarize  the  results  of  simulations  that  were  performed  by  using  the 
method  of  Section  6  to  simulate  US  motion  and  the  equations  of  Sections  7  and  8  (as  modified,  for 
efficiency,  in  Appendix  C)  to  simulate  RCS  and  glint.  The  source  code  that  was  used  is  presented  in 
Appendix  D.  Because  of  the  number  of  parameters  affecting  the  scattering  process,  a  parametric  study 
of  the  model  was  beyond  the  scope  of  this  report.  Among  those  parameters  are: 
ship  structure 
ship  maneuvering 
ship  speed 

ship  response  to  the  sea  (RAOs) 
sea  state 

sea  structure  (sea,  swell,  etc.) 

ship  heading  relative  to  the  dominant  wave  direction 

radar-to-ship  aspect  angle 

radar-to-ship  elevation  angle 

radar  frequency 

radar  pulsewidth 

relative  motion  between  the  radar  and  the  ship 
US  amplitudes 
US  phases 
US  locations. 

We  note  that  when  measurements  of  ship  scattering  (e.g.,  RCS  or  glint)  are  made,  measurements  of 
the  first  twelve  parameters  listed  should  also  be  made. 

In  the  simulation  studies,  we  have  investigated  the  qualitative  effects  of  variations  in  ship 
maneuvering,  radar-to-ship  aspect  angle,  and  multipath  (which  is  primarily  a  function  of  radar-to-ship 
elevation  angle,  radar  frequency,  and  sea  state). 

9.1  Simulation  Parameters 

In  the  simulations,  the  radar’s  wavelength  was  3  cm  and  the  radar  was  located  at  a  fixed  point  that 
was  40  m  above  the  sea  surface.  The  ship’s  position  was  initialized  at  10  km  from  the  radar.  This 
yields  grazing  angles  on  the  order  of  0.25°. 

Two  maneuvers  were  considered  for  the  ship:  straight-line  motion  and  turning.  For  the  straight- 
line  motion  the  ship’s  speed  was  10  kn,  and  for  the  turning  motion  the  turn  rate  was  0.4°/s  and  the 
ship’s  speed  was  6  kn.  The  maneuvers  were  simulated  with  and  without  multipath  effects.  For  the 
straight-line  motion  an  rms  wave  height  of  0.762  m  was  used  because  this  corresponds  to  the  sea  state 
associated  with  the  ship-motion  spectra  that  were  used  (Appendix  A).  For  the  turning  maneuver  an 
rms  wave  height  of  0.135  m  was  used  and  the  ship-motion  spectra  were  changed  from  the  spectra  of 
Appendix  A  by  reducing  the  rms  values  by  a  factor  of  10  for  each  of  the  six  motion  processes;  this  was 
done  to  approximate  the  ship’s  motion  in  slight  to  moderate  seas.  For  the  geometry  of  the  simulations 
the  sea-roughness  factor  (a  fundamental  parameter  characterizing  multipath,  defined  by  Asm  tfi/k  where 
h  is  the  rms  wave  height,  0  is  the  grazing  angle,  and  X  is  the  radar’s  wavelength;  see  Ref.  49)  is  on  the 
order  of  0.025  for  the  turning  motion  and  0.133  for  the  straight-line  motion.  The  simulation  time- 
increment  was  0.002  s,  and  the  simulated  time  was  2. 1  s. 

The  USs  used  in  the  simulation  were  chosen  based  on  visual  identification  of  scatterer  type  and 
approximate  projected  area.  For  X  —  0.03  m  the  sphere-like  (constant  RCS)  scatterers  have  amplitudes 
that  are  small  relative  to  the  flat-plate  and  comer-reflector  types,  assuming  comparable  projected  areas. 
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Therefore  only  flat-plate  and  corner-reflector  types  were  used  in  the  simulations  because  such  scatterers 
dominated  the  ship’s  structure.  A  ship  was  chosen  for  which  photographs  were  available  and  the  USs 
were  chosen  for  three  aspect  angles:  broadside,  aft-quarter,  and  stern.  Table  9.1  lists  the  location,  mag¬ 
nitude,  and  type  (flat-plate  is  type  2  and  corner-reflector  is  type  3)  for  the  USs  that  were  used.  Figure 

9.1  is  intended  as  an  aid  in  visualizing  the  ship  representation:  it  contains  two-dimensional  plots  of  the 
US  locations  relative  to  the  ship’s  eg  using  the  location  data  from  Table  9. 1 . 

9.2  Discussion 

The  output  of  each  simulation  run  is  presented  as  a  figure  composed  of  four  plots:  RCS,  azimuth 
glint,  and  elevation  glint  as  functions  of  time,  and  azimuth  glint  vs  elevation  glint  for  each  time  instant 
(connecting  the  points  from  instant  to  instant).  RCS  is  plotted  in  dB  above  one  square  meter  (dBsm), 
and  all  RCS  plots  use  the  same  scale.  The  scales  for  the  glint  plots  vary  (due  to  the  nature  of  glint). 
These  plots  are  Figs.  9.2  to  9.15. 

Figures  9.4  through  9.15  display  the  output  of  twelve  simulations.  The  parameters  varied  in  these 
simulations  were  maneuvering  (straight-line,  turn);  aspect  (broadside,  aft-quarter,  stern);  and  multipath 
(yes,  no) .  Although  the  parameter  variations  are  limited,  the  outputs  indicate  the  qualitative  effects  of 
the  parameters  on  RCS  and  glint.  In  the  remainder  of  this  section  we  discuss  some  of  the  most  salient 
features  of  the  output  data. 

9.3  Multipath  Effects 

Figures  9.2  and  9.3  illustrate  the  effects  of  (point- source)  multipath  for  the  two  multipath  condi¬ 
tions  chosen.  In  each  case  the  only  motion  is  that  of  the  sea  surface;  the  ship  scatterers  do  not  move. 
We  note  the  following  implications  of  the  sea-roughness  factor:  the  rms  waveheight  of  0.762  m  implies 
that  diffuse-component  multipath  effects  (random)  dominate  those  due  to  the  specular  component  and 
the  rms  waveheight  of  0.135  m  implies  that  the  specular  effects  (deterministic)  dominate.  We  also  note 
that  the  diffuse  multipath  effects  are  essentially  independent  (except  through  variations  in  the  value  of 
the  roughness  factor)  of  ship  motion  for  our  scenarios,  but  the  specular  multipath  effects,  which  are 
deterministic,  are  of  course  correlated  with  that  motion.  Clearly,  both  specular  and  diffuse  multipath 
modulate  the  scattered  fields  further  complicating  the  problem  of  US  identification  in  the  presence  of 
multipath. 

9.4  RCS 

It  is  clear  that  the  correlation  properties  of  the  RCS  are  strongly  dependent  on  aspect  angle.  We 
see  this  in  Figs.  9.4(a)  through  9.15(a)  and  it  is  quantified  by  the  normalized  spectra  shown  in  Figs. 
9.16  through  9.18.  (These  spectra  were  calculated  using  the  TSAP  program,  which  is  described  in  Ref. 
96.  A  time  series  of  1024  points,  2.048  s,  was  used.  First  the  autocovariance  function  was  calculated 
using  a  12.5%  lag.  This  result  was  windowed  and  then  an  FFT  was  used  to  obtain  the  spectral  esti¬ 
mates.) 

Examining  the  quartering-aspect  RCS  in  Fig.  9.5(b)  we  see  that  as  shown  in  Chapter  7,  even 
assuming  constant  amplitude  functions  (Table  9.1(b))  the  slow,  oscillatory  motions  of  the  ship  can 
yield  RCS  realizations  that  appear  nonstationary  over  short  time-intervals.  (Note  that  there  is  an 
apparent  decreasing  trend  in  the  variance  of  that  data.) 

9.5  Glint 

From  the  data  we  see  that  glint  is  a  process  with  a  wider  bandwidth  than  the  RCS  process.  Figure 
9.19  contains  examples  of  glint  spectra  (calculated  by  using  the  procedure  described  in  Sec.  9.3)  for 
data  (from  Figs.  9.5  and  9.8)  that  exhibit  no  obvious  nonstationary  characteristics  in  the  mean  or  vari¬ 
ance. 
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Fig.  9.8  —  Straight-line  motion,  ship  speed  “  10  kn,  aft-quarter  aspect,  rms  wave  height  —  0.761 
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Fig.  9.9  —  Straight-line  motion,  ship  speed  -  10  kn,  stern  aspect,  rms  wave  height  -  0.762  m 
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Fig.  9.12  -  Ship  turn  rate  -  0.4°/s,  stern  aspect,  no  multipath 


Fig.  9.13  -  Ship  turn  rate  -  0.4°/ s,  broadside  aspect,  rms  wave  height  -  0.135  m 
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Fig.  9.14  —  Ship  turn  rate  —  0.4°/s,  aft-quarter  aspect,  rms  wave  height  —  0.135  m 


Fig.  9.15  —  Ship  turn  rate  —  0.47s,  stern  aspect,  rms  wave  height  —  0.135  m 
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Fig.  9.16  —  RCS  spectra  for  straight-line  motion,  ship 
speed  “  10  kn,  no  multipath 

Fig.  9.17  —  RCS  spectra  for  straight-line  motion,  ship 
speed  “  10  kn,  rms  wave  height  ~  0.762  m 
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Fig.  9.18  —  RCS  spectra  for  turning  motion,  ship 
speed  —  6  kn,  rms  wave  height  —  0.135  m 
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Fig.  9,19  —  Glint  spectra  (straight-line  motion,  ship  speed  -  10  kn,  aft-quarter  aspect) 
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A  well-known  property  of  glint  is  also  illustrated  in  the  data:  sharp  decreases  in  the  magnitude  of 
RCS  ("fades")  produce  sharp  increases  ("spikes")  in  the  magnitude  of  glint.  Recalling  from  Section  8 
that  glint  is  the  ratio  of  two  functions  where  the  denominator  approximates  RCS,  we  see  that  this  fade- 
spike  correlation  is  to  be  expected. 

If  Figs.  9.4(b)  through  9.15(b)  are  interpreted  as  scatter  diagrams,  then  we  see  that  azimuth  and 
elevation  glint  are  in  general  not  strongly  correlated.  Only  the  stern  aspects  show  strong  correlation 
(due  to  the  associated  scenario  producing  very  little  signal  fluctuation).  To  further  confirm  this  obser¬ 
vation,  cross-covariance  functions  for  the  glint  data  of  Figs.  9.4,  9.5,  9.7,  and  9.8  were  calculated  (using 
TSAP)  and  are  plotted  as  Fig.  9.20.  We  suspect  that  the  relatively  strong  but  brief  correlation  that  does 
occur  (at  zero  lag)  is  due  to  the  simultaneous  occurrence  of  spikes.  We  also  note  that  multipath 
appears  to  contribute  to  decorrelation  between  azimuth  and  elevation  glint. 
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Fig.  9.20  —  Cross-covariance  functions  for  the  glint  data  of  four  simulations 
(straight-line  motion,  ship  speed  -  10  kn) 


An  examination  of  the  azimuth-glint  plots  shows  that  the  data  are  often  asymmetrical  about  a 
mean  near  zero  (e.g.,  in  Fig.  9.4(c)  the  positive  values  are  more  correlated  (smooth)  than  the  negative 
values).  This  suggests  that  such  data  may  not  be  Gaussian  distributed.  (It  is  known  that  some  radar 
targets  display  glint  statistics  that  are  approximately  Gaussian  [97].)  This  observation  was  confirmed  by 
using  a  Kolmorgorov-Smirov  test  of  fit  (from  TSAP)  to  test  the  hypothesis  that  the  azimuth-glint  data 
is  Gaussian  distributed:  for  each  time  series  the  result  was  rejection  of  the  hypothesis  at  an  a -level  less 
than  0.01.  This  result  also  holds  for  the  elevation-glint  data  except  for  one  time  series:  that  in  Fig. 
9.10(d).  We  comment  here  that  what  we  simulated  was  "ideal"  glint  in  the  sense  that  it  cannot  be 
exactly  measured  by  radars.  Constraints  on  radar  system  performance  (antenna  pattern  effects,  finite 
bandwidths,  system  noise,  etc.)  may  result  in  error  signals  that  are  more  symmetrical  (due  to  smooth¬ 
ing  of  sharp  variations)  and  hence  more  likely  to  appear  Gaussian. 
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Figure  9.21  illustrates  the  correlation  between  RCS  and  azimuth  glint  for  four  of  the  simulations. 
We  note  that  a  weak  but  definite  correlation  exists  between  these  two  processes. 
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(c)  Aft-quarter  aspect,  straight-line  motion,  ship  speed  -  10  kn  (d)  Aft-quarter  aspect,  turning  rate  -  0.47s,  ship  speed  «*  6  kr 


Fig  9.21  —  Normalized  cross-covariances  between  RCS  and  azimuth  glint  in  the 
absence  of  multipath 


10.  CONCLUSIONS  AND  FURTHER  RESEARCH 


The  major  result  of  this  research  is  the  development  of  a  unified  stochastic  model  and  associated 
simulation  for  the  EM  scattering  from  complex  objects  in  the  presence  of  overwater  multipath.  By 
modeling  each  component  of  this  very  complex  problem  using  an  N-source  formulation,  we  have 
obtained  a  unified  RCS-glint  model  that  is  useful  in  applications.  We  have  introduced  the  concept  of 
unit-scatterer  and  shown  that  it  is  a  generic  element  for  use  in  modeling  scattering  effects.  Though  by 
defining  the  unit-scatterer  as  our  generic  scattering  element  we  have,  in  a  sense,  only  shifted  the 
scattering  problem  from  the  entire  object  to  sub-parts  of  the  object,  we  suggest  that  the  unit-scatterer 
concept  is  fundamental  to  the  practical,  pulse-by-pulse  modeling  and  simulation  of  scattering  from  com¬ 
plex  objects  We  believe  that  it  may  not  be  necessary  to  precisely  characterize  unit-scatterers  to  obtain 
good  approximations  to  the  total  scattering  If  this  is  the  case,  then  having  determined  classes  of  US 
types  one  could  quickly  build  models  and  efficient  simulations  for  the  scattering  from  complex  objects 
by  using  elements  from  these  classes  as  building  blocks. 

We  suggest  four  major  areas  of  research  for  extension  of  this  work: 

(1)  The  primary  area  concerns  the  identification  and  characterization  of  unit-scatterers. 
Particularly  useful  (and  cost  effective)  would  be  simulation  studies  using  the  geometric 
simulations  such  as  those  at  the  NRL  and  the  Georgia  Institute  of  Technology.  It  would  be 
important  to 
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(a)  identify  US  types, 

(b)  characterize  US  amplitude  functions,  and 

(c)  characterize  US  phase  functions  with  emphasis  on  determining  the  rate-of-change  of  phase 
with  US  rotation. 


(2)  Analysis  of  existing  ISAR  data  for  identifying  and  characterizing  unit-scatterers  as  suggested 
in  the  first  area. 

(3)  Second-order-statistical  comparisons  of  simulation  data  with  measured  data  where  the 
simulation  data  is  generated  using  the  parameters  associated  with  the  measured  data. 

(4)  A  study  of  the  effects  of  multipath  on  objects  that  are  distributed  to  determine  a  better 
model  for  the  multipath  effects  on  unit-scatterers  than  the  existing  point-source  models. 
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Appendix  A 

AN  EXAMPLE  SIMULATION  OF  SHIP  MOTION 

A  simulation  of  the  6  dof  motion  for  a  ship’s  eg  is  described  in  Section  6.  This  appendix  presents 
the  output  of  a  simulation  that  was  implemented  according  to  that  description.  The  ship  chosen  for 
simulation  is  the  DD963  because  that  is  the  ship  used  as  the  example  in  Ref.  78.  A  ship  speed  of  10 
knots  and  a  relative  ship-to-dominant-wave  heading  of  30°  were  chosen  for  the  simulation.  Figure  A1 
illustrates  the  spectra  calculated  by  DTNSRDC  for  the  DD963  at  10  knots.  The  center  frequency, 
bandwidth,  and  rms  values  used  in  the  simulation  were  estimated  from  this  figure  and  are  listed  in 
Table  Al.  The  frequency  increment,  Aw,  was  chosen  to  be  constant  and  equal  to  BW/10  for  each  com¬ 
ponent  process.  The  simulation  runs  used  a  time-step  size  of  1.0  s  in  generating,  six,  2048-point  sam¬ 
ples.  The  first  300  points  of  each  of  these  samples  are  plotted  as  Fig.  A2.  The  samples  were  analyzed 
using  the  time  series  analysis  package  described  in  Ref.  79  to  estimate  the  associated  spectra  and  rms 
values.  The  resulting  rms  values  are  tabulated  in  Table  A2  and  normalized  spectral  estimates  are  plotted 
in  Fig.  A3.  Figure  A4  illustrates  2-s  samples  of  the  motion  processes  generated  by  using  a  time-step 
size  of  0.002  s. 


SWAY  ly) 


Fig.  Al  —  DTNSRDC  motion  spectra  for  the 
example  simulation 


(a)  Angular-position 
spectra  (deg) 


n  (b)  Linear-position 
spectra  (ft) 

Table  Al  —  Spectrum  Parameters  That  Were 
Input  to  the  Example  Simulation 


Process 

wo  (Hz) 

BW  (Hz) 

RMS  Value 

X 

0.061 

0.018 

0.344 

Y 

0.076 

0.029 

0.357 

Z 

0.069 

0.025 

0.411 

PHI 

0.083 

0.022 

0.0539 

THETA 

0.065 

0.018 

0.0089 

PSI 

(normalized)  of  the  simulated  cg-motion  processes 


Appendix  B 

AN  ESTIMATE  OF  RANGE  VARIATION  DURING 
A  SHORT  TIME  PERIOD 


Let  Rj(t)  be  the  range  from  the  radar  to  a  scatterer  on  the  ship.  The  phase  of  the  received  field 


due  to  range  is 


<M')  =  !*/('> I- 

To  estimate  the  variation  in  |R,  (f)|  over  short  periods  of  time,  we  use  the  measure 


£{AR,2}1/2  =  £{(|£,(r  +  r)  -  £,(r)|)2}1/2. 


From  Eq.  (6.5b), 


bR,  =  | rs(f  +  r)RfU  +  r)  +  Rf(t  +  r)|  -  \Ts(t)Ri(t)  +  Rf(t)l  (B3) 

For  simplicity  we  assume  that  there  is  no  maneuvering,  hence  Ts(t)  =  /  and  RfU)  “  Rf( 0).  There- 


From  Eq.  (6.5), 


^R:  -  \Rf(t  +  r)  +  Rf( 0)|  -  'Rf(t)  +  Rf( 0)|. 


A£,  =  | Rj (0)  +  7;{3cc*(r  +  r)  -  xc,(0)]  +  iff (0)1 


-  |  R,  (0)  +  7;  [*<*(/)  -  *c*(0)]  +  j?f(0)|. 


Collecting  the  initial  terms, 


A/?,  -  | R?  +  T^U  +  r)|  -  I R?  +  7}xc*(f)| 
<  I  T,  [xfg  (/  +  t)  -  xcg  (/)  1 1. 

Expressing  the  difference  in  Eq.  (B6)  in  terms  of  components, 

A/?,-  <  1 7)  (Ax,  Ay,  A r,  A <f>,  AS,  0<|r]r I 


A R,  <  [(Ax  +  z, A0  -  y,A</»)2  +  (Ay  -  z,A4>  +  x,A« p)2  +  (A z  +  y,A<£  —  x,A0)2]l/2 


‘i* 
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so  that 

£|  A R2  <  Ax2  +  (ztA9)2  +  0», A ifi)2  +  A>>2  +  (ztA<t>)2  +  (x,Atp)2 

+  Az2  +  (y,A<f>)2  +  (x,A#)2  +  cross-product  terms.  (B8) 

Because  the  component  processes  were  assumed  to  be  independent,  we  have 

£{A£(2}  -  £{Ax2}  +  £(A.y2}  +  £{Az}2  +  (y2  -I-  z,2)£{A^2} 

+  {x2  +  z2)E{A92)  +  {x2  -+-  >v2)£{A02}  (B9a) 

where  each  expectation  has  the  form 

£{Ax2}  —  £{lx(r  +  r)  —  x(r)]2}  “  2£,(0)  —  2Rx(t).  (B9b) 

As  described  in  Section  6,  a  useful  approximation  to  the  covariance  function  for  each  process  has  the 

form 

f"  Rj(r)  =  (r2e  cos  co0(r  (BIO) 

where  a  is  1/2  the  bandwidth  and  <u0(  is  the  center  frequency  of  the  ith  process.  Using  these  approxi¬ 
mations  we  have,  for  example, 

|j|  £{Ax2}  =■  £{(x(f)  -  x(t  +  r)]2}  »■  2crl[\  —  e  cos  w^t].  (Bll) 


For  a  ship,  three  scatterer  locations  yield  the  largest  values  for  the  position  coordinates:  the  loca¬ 
tion  highest  above  the  ship’s  eg  and  the  locations  at  the  bow  and  stern.  For  the  example  ship,  the 
DD963,  Xm,,  ~  50  m,  ymtx  ~  5  m,  and  zmax  =  15  m.  Because  of  ship  structure,  each  maximum  does 
not  occur  at  the  same  location.  The  broadside  aspect  presents  the  largest  range  variations  because  the 
largest  angular  rate,  roll,  has  the  maximum  effect  on  the  z-axis  components  and  the  largest  scatterer- 
distance  from  the  eg,  along  the  x-axis,  also  yields  its  maximum  range  effect.  Even  if  the  scatterer  was 
located  such  that  the  location  maxima  occurred,  we  find  that  using  the  position  values  above,  the  pro¬ 
cess  values  from  Table  Al,  and  a  typical  PRI  of  r  *  0.002  s,  £{A£2},/2  =  2.1  cm,  and  for  r  —  1.0  s, 
£{A£2),/2  =  0.42  m. 
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Appendix  C 

SIMULATION  EQUATIONS  FOR  RCS  AND  GLINT 

We  use  Eqs.  (7.7a)  and  (7.9)  to  simulate  RCS  where,  in  anticipation  of  the  glint  simulation,  we 
modify  them  to  be  in  the  form 

IN  2N 

c r  “  £  £  cos  OF,  —  V j)  (Cla) 

<-i  J- 1 


where 


Vj  -  'Ey  -  (|R, I  -  |R,I)  +  Pi(5,)  ~  /3 j(Sj)  +  y,(R„  0,)  -  yj(Rj  -  Oj)  (Clb) 

A 


\\Fk\  ,  1 

\Fk(Fk  —  1)1  ,  N  +  l^k<2N 


|<  Fk  ,  1  <  k  <  N 

r*“|<  Fk(Fk-  1)  ,  N  +  HZk^2N. 

We  describe  in  detail  the  simulation  of  the  multipath  coefficients,  Fkt  in  Refs.  49  and  98. 


(C2b) 


The  glint  equations  used  for  simulation  are  Eqs.  (8.50).  Next,  we  make  an  assumption  that  sim¬ 
plifies  those  equations  for  use  in  the  simulation. 

Because  we  have  limited  our  studies  to  small-grazing-angle  illumination  of  the  ship  and  because 
the  ship’s  eg  is  near  Earth’s  surface, 


which  implies  that  (see  Fig.  8.5) 


Equations  (C4)  imply  that 


R?  «  RJ  +  R‘ 

ii  xt  * 


14-1=  I*,  i 


I4J  =  14. 


sin  i l>s 


sin  </>, 
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cos  0,  =  ~r=rr-  (C5d) 

\Ri  \ 

By  use  of  Eqs.  (C5),  the  transformation  matrix  in  Eqs.  (8.50)  becomes 

cos  0S  cos  0S  cos  0S  sin  0S  —1 

Tet~  sin  0S  —  cos  (Jr*  0  (C6) 

cos  ips  sin  ifis  cos  0S 

and  the  cross  product  term  becomes,  from  Eq.  (8.29b), 

(m,  x  lk)e  -  cos  0,,'  +  sin  0,7  +  cos  9k  cos  (0*  -  t f>/)k.  (C7) 

Therefore, 

(rh/  x  jk)s  =■  TEj ( ntj  x  jk)e 

=  [cos  9S  cos  if)s  cos  0,  +  cos  9 s  sin  i fts  sin  0,  -  cos  9k  cos  (0*  -  0,)] 4 

+  [sin  0S  cos  i />,  —  cos  0S  sin  0,]mj 

+  [cos  05  cos  0,  +  sin  0S  sin  0,  +  cos  0S  cos  9k  cos  (0*  -  « 0,)]«j  (C8a) 

-  [COS  0S  COS  (0S  -  0,)  ~  COS  0*  COS  (0*  -  0,)]/j 

+  sin  (0*  -  0,)m,  +  [cos  (0s  -  t 0,)  +  cos  9S  cos  0*  cos  (0*  -  0/)ns.  (C8b) 

By  use  of  Eq.  (C8b)  in  Eqs.  (8.50),  we  have 

2 N  IN 

II  MidiMjdj  cos  (0,  -  <f>j )  sin  (0S  -  0,) 

C  ss  |j?  I  '-■*  ^ _ _ _  (C9a) 

~  2N  2 N 

II  MidiMjdj  COS  (0,  -  0; )  COS  (0S  -  0,) 
and 

W  21V  .. 

II  MidiMjQj  COS  (0,  -  0y )  [COS  9j  COS  (0 j  -  0,)  -  COS  0J  COS  (0S  -  0/)J 


<7,  =  |«sl 


2/V  2V 

II  MjdjMjdj  cos  (0,  -  <t>j)  cos  (0j  -  0,) 

<-i  j- 1 


iN  2N 

II  MjdjMjdj  cos  (0,  -  0y)  cos  0 j  cos  (0 j  -  0/) 

1^  I  —  ■/r1 - - - cos  0 j  . 

|n!'  2M  2Af  5 

II  MjdjMjdj  cos  (0/  —  0y)  cos  (0j  —  0,) 


/-i  y-i 


Because  sine,  cosine,  and  arctangent  functions  are  time-expensive  on  digital  computers,  we  use  Eqs. 
(C5)  in  Eqs.  (C9)  to  obtain 


RXR*,  Ry, 

' 

R*kR*.  ~  RykRy* 

lif,l 

lif,l  •  lif,l 

\Rk\ 

lif,l  •  lifkl 

Rv  Rx  -  Rx  Ry  Rx  Rx  +  Ry  Ry 

y%  xi _ xs  yi  «  ,  *1 _ ys  yJ_  * 

I  if,  I  •  I  if,  I  m*  l*,l  •  I  if,  I 
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where  we  recall  from  Section  8  that 


Rx,‘  +  Ry,J  +  ,  i  <  i  <  N 

RX1  +  Ry  j  -  (hK  +  h,)k  ,  N  +  1  <  i  <  2N 


(Cl  la) 


=  Rxi  +  Ryj  +  Rz'k. 

Therefore  the  simulation  equations  for  glint  become 

Ga  =  \RS\ 


2 N  2N 

/i  1/^  si  ( s4\  -  s4\  ) 

“  R*,Ry, 

2^  2*  MiQjMjUj  COS  \<pj  y>jJ 
i-t  j-t 

\R,\ 

IN  IN 

V  V  Kf  n  M  a  cos  (th,  —  A.) 

R*, Rxt  ~  Ry,Ryi 

7  -  7  .  ;Mi  lYAj  uj  V'U3  'Ky/ 

/-I  y-l 

\R,\ 

(Cl  lb) 


(C12a) 


Ge  =  \RS\ 

IN  2 N  RJ  RXRX,  +  Ry  Ry, 

IN  2  N  RX  Rx  -  Ry  Ry 

II  COS  (<j>i  -  0y)l - J  --■■I 


(Cl  2b) 


Rx  Rx  +  RyRy' 

'  J  '  yi  y< 


I  5  CM  <+,  -  *,>*,' 


IN  2N  RX  Rx  —  Ry  Ry 

l  l  M,  a,  Mj  Qj  cos  (<t>i  —  )  — - 


(C12c) 


Finally  we  note  that  because  of  the  long-range  assumption,  (<//s  -  i />,)  and  («(»,  -  fly)  are  small  angles  so 
that 


cos  (0  s  —  0 , )  ~  1 

COS  (\Jl,  -  l/lj)  =  1. 


(Cl  3a) 
(Cl  3b) 
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fi  " 

, 

K  . 


By  use  of  Eqs.  (Cl 3)  in  Eq.  (Cl 2a)  and  Eq.  (Cl 2b),  we  have 


IN  IN  y ,nx, 

I  I  M^MjOj  cos  (<*>,  -  4>j) - r= 

i=>i  w  y-t 

Gfl  ~  l)*jl  2V  2iV  ~  ” 

£  ^  MiaiMjOj  cos  (<t>,  -  4>j> 
-l  y-i 


(Cl  4a) 


RxjRx,  Ryfiy/ 


X  £  MiaiMjOj  cos  (<£/  4>j)Rij  |£  |  |^|2 

Ge  ~  |/fj  -  2V  2V  ~  ~  ~ 

£  £  cos  (<£,  -  <^>) 


-  .  (C14b) 


/-l  y-l 


-SrS  vj 


Appendix  D 

SIMULATION  DESCRIPTION  AND  SOURCE-CODE  LISTINGS 


The  simulation  was  written  in  FORTRAN  and  run  on  a  PRIME  computer  system.  A  flow  chart  is 
presented  in  Fig.  Dl.  Parameters  are  input  (1)  directly  by  the  user  and  (2)  via  data  files.  An  example 
of  user-interactive  data  is  shown,  in  the  proper  sequence,  in  Table  Dl.  The  first  input  data-file  defines 
the  ship-motion  spectra  and  the  second  defines  the  US  parameters;  example  files  appear  in  Figs.  D2 
and  D3.  The  ship-motion  input-file  is  called  SHIP-PARMS  and  the  US-parameters  input-file  is  called 
SCATT-PARMS.  Data  are  output  in  binary  format,  one  record  per  time  increment.  The  output  file¬ 
names  are  defined  in  a  DATA  statement  in  subroutine  INIT.  We  note  that  the  compiler  replaces  the 
PRIME  statement 

SINSERT  filename 

with  the  contents  of  the  file  that  is  named  'filename'. 


Fig.  Dl  —  Simulation  flow  chart 
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Table  D1  -  The  User 
Input-Sequence 


0.0300 

Radar  Wave  Length 

0.0 

Pet.  Range  Error 

0.002 

Step  Size 

0.135 

Wave  RMS 

Yes 

Multipath  Indicator 

10.00 

Run  Time 

10000,7,0 

Initial  Ship  Position 

0.0,3.19,0 

Initial  Ship  Speeds 

88.0 

Initial  Ship  Orientation 

0.0070 

Ship  Turn  Rate 

SCATTERER  DEFINITIONS 


LOCATION 

(WRT  CG)  RCS 

SCATT.  PLATE 

X 

Y 

Z  AMP 

TYPE  1 

DIAM 

-50. 

-4. 

2.  4. 

3 

0. 

-5. 

-4. 

7.  0. 

2 

5. 

3. 

0. 

20.  3. 

3 

0. 

10. 

-4. 

7.  0. 

2 

5. 

20. 

-4. 

7.  0. 

2 

3. 

30. 

-3. 

15.  5. 

3 

0. 

Fig.  D2  —  A  £/S -parameter  input-file 

CG-PROCESSES  SPECTRUM-DEFINITION 

FREQ  (RAD/S)  BW  (RAD/S) 

RMS  VALUE 

PROCESS 

0.386 

0.114 

0.344 

X 

0.477 

0.182 

0.357 

Y 

0.432 

0.159 

0.411 

Z 

0.523 

0.136 

0.0539 

PHI 

0.409 

0.114 

0.0089 

THETA 

0.455 

0.159 

0.0066 

PSI 

Fig.  D3  —  A  ship-motion  input-file 


The  remainder  of  this  appendix  lists  the  insert  files  followed  by  the  source-code  except  for  that 
used  to  simulate  multipath.  The  multipath  source  code  is  presented  in  Ref.  98. 
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SINSFPT  FILES  (PRIME  FOETPAN) 


SINSFRT  SIZE 

IHTFGFP  DIMES ,DIMEN2 
PA  RAMETER  (DIMFN=25> 

PARAMETER  (DIMFN2=DIMEN*2> 

SINSERT  BLOCK  1 

INTFGER  **»  ISEED 

COMMON /BLOCK  1 /I SFFD.DT, TIME 

SINSERT  BLOCK 2 

COV¥ON/BLOCK2/XYZ(3,niMFN)»AMP(OTHFN),T(3»6,DIMEN), 

*  CNTFRC(« ),PNDWTH(6),PMS(6)  , RNGERR ,VAVRHS,RADHGT 


FILE:  MAIN 


CALL  PPIVFP 
CALL  EXIT 
FNP 


FILE:  DRIVER 


FURFOSE: 

TO  GENERATE  A  TIME  SERIFS  OF  RAPA. R  SIGNALS  RESULTING  FROM 
ILLUMINATION  OF  A  SHIP  AT  LONG  RANGE.  THIS  ROUTINE  IS  THE 
DRIVER  FCR  THF  PACKAGE. 

IN  FUT : 

ARGUMENTS: 

NONF, 

OUTPUT: 

ARGUMENTS: 

NONF. 

AUTHOR:  P.Y.NOPTHAM. 

DEVELOPED:  <»/82. 

MODIFIED:  1982,  1983. 


SUBROUTINE  DRIVFP 
INSFRT  S YSCOM>A$KEYS 
INSFPT  RCS>ROUTINpS>SI7p 
INSFRT  RCS>ROUTINES>nLOCK  1 
INSERT  RCS>ROUTINES>RLOCK2 

DOUBLE  PRFCISION  TMP ,PF LRI { 3 , DIMEN  ) 

DOUPI.F  PPFCISION  DV,DY  ,P7,PR 

LOGICAL  UPDATE (DIMES) ,FTPST(DIMEN) ,MFFLAG 

INTEGEP  SCTTYP(DIMpN) 

PARAMETER  ( PION2= 1 . 57 1 , TVOF 1=6 . 2 83 2 ) 

PI" ENS ION  LUNITS(9),XSpIF(6)  ,XCG(6) ,RT(3,PTMEN),PSI(pIMEN) 
*  DIAM(DIMEN) 

INITIALIZE  SIMULATION 


oo  on  no  no  o n  on  o  o 
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CM  I  INIT<LUNITS,NSCATT,SCTTvP,DIAM,VAVLEN,MFFt*.C) 

LUN  1  =LUNITS (  1  ) 

LUN2=LUNITS(2) 

LUB’*Ll*NITS(?) 

LUN4=LUNITS(4) 

LUN5=Ll’NITS<5) 

LUN6=LUNITS ( 6 ) 

IUN7=LUNITS(7) 

LUNP=LUNITS<8> 

LUN9=LUNITS<9) 

10  CALL  TNOUA ( *  ENTER  MAX  TIKE:  *,16) 

FEA  D { 1  ,*,EPF=10)  TMAX 
MAX=TMAX/DT+1 
DO  15  1=1, DIMEN 
FIRST(I)=.TPUE. 

UPDATE (I )= .TRUE. 

15  CONTINUE 

BEGIN  SIMULATION  LOOP 
DO  90  1=1, MAX 
?IME=DT*I 

UPDATE  6  DOF  POSITION  VALUES  FOR  SHIP'S  CG 
DO  2  0  J= 1  ,6 

CALL  PPOCES (J,CNTFPQ(J),PNDWTH(J),0, .FALSE » ,XCG (J),TI*'F) 
XCG(J)=RKS(J)*XCG(J) 

2C  CONTINUE 

UPDATE  SCATTEFFP  POSITIONS 
DO  35  K  =  1 ,NSCA TT 

CALL  LINHAP<?< 1, 1,K) ,3,«,XCG,RI( 1,*>) 

DC  30  .7=1,3 

RI(J,K)=RI(J,K)+XY7(J,X) 

30  CONTINUE 

35  CONTINUE 

UPDATE  SHTP  "ANEUVFPING  POSITION 
C*  LL  SHIP(TIMF ,XSHIP) 

UFDATF  PADAP-TO-SCATTFREP  RANGE 
DO  40  r= 1,NSCATT 

CALL  TRNSFF(FI(1,K),XSHIP,PADHGT,DBLPI(1,K)) 
PSI(K)=ASPECT(ySHIP( 1) ,XSKIP(2),XSHIP(6) ,K,FIPST(K))+XCG(6) 
4C  CONTTNUF 

UrCATE  PCS  ANT  GLINT  FPOCESSES 
PHI=XCG(U) 

CALL  e'IGNIS(PSI  ,  ISC  ATT  ,SCTT  YP  ,DI A  M  ,PH  I  ,WA  VLEN  ,DB  LPI ,  X  SHI®  , 

*  PCS,AZGLNT,ELGLNT ,MPFLAC) 

OUTPUT  DESIRED  DATA 
WRITS ( I UN7 )  PCS 
WRITE  (  LUK  P  )  f.  ZGLNT 
W  PITF (LUN  9 )  FI.CLNT 

C  TURN  OFF  INITIALISATION  FLAGS 
DO  PC  F=  1  ,NSCA.TT 
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I 


IF( .NCT.FIFST(K))  GO  TO  »0 
PIPST(K)=.EAISF. 
UFDATE(K)=. FALSE. 

60  CONTINUE 
90  CONTINUE 

r 

C  CICSF  OUTPUT  PII ES 
DO  100  1=1,9 

CALL  CL0S$A(LUNITSm-4) 

100  CONTINUE 
r 

C  END  OF  SIMULATION 
PETURN 
END 


C  FILE:  I  NIT 

C - - - - - 

C  PURPOSE: 

C  TO  INITIALIZE  THF  SIMULATION  PASANFTERS. 

C 

C  INFUT : 
r.  ARGUMENTS: 

C  NONE. 

C 

C  OUTPUT: 

C  ARGUMENTS: 

C  LUNITS  =LOCICA L  UNIT  NUMBERS  FOR  THE  OUTPUT  FILES. 

C  NSCATT  =NUMBER  OF  SCATTFPERS . 

C  SCTTYF  =LNDEX  INDICATING  TYPF  OF  SCATTERER. 

C  =1  FOR  CONSTANT  TYPE. 

C  =2  FOP  FLAT-FI.ATF  TYPF. 

C  =?  FOR  CORNER  TYPF. 

C  DIAM  =DIAMETERS  OF  THE  FLAT  PLATF  SCATTEREPS. 

C  WAVLFK  =PATA R  HA VF  LENGTH. 

C  MPFLAC  ="ULTI A  PTH  FLAG. 

C  =  .TPUF .  IF  MULTIPATH  EFFECTS  ARE  TO  BE  ACCOUNTED  FOR 

C  =.FALSF.  OTHEPHISF. 

C  /BL0CK1/: 

C  DT  =STMULATION  STFP  SIZE. 

C  ISEED  =P ANDOM  DUMBER  GENERATOR  SEEDS  (1*4). 

C  /PL0CK2/: 


RADHGT  =HEIGHT  OF  THF  RADAR  ABOVE  THF  SEA. 

CNTFPO  =CFNTER  PRFCUENCIPS  OF  THE  SHIP-MOTION  SPECTPA. 

BNDWTH  =P ANDHITTH  OF  THF  SHIP-MOTION  SPECTRA. 

RMS  =RMS  VALUES  OF  THF  SHIP  MOTION  PROCESSES. 

PNGEPP  =P ANGF  FPROP  IN  THF  PHASE  CENTER  LOCATION. 

XY?  =CAFTFSIAN  COORDINATES  OF  THF  SCATTERING  CENTERS 

IN  THE  SHIP-CG  AXIS  SYSTEM. 

AMP  =ANPLITUDES  OF  THE  SCATTERERS. 

T  =TRANSFCRHATION  FFOM  SKIP-CG  6-DOF  POSITION  TC 

SCATTEFFR  f.-rOF  POSITION. 


AUTHOR:  D.Y.NCPTHAN. 
DEVELOPED:  U/P2. 
MCDIEIFD:  1962. 
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SUBROUTINE  I N IT< LUK ITS , NSC ATT , SCTTYP ,DIAN , NA VLFN , MFFLAG ) 
SINSEPT  S YSCCW>AS  KFYS 
SINSFRT  RCS>POUTINFS>SIZE 
SINSFRT  FCS>POUTINES>PLOCK1 
SI NSFPT  ECS>FOUTINES>BIOCK2 

INTECEP  LUKITS(1),FNAf*(8,<)) 

INTEGER  SCTTYP(NSCATT) 

LOGICAL  FUNIT,YESNO,RPFLAG 
INTEGER  *4  11,12 
PA  PAHFTEP  (TWCPI  =  6  .2832) 

TATA  FNAN/'X-CG* ,6** 

*  * Y -CG  *  ,  6  *  '  *, 

*  '  Z-CO  ’  ,  6 *  *  ’, 

*  ' FHI -CG ' , 5*  * 

*  'THETA-CG' ,4**  •, 

*  •PSI-CG',5*'  ’, 

*  ' RCSHAG '  ,  5*  ' 

*  * AZ-OLTNT'  ,4*  •  *, 

*  *  FL-GLINT  *  ,4*  •  •/ 

C 

ISEED=  1 
F ADHGT=40 • 0 

?  CALL  TNOUA ( ' FNTEP  RADA  B  HA VLENG^H  :  ’,23) 

READ ( 1 ,  *  ,FFF  =  3 )  WAVLEN 

CALL  GFTSCT (XYZ,NSCATT,AEP, SCTTYP,DI AN) 

DO  5  I»1,NSCATT 

IF(SCTTYP{I).EC.3)  A*P(I)=AKP(I)/VAVLEN 
5  CONTINUE 

CALL  GETPRF  (CNTFRO  , PNDWTH ,RNF) 

CAL!  GPTPAP(XYZ,NSCATT,T) 

7  C AIL  TNOUA ( * FNTFP  FERCFNT  PANGF  FRROR :  ’,27) 

RFADO  ,*,ERR  =  7)  FC7 
PNGEFF=WAVLEN*(PCT/100.) 

10  CALL  TNCUA ( ' * NTFP  FTP  SIZF:  ’,17) 

PEAD( 1  ,*,EFR=1C  )  DT 

OPEN  OUTPUT  FILES 
DO  45  1=1,9 

I F ( FUNIT  <  NFU  )  )  GO  TO  40 
STOP 

40  CALL  OPE NS A ( A SWRIT+A $SAWF,FNAR(1,I),16,NFU) 

LUNITSd  )=NFU+4 

45  CONTINl'p 

C 

C  SFT  UP  FOR  NULTIPATH 

46  CAII  TNOUA ( * FNTFP  WAVF  RMS :  *,16) 

REA D ( 1 , *  ,EPR  =  46  )  WAVR*S 

IF  (  .NOT.YESNO!  ‘USE  EULTIPA.TH?  ’))  CO  TO  55 
t!PFLAG  =  .TpUF. 

T  1  =  10  1 
12=1001 

DO  50  1  =  1 , NSC  ATT 

CALL  NFINIT(L,C.0,WAVP!T,WAVLFN,0,L1,T2,DT) 

1 1  =  1 1*5 
T2=I2*6 
CCNTINUP 
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55  MPFLAG=.FALSF. 

60  CONTINUE 

r» 

C  CALCULATE  AND  STOFE  SPECTFA 
DO  70  1=1,6 

CALL  PROCES(I,CNTFPQ(I),ENDHTH(I),INTL(I) , .TRUE. ,P«DVAR,0 .0 > 
70  CONTINUE 
BFTUPV 
END 


C 

c 

c 

c 

c 

c 

c 

c 

r 

C 

C 

c 

c 

c 

r 

C 

c 

c 

c 

c 


SU5PCUTINF  GETSCT(XYZ,NSCATT,AHP,SCTTTP,DIA«> 

*1 NSFPT  SYSCOM>ASKEYS 

INTFGEP  FN AME(  16  ) 

LOGICAL  FUNIT 
INTFGEP  SCTTYF(NSCATT) 

DIMENSION  XY7 ( 3 ,NSCATT ), AMP ( NSCATT) »DI AH (NSCATT) 

DATA  FNAME/'SCATT-PARHS  '.lO**  '/ 

C 

C  OPFN  I N FUT  FILFS  AND  PEAD  IN  SCATTFPING  PARAMETERS 
IF  <  FUNIT (NFU ) )  GC  TO  10 
STOP 

10  CALL  CPEMS*(»$°FAD+*?SAMF,FXAMF,32,NEU) 

Ll’N  =  NFl'  +  U 

PF  AD (LUN  ,11)  DUM 1 ,DUH2 ,DUM3 ,DUN4 

11  FCFMAT(A4/A4/A4/A4) 
t'SCATT  =  0 

DO  15  1=1,100 

PEAD (LUN ,*,FND  =  20,E??  =  20)  XY Z ( 1 , 1 ) , X YZ (2 , I ) , X YZ ( 3 , 1 ) , AKP< I ) 
*  SCTTYP<I),DIAM<I) 

NSCATT=NSr ATT+1 
1C  CONTTNU" 

20  CCNTINUF 


FILE:  GETSCT 


PURPOSE: 

TO  PEAD  IN  THE  SCATTFPING  PARAMETERS  FOP  THF  INDIVIDUAL 
UNIT-SCATTEPEPS. 


INPUT: 

ARGUMENTS: 

NONE. 


OUTPUT: 

ARGUMENTS : 
XYZ(3,I) 

NSCATT 

AMP 

SCTTYP 

DIAM 


=CAPTESIAN  COOPDS.  OF  SCATTERFR  LOCATIONS 
IN  THF  SHIP-CG  AXIS  SYSTEM. 

=NUH3FR  OF  SCATTFPFRS. 

=AHPLITUDF  OF  THE  SCATTFPEPS. 

=  1 NDIC ATES  TYPE  OF  SCATTFRER. 

=FLAT  PLATE  DIAMETER,  IF  APPROPRIATE. 


AUTHOR:  D  .  Y . NORTHAM . 
DEVFLOPED :  4/82. 
MODIFIED:  1982. 


**  o  :■>  o  n 


NRL  REPORT  8IS7 


C  WPITF  "ESSAGE  AND  CIOSR  INPUT  FILE 
WRITE (1,25)  NSCATT 

25  FORMAT (  1  THESE  APE', 13,*  SCATTERERS. * ) 
CALL  CLOSSA(NFU) 

FFTUPN 

END 


C  FILE:  GETFRM 

C - 

C  PURPOSE: 

C  TO  READ  IN  THF  SHIP-MOTION-PROCFSS  PARAMETERS. 

C 

C  INPUT: 

C  ARGUMENTS: 

C  NONE. 

C 

C  OUTPUT: 

C  ARGUMENTS : 

C  CNTFRO  =CENTEP  FREQUENCIES  OF  THE  SHIP-MOTION  SPECTRA. 

C  BNDHTH  =E ANOWIDTH  OF  THF  SHIP-MOTION  SPECTRA. 

C  RMS  =RMS  VALUES  OF  THE  SHIP  "OTION  PROCESSES. 

C 

C  AUTHOR:  D.Y.NOPTHA". 

C  DEVELOPED:  4/82. 

C  MODIFIED:  1982. 

C - 

SUBROUTINE  GETPRM < CNTFPQ ,BNDVTH , RMS ) 

SINSFRT  SYSCOM>A$KEYS 

INTEGER  FNAMFOO 
LOGICAL  FUNIT 

DIMENSION  CNTFRQ(6),BNPWTH(6),RMS(6) 

DATA  FNAME/' SHIP-FARMS' ,11*’  V 

r 

C  OPEN  INPUT  FILE  A"D  RFAD  IN  MOTION  PAR  AMFTE°S 
IF(FUNIT(NFU)  )  C-C  TO  10 
STOP 

10  CALL  OP*NSA (ASREAD+ASSAMF, FRAME, 32, NFU) 

I  UN=NFl'  +  4 

REA.D(LUN  ,13)  DUM1  ,DUM2  ,DUM3 
13  FOR MAT ( A4/A4/A4 ) 

NFA  RMS  =  0 
TO  15  1=1,6 

NPARMS=NFARMS+1 

READdUN,*)  CNTFPQU ) ,PNDWTK(I> ,RM?(I) 

15  CONTINUE 

CLCSF  INPUT  FILE 
20  CALL  CLCSS A  < NEU ) 

F.ETURN 
FND 


FILF:  GFTMAP 


FUFPOSr: 

TO  DEFINE  TUE  LINEAR  MOTION -MAFS  FOR  THE  INMVIDUAJ 


C 
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UNIT-SCATTEPEF  LOCATIONS. 


INPUT: 

ARGUMENTS: 
XYZO ,  I ) 

NSC ATT 

OUTPUT: 

ARGUMENTS: 

T 


'CAPTFSIAH  COORDS.  OF  SCATTERER  LOCATIONS 
IN  THr  SHIF-CG  AXIS  SYSTFN. 

•  NUHPFF.  OF  SOATTFPFRS. 


=TPANSFOPNATION  FPCK  SHIP-CG  f -DOF  POSITION  TO 
SCATTFRER  6-DOF  FOSITIO*. 


AUTHOR:  P .  Y .  NORTHA  N  . 
DEVELOPED:  4/82. 
MODIFIED :  1982. 


SUPPCUTINP  GFT¥AP(YYZ,N?CATT,T) 
DIMENSION  XYZO,NSCATT),TO,f  ,NSCATT) 

DO  20  K  =  1 , NSCATT 
DO  10  1=1,3 
DC  10  J=1 ,3 
T(I,J,N)=0.0 
IF(I.EQ.J)  T<I,J, 10*1.0 
10  CONTINUE 

T ( 1 , 4  ,K )  =  0 .0 
T(2,4,K)=-XYZ(3,K) 

T(3,4,K)=XY?(2,K) 

T(1,5,F)=XY”(3,K) 

T(2,c,K)=C.O 
T(3,5,K)*-XY7( 1,N) 
T(1,6,K)=-XYZ(2,K) 

T(2,6,K)=XY7(1,X) 

TO, 6, 10=0. C 
20  CONTINUE 
RETURN 
FND 


FILE:  TRNSFR 


PURPOSE: 

TO  OPTAIN  THE  COORDINATES  OF  THE  ITH  SCATTERER  IN 
THE  RADAR  AXIS-SYSTEF. 

INFUT: 

APGU'TNTS: 

PI  =SCATTEPFR  COORFIN ATFE  IN  THE  SKIP-CG  AXIS  SYSTFN. 

XSHIP  =SHIP-CG,  8-DOF  POSITION  DUF  TO  MANEUVERING. 

RATHGT  =HFIGHT  OF  THF  RADAR  A°OVF  THF  SEA. 

OUTPUT: 

A  RGU“FNTS : 

DBLRI  =F ADA.R  AXIS-SYSTEH  COOPDS.  OF  THF  ITH  SCATTERFP  (1*4), 
AUTHOR:  D.Y. NORTHA* . 
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C  DFVELOPFD :  1982. 

C  MOTIF  TrD :  1982. 

SUBROUTINE  TPNSF P  (  P I  , X SH IP ,  P APHC,", TBLRI ) 
COUPLE  PRECISION  DPLPI 

DIM  F  NS  TON  PI(3),ycL'Tp(6  )  ,  A < 3 , 3 >  ,DPLP I ( 3 ) , Y ( 3 ) 

SPHI=SIN(XSHIP(4 )) 

CFHI=COS(XSHIP(H)) 

STHr=?IN(XSHTP(5) ) 

CTFF=COS(XSHIP(S) ) 

SPSI=SIN (XSHIP(6 ) ) 

CPSI=C0?(XSHIP(6 ) ) 

A  (  1  , 1  )=CPSI*CTHF 
A  ( 1  ,2)=CPSI*STFF*SCHI-SFSI*CPHI 
A (  1 ,3)=CPSI*STKE*CPHI+SPSI*SPHI 
A(2,1)=SPSI*CTFF 
A(2,2)=SFSI*STHF*SFHI+OPSI*CPHT 
A(2,3)=SPSI*STHF*CPHI-CPSI*FPHI 
A  (  3 , 1 )  = -STHE 
A(3,2)=CTHF*SPUI 
A(3,3)=CTHE*CFHI 

POSITION  CUE  TO  SFA  MOTION 
TO  20  1=1,3 
y  ( i ) =0 . o 
DO  1C  J  =  1 

Y(I)=Y(I)+A(I,J)*PI(J) 

1C  CONTINUF 
20  CONTINUF 

POSITION  TUF  TO  K  »  N  r U  VF.PI  NO 
TO  30  T=  1 , 3 

D0LFI(I)=Y(I)+XSHIP(I) 

30  CONTINUF 

CORprcT  FCP  PATAF  BEING  ABOVE  THE  rA PTH 
DFLFI(3)=DPI?I(3)-PACFCT 

PFTUPN 
FN  0 


FIIF:  LINMAF 


PURPOSE: 

TC  I,INEBPLY  TP  A  N  SFOc  K  A  V  INPUT  VFCTCP  INTO  AN  OUTPUT  VECTOR 


I N  E  U  ^  : 

APGl’MFNTS: 

A 

N'OKA 
NO  OLA 
X 


=  TP.*vSFGPy*‘FION  K  ' TDI v . 

=  N 'J v 0 F F  Or  POKS  IN  ’  A  •  . 
=NUMnrP  Or  COLUMNS  IN  ’A’. 
=  V  EOTOP  TO  ft  tfa:jcfc°mft. 


OUTPUT  : 
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C  ARGUMENTS: 

C  Y  =FESITLTANT  VECTOR . 

C 

C  AUTHOR:  D. Y.NCPTHA*. 

C  DEVELOPED:  4/62. 

T  MODIFIED:  1982. 

c - - - 

SUBROUTINE  LINE AP ( A ,NRCW A , NCOL A , X , Y ) 

C  DIMENSION  A(NROVA,NCOLA),X(NCOLA),Y{NROWA) 

DI“FNSICN  A(3,6U(6),Y(3) 

C 

DO  20  I=1,NF0WA 
Y(I)=0.0 
DO  10  J=  1  ,NCOLA 

v(I)=Y(I)  +  A(I,iT)*X(J) 

10  CONTINUE 
20  CONTINUE 
RETURN 
END 


C  FILE:  SIGNLS 

C - - - 

C  PUPPOSF : 

C  TO  CALCULATE  THE  RECEIVED  RA D AR -SIGN ALS . 

C 

C  INPUT: 

C  ARGUMENTS: 

C  PSI 

C  NSCATT 

C  SCTTYP 

C 
C 

C  DI  AM 

C  PHI 

0  VAVIEN 

C  PPLRI 

C  XSHIP 

C  KPFLAG 

C 

c 

C  AMP 

C 

C  OUTPUT: 

C  A  PGUMENTS : 

C  FCS 

C  AZGLNT 

C  ELGLNT 

C 

C  AUTHOR:  D. Y.NCFTHAM. 

C  DEVELOPED:  4/82. 

C  MODIFTFD :  1982,  198?. 

C . . - - - - - 

SUBROUTINE  SIGNLS (FFI , NSCATT, SC’rTYF,DI AH, PHI, WAV LEN, DP LPI,X SKIP, 
*  pCS,A7GLNT, ELGLNT, MPFLAG) 

SINSFPT  PCS>P0UTTNFS>SI7F 
SINSERT  FCS>F0UTINES>RL0C’f1 


“SHIP  * S  ASFEC  ANGLE. 

“NUMBER  OF  UNIT-SCATTFRERS 

“FLAG  INDICATING  THE  TYPE  CF  SCATTERED 

=0  FOP  CONSTANT. 

=1  FOP  FLAT  PLATE. 

“DIAMETER  OF  THF  ELAT  PLATES. 

“POLL  ANGLF  OF  THE  SHIP. 

“RADAP  WAVELFNGTH. 

“POSITION  VFCTCR  CF  THE  ITH  SCATTERER,  IN  EARTH  COORDS. 
“POSITION  VECTOE  CF  THE  SHI*,  IN  EAPTH  COORDINATES. 
“MULTIPATH  FLAG. 

=.TF"r.  IF  MULTIPATH  EFFECT*  ARE  TO  BE  CALCULATED. 
“.FALSE.  OTHERWISE. 

“AMPLITUDE  OF  THE  ITH  SCATTFRER. 


“TA.FGFT  PCS 

“TAPGFT  AZIMUTH-GLINT 

“TAPGET  EIEVATION -GLINT 
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SINSERT  RCS>P0UTINFS>°10CK2 

TCUBLE  FRECISICN  r P  ,RA TIC , F ANGF  , DX , P * , P » 

DOUBLE  PRECISION  PnI.RI 
INTEGER  *U  ISFEP 
INTFGFP  SCTT VP ( DIMFN  ) 

LOGICAL  MPFIAG 

DIMENSION  RANGF<DIMEN2),DRL?I< 3 , PI HFN ) , A ( D INEN7  )  , XSHIP < 6 ) , 

*  PSI(PIKFN),DIAM<PIHRN),PHASE(DINEI»2>,GMI>RNe<DIPEI»> 
PARAMETER  ( TKCP 1=6 . 763 2 ) 

PATA  ISF.ED/10  .0/ 

CALCULATE  RADAR-TO-SCATTFRER  RANGES 
WAVRUK=TWOPT/WAVLFN 

FSHIP=SQRT(XSHIP< 1 )  * *2*XSHI F < 2 ) * *2* X SKIP< 3 > »*2 > 

PO  10  1= 1 , NSC ATT 
PX=DELPI( 1,1) 

P Y=DBLPT (2,1) 

P7=DBLRI (3,1) 

DW=RNGERR*RANP?A(ISEFP) 

GNPRNGd  )=DSQRT(PX**2  +  DY**2)-»DW 
RANGrd)=rS0RT(DX**2*DY**2+D?»*2>*DW 
IF (HPFLAG )  RANGE (I+NSCATT )  = 

*  PSORT ( DX*  *2*DY*  *2*(XY7 (3,1 )*XSHIP( 3)-FADHGT)**2)*DH 
10  CONTINUE 

SUN  1=0.0 
SUM2=0 .0 
SUP  3 =0.0 

IF ( MPFLAG )  NPTS=2*NSCATT 
IF( .NOT.MPplAG)  NPTS=NSCATT 

CALCULATE  SCATTERING  AMPLITUDES 
DO  15  J=1 .NSCATT 
A (J)=AMP(J) 

IF(SCTTYP(J).FC.2)  CALL  FLTPLT ( PI  AM ( J ) ,W A VNUH , PSI ( J > ,PHI , A < J > > 
IF ( .NOT. MPFLAG  )  GO  TO  13 


CALCULATE  MULTIPATH  AMPLITUDE  AND  PHASF  EFFECTS 
TAPHGT=VSHIP (3 )+XY7 ( 3 ,J ) 

CALL  MPMAIN (J.TARHGT  ,R A DHGT , GNPPNG ( J) , TINE# PSISPC, XRF AL » X IRAG ) 
F 1  =SQRT ( ( 1 .0  +  XRFAL ) * *2+ XIM AG* *  2 ) 

F2=F1*SCPT(XRFAL**24XIMA.G**2) 

A (J+NSCATT)=F2*A (J) 

A(J)=F1*A (J) 

PHASE( J)=ATAN2(XIMAG,1.0+XRFAL) 

PH »SF( J+NSCATT )=ATAN2(y TMAG*( 1.0*2. 0*XPEAL),XRF AL+XSEAL**2- 

*  XIMAG**?) 

GO  TO  15 

13  PHASF ( J ) =0 . 0 
If  CONTINt!p 

CALCULATE  FCS  AND  CLIFT 
DO  25  1=1 , KPTS 

TP(I.IE.WSCAT"')  T  ,=I 
IF(I.GT.NSCATT)  H=I-KSCATT 
AZCCEF=(XSHIP(2)*P3LFI<1,T1))/PANGE(I) 

*  -<XSHIP(1 )*PBLFI(2,I1 >)/RANGE(I) 
rc  25  J= 1  . NpTS 


»  *  .  ’j  i 


-.-.S  .\*N>- 
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IF( J.LE.NSCATT)  C-0  TC  21 
J  1=J-N*>C  R  TT 

P  ZI  =  - ( XSHIP ( 3 )+PPLFI ( 3 , J )+R ADHGT ) 

GO  TO  22 

21  J1=J 
9ZI=PBLPI(3, J) 

22  CONTINUE 

ELC0EF=RZI*(DBLRI(1/J1)*DBLPI(1#I1)+DBL8I(2#J1 )* 
*  DBLRI(2,I1 ) )/(RANGE(I)*(RANGE(J)**2)) 

CALCULATE  PHASF  DUE  TO  RANGF 
DR=RANGE(I )-RANGF< J) 

RATIO=DR/WAVLEN 

ANGIE=2.0*TW0PI* (RATIO-DINT (PATIO) ) 

n 

C  ADD  MULTIPATH  PHASE 

ANGLE=?NGLE+FHASE(I)-PHASE(J) 

n 

C  C » ICUT  ATE  SUMMATION  T^RFS 

TEPM=A(I )*A(J)*COS< ANGLE) 

SUF1=SUNH-TFRM*ELCOEF 
SUM2=SUF2+TERM*A7CCEF 
SUM3=SUM3+TERH 
25  CONTINUE 

SUM3=SUN3+1 .OF-IO 
RCS=SUF  3 

r 

C  CALCULATE  GLINT  IN  NETERS 
AZGLNT=SUH2/SUF3 

ELOLNT=PSHIP*SUH1/FUM3  -  (XSHIF ( 3 ) -R ADHGT) 

RETURN 

END 


C  FILE:  FLTPLT 


C  FVPPOSF: 

C  TO  CALCULATE  A  FLAT -PL A TF -TYPE  AHPLITUDE  AT 
C  NEAR-PEPPENDICL’LAP  INCIDENCE. 

r 

C  INPUT: 

C  ARGUMENTS; 

C  DIAM  =DIAMFTEES  OF  THE  PLATES. 

C  WAVNUM  =HA VF  SUFFER  FOP  THE  PADAR. 

C  AZIMTF  =PACAR-TO-PLATE  AZIFUTK  ANGLE. 

C  FLFV  =PADAR-TO-PLATE  ELEVATION  ANGLE. 

C 

C  OUTPUT : 

C  APGUVFNTS : 

C  AFP  =FESULTING  SCATTERED  AFPLITUPF. 

C 

C  AUTHOR:  D.Y  .NORTH*  F . 
r  DEVELOPED:  U/R2. 

C  MODI  FI FD :  1982. 


no  no 
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C1  =  HAVVl,H*(CTA?t**2)/‘?0PTFI 

CALCUI ATE  AZIMUTH  COMPONENT 
IF(AZIMTH.NE.O.O)  GO  TC  15 
A?MAMF=1.0 
GO  TC  20 

15  C2=VAVNUH*DIAH*SIN<AZIMTH-PI0N2> 
AZMAMP=ABS(STN(C2)/C2) 

CALCUIATE  FLEVATION  COMPONENT 
20  IF(ELEV.NE.O.O)  GO  TO  25 
ELEAMP=1 .0 
GO  TC  30 

25  C3*HAVNUM*DIAK*SIN(ELEV) 

ELE AMP=ABS (SIN(C3)/C3) 

CALCUT ATE  TOTAL  AMPLITUDE 
30  AMP=C1*SQRT(ELEAMP*AZMAMF) 

RETURN 

END 


C  FILE:  SHIF 

C - 

C  PURPOSES 

C  TC  CALCULATE  THE  EANFUVEPING  TP AJFCTORT  CE  THE  SHIP 
C  IN  EARTH  COORDINATES. 

C 

C  INFUTs 
r  ARGUMENTS s 

C  TIME  ^CURRENT  SIMULATION  TIME. 

n 

C  OUTPUTS 
C  ARGUMENTS  5 

C  YSHIP  =EARTH  COORDINATES  CF  THE  SHIP  DUE  TO  MANEUVERING • 

C 

C  AUTHORS  D.Y.NORTHAM. 

C  DEVELOPED  8  1982. 

C  MCCIFIECs  1982,  1983. 

C - 

SUBROUTINE  SHIF(TIMF,XSHIP) 

LOGICAL  FIRST 
DIMENSION  XSH IP ( 6  ) 

DATA  XO,  YO,  Z0,VX,VI,VZ/6  *0.0/,  FIRST/.  T»>UE. / , CKEGA  ,  YSI0/2*0 . 0/ 

C 

IF ( . NOT .FIPST )  GO  TO  10 
FIRST®. FAt'E. 

5  CALL  TNOUA( ' ENTER  INITIAL  SHIF  CCOPDSs  *,27) 

PFAD ( 1  ,* ,EPP*5)  XO ,Y0, ZO 

7  CALL  TNOUA ( ' rNTEP  SHIP  SFEFDs  ’,18) 

PFAD( 1 ,  *  ,  FPR=7 )  V 

8  CALL  TNOUA ( * FNTFE  INITIAL  SHIP  POTATION  ANGLE  (DEG)s  ’,41) 
READ ( 1 , * , FPP=8 )  FSIO 

PST0=FSI0/57.3 
VX=V  *CCS ( FSIO ) 

VY  =V  *SIN ( PSIO ) 
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Q  CAIL  TNCUA(  '  ENTER  SHIP  TURNING  PATE  (RAP) 
PEAD(1,*,ERR=9)  ONEGA 
10  XSHIP( 1 )=X0+VX*TIWE 
XSHTF(2)=Y0*VY*TINE 
XSHIF(3)  =  Z0-fVZ*TTNF 
XSHI P ( 4  )=0  •  0 
XSH IP ( 5 ) =0 . 0 

XSHIF(6 )=PSI0+CHEGA*TI*F 

RETURN 

END 


FILE:  ASPECT 


FUPPOSE: 

TO  CALCULATE  A  SCATTERFP ' S  ASPECT  ANGLE  RFLATIVE  TO  THF  RADAR 

INPUT: 

AFGUHENTS: 

X  =X-COO?DINATE  OF  THF  SCATTEREP. 

Y  =Y-COORPIN ATE  OF  THE  SCATTEREP. 

PSI  =SH TP  ROTATION  ANGLE. 

SCTNUH  =INDEX  TO  THE  SCATTFRFPS  (I). 

FIRST  =FIRST-PASS  FI  AG . 

=.TPUF.  FOR  THE  FIRST  PASS  THROUGH  THE  SIHULATION. 

=  . FALSE.  OTHERWISE. 

OUTPUT: 

AFGUWENTS: 

ASPFCT  =SCATTEREF 'S  ASPECT  ANGLE  (RAD). 

AUTHOR:  I*.  Y.NOPTHAH. 

DFVELCPED:  4/82. 

WCCIFIED :  1982. 


FUNCTION  ASPECT(X.Y,PSI, SCTNUH, FIR  ST) 

INSERT  PCS>ROUTINES>SIZE 
LOGICAL  FIRST 
INTEGER*4  ISFED 
INTEGER  SCTNUH 
PIHFNSION  ASPCTO (DIHEN ) 

PA RAWFTR  (  CFGPFE=2 .0) 

PAT  A  ASPCT0/EINEN*0 .0/,ISEED/1/ 

ASSIGN  RAN DON  INITIAL  ASPFCT 
IF(  .NOT. FIRST)  GO  TO  10 

ASPCTO  (SCTNUH  )«2.0*(  RANDS  AdSFED ) -0  .  S  )*  (DFGREE/57 .296  ) 

CAICUI ATE  ASPECT 

10  A SFECT=ASPCTO (SCTNUH )*FSI+ATAN2(Y,X) 

20  PFTURN 
END 


FILE:  PROCES 


on  no  o  n  on 


NHL  REPORT  ttt7 


TO  REALIZE  THE  SHIP-CG  MOTION  PROCESSES. 


INPUT: 

ARGUMENTS ! 
II 

CNTFRQ 

BNDNTH 

ISFED 

FIRST 


OUTPUT: 

ARGUMENTS: 

RNDVAR 


'INDEX  TO  THE  SHIP-CG  PROCESSES. 

=CENTF5  FPPCUFNCIFS  OF  THP  SPFCTPA  (RA P/S). 

'P  A  NDWIDTFS  OF  THE  SPFCTPA  (R»C/S). 

'RANDOM  NUN PEP  GENERATOR  SEEDS  (1*4). 

'FTRST-FASS  FLAG. 

=.TRUp.  FOR  THF  FIRST  PASS  THROUGH  THE  SIMULATION, 
'.FALSE.  CTHEFHISE . 

'CURRENT  SIMULATION  TINE  <S>. 


*CUPRFNT  VALUE  CF  THE  DESIRED  PROCESS. 


AUTHOR:  D.Y.NORTHAN. 
DEVELOPED:  4/P2. 
MODIFIED :  1982. 


SUBRCUTINF  PRCCES( II, CFTFRO # BN DWTH,ISEED, F IRST .RNDVAR, TIHF) 
PAPAHETFP  (THOPI=6.2fl3?,PHSFAC*0.4449,NAXFPQ*31> 

LOGICAL  FIRST 
INTFGFP*4  I FEED 

CO**’ON/RVPLr/S(NAXFF0.6),H(NAXFR0,8),DH(6),PHI(NAXFR0,6), 

*  PARAMS(10),NFPEQ(6> 

IF(.NOT.FTRST)  GC  TO  100 

INITIALIZE  RANCOR  PHASES 
DO  5  1*1 ,NA  XFRC 

PHI(I,II)=TMOPI*RAND$A( ISEED) 

*  (~ORTTNUE 

INITIALIZE  PARAHFTEPS 
PARA  ESC 1)=CNTFR0 
PARAMS(2)=PNPWTH 
WNIN=CNTFRC-PNrWTH 
PARAMS(3)=HMIN 
NFREC(II)*EAXFRC 
DW (IT )=ENrVTH/10 . 

CALCULATE  THF  SFECTFA 

CALL  SPCTFM(1,PAFANS,NFREQ(II),DV(II) ,W( 1,II),S( 1,11)) 

STORE  NOPMALI7EP  SPECTNIIN 
MAX=NFPEC(II ) 

DO  20  I  =  1  ,  K  A  X 

S(T,II)*SC?T(2.0*S(I,II)*»'V:(TI)  ) 

20  CONTINUE 
rO  TO  200 

CALCULATE  THF  CURPENT  VALUE  OF  THE  PPOCFSS 
100  SU"=C.C 

KAX=NEPFC(II) 

DO  110  1*1,  NAX 
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CCSIN'E=COS(rHI(I,II)+W(T<TI)*TISE) 
SUM=SUM+S(I,II )*O0SINE 
110  CONTINUE 

RNrVAR=FMSFAC*SUM 
200  RETURN 
FND 


r.  FIIE :  SPCTPM 

c - 

C  PURPOSE! 

C  TO  GENERATE  THE  SFECTRUM  OF  A  PROCESS  OVER  A  RANGE  OF 

C  POSITIVF  FREQUFNOIFS. 

r 

C  INPUT: 

C  ARGUMENTS : 

C  TYPE  =FL  AG  INDICATING  TYPE  OF  SPECTRUM  TO  BE  GENERATED  (I). 

C  PAPAMS  ( 1 )  =CFNTEP  FRFQUENCY  OF  THE  (DNTRCDAL)  SPECTRUM  (RAD/S). 

C  PAPAMS (2)  =BANPWIDTH  OF  THE  SFECTRUM  (RAD/S). 

C  PAPAMS (3)  =LOWEST  FREQUENCY  AT  WHICH  THE  SPECTRUM  IS 

C  TO  PE  CALCULATED. 

C  NFREQ  =NUHPER  OF  FREQUENCIES  AT  WHICH  THE  SPECTRUM 

C  IS  TO  PE  CALCULATED. 

C  DW  =FREQUFNCY  INCREMENT  (EAD/S). 

C 

C  OUTPUT: 

C  ARGUMENTS: 

C  H  =FREQUENCIES  AT  WHICH  THE  SPFCTRUM  IS  CAICULATED  (RAD/S). 

C  S  =V A LUES  OF  THE  SPFCTRUM. 

C 

C  AUTHOR:  D . Y. NORTHAM. 

C  DEVELOPED:  U/82. 

C  MODIFIED:  1982. 

C - 

SUPPOUTINE  SPCTRM ( TYPE ,P ARAMS, NF PEQ ,DW , W ,S > 

INTEGEP  TYFE 

DIMENSION  PARAMS(1),S(1),V(1) 

C 

GO  TO  (  10)  ,  TYPE 
WPITE(1,1)  TYPE 

1  FCPMATCNO  SPECTRUM  PFFTNFD  FOR  TYPE*’, 12,’  (SPCTRM).’) 

GO  TO  20 

BANDPASS  SPECTRUM 
10  CNTFRC=PARAMS( 1) 

PNPWTH=PAFAWS<2> 

VEIN  =P  AP  AH S  ( 3  ) 

ALPHA=BNDVTH/2.0 
DP  15  1=1 , MFREP 

W(I)=WMIN+(I-1)*rW 

S(I )=ALPHA/(ALPHA**2+(W(I )-CMTFRQ)**2) 
S(I)=S(I)+ALPHA/(ALPHA**2+(W<I)+CNTFR0)**2) 

S ( I )=2. 0*S ( I ) 

15  CONTINUE 
20  RETURN 
END 
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